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ABSTRACT 


A  detailed  description  of  the  digital  computer  program 
developed  for  the  computation  and  plotting  of  Noise  Exposure 
Forecast  contours  is  presented.  All  routines  including  flow  • 
charts,  program  listings,  and  internal  data  organizations  are 
provided.  A  orief  summary  of  the  basic  Noise  Exposure  Forecast 
concepts  and  equations  aj»e  included.  Additional  information 
and  examples  of  the  results  of  the  computation  program  can  be 
found  in  the  following  reports  prepare^  in  performance  of 
Contract  FA68WA-1900:  FAA-NO-69-2,  FAA-NO-70-7 ,  FAA-NO-70-8, 
and  FAA-NO-70-9. 
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I .  INTRODUCTION 


This  report  presents  a  Ueta.tj.ed  description  of  a  digital 
computer  program  which  was  originally  developed  by  Bolt,  Beranek 
and  Newman,  Inc.  under  Federal  Aviation  Administration  Contract 
FA67WA-1705,  and  refined  under  Contract  FA68WA-1900,  as  an 
aid  in  producing  contours  of  equal  noise  exposure  using  procedures 
set  forth  in  Reference  1.  The  program  documentation  and  the 
preparation  of  this  report  was  conducted  in  performance  of 
Task  I  under  Phase  II  of  Contract  FA68WA-190Q.  The  material 
in  this  report  is  arranged  to  provide  both  the  user  and  the 
programmer  with  the  information  necessary  to  use,  or  modify 
the  routines.  It  is  helpful  if  the  reader  has  some  knowledge 
of  the  Noise  Exposure  Forecast  (NEF)  concept;  however,  this 
is  not  absolutely  necessary. 

Section  II  presents  a  users  guide  to  the  program  and  includes 
a  description  of  data  input  and  output  formats  and  the  organiza¬ 
tion  of  the  data  elements  for  program  execution.  An  alphabetic 
listing  of  error  messages  and  their  probable  causes  is  also 
presented. 

Section  III  outlines  the  general  computational  procedures 
for  determining  the  contour  locations. 

Appendix  A  provides  samples  of  the  input  data  and  associated 
program  output  for  the  NEF  contours  calculated  for  a  single 
flight  path. 

Appendix  B  provides  a  listing  of  the  mathmetical  relation¬ 
ships  pertinent  to  the  program. 

Appendix  C  provides  a  set  of  flow  diagrams  and  the  program 
listings  for  all  routines. 


II.  PROGRAM  USAGE 


A.  Discussion  of  Data  Required 

This  section  of  the  report  outlines  the  types  of  data  re¬ 
quired  In  order  to  compute  the  NEF  contours  for  each  flight 
path.  The  general  flow  of  control  within  the  program  and  the 
format  required  for  data  input  is  explained.  (A  complete  set 
of  flow  diagrams  and  program  listings  may  be  found  in  Appendix 
C.)  To  compute  the  NEF  contour  distances,  data  concerning 
aircraft  noise,  performance,  and  volume  of  activity  for  each 
flight  path  must  be  read  by  the  program. 

The  aircraft  noise  data  is  described  in  terms  of  the 
Effective  Perceived  Noise  Level  (EPNL).  An  EPNL  vs  Aircraft- 
to-Observer  Slant  Distance  function  is  specified  in  terms  of 
a  list  of  EPNL  values  at  various  slant  distances.  This 
function  varies  only  with  slant  distance  and  is  not  dependent 
upon  the  position  on  the  flight  path. 

The  aircraft  performance  data  is  described  in  terms  of  an 
Altitude  Profile  and  a  Delta-EPNL  Profile.  The  Altitude 
Profile  consists  of  a  list  of  X,  Y  values  (assumed  to  be 
connected  by  straight  line  segments)  which  specify  the  altitude 
of  the  aircraft  (Y-coordinate)  as  a  function  of  the  aircraft's 
distance  along  the  flight  track  (X-coordinate)  from  start  of 
takeoff  roil. 

The  Delta-EPNL  Profile  also  consists  of  a  list  of  X,  Y 
values.  This  profile  allows  the  user  to  modify  the  overall 
noise  level  characteristics  of  the  aircraft  as  a  function  of 
distance  along  the  flight  track.  This  capability  is  useful 
for  introducing  power  cutbacks  at  selected  points  along  the 
flight  path.  The  program  assumes  that  distances  are  expressed 
in  feet  and  noise  levels  in  EPNdB. 

Aircraft  classes,  based  on  physical  attributes  of  the  air¬ 
craft,  are  defined  in  terms  of  noise  and  performance  data. 
Classes  such  as  "four-engine  turbojet”  and  "two-engine  turbo¬ 
fan”  are  used.  Each  aircraft  class  is  subdivided  into  trip 
length  categories.  These  categories  allow  for  varying  per¬ 
formance  characteristics  within  a  single  aircraft  class.  The 
trip  length  category  for  each  aircraft  class  must  be  defined 
in  terms  of  an  altitude  profile,  Delta-EPNL  profile,  and  EPNL 
profile.  These  profiles  heed  not  be  unique  to  a  particular 
aircraft  trip  length.  Usually  a  modest  number  of  profiles  are 
sufficient  to  define  several  different  classes  of  aircraft. 
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The  level  of  flight  activity  is  the  volume  of  operations 
occurring  during  each  of  two  daily  time  periods.  These  time 
periods  include  the  daytime  hours  (0700-2200)  and  the  night¬ 
time  hours  (2200-0700).  Each  time  period  must  have  a  level 
of  flight  activity  specified. 

A  generalized  flow  diagram  of  the  input ,  output,  and  com¬ 
putation  stages  of  the  program  is  shown  in  Fig.  1.  The  data 
is  read  and  checks  are  made  to  determine  that  all  necessary 
data  ha«  been  entered.  If  no  input  errors  have  been  encountered 
the  program  will  proceed  to  compute  the  desired  contours. 
Processing  of  the  flight  paths  is  sequential  until  all  data 
has  been  processed. 

B .  Input  Data  Specifications  and  Formats 

Data  is  prepared  on  punched  cards  for  program  input.  A 
unique  format  exists  for  each  type  of  data.  These  formats 
are  detailed  in  the  paragraphs  below.  Examples  are  shown  in 
Fig.  A-l . 

1. )  AIRCRAFT  DESCRIPTOR 

An  aircraft  class  is  defined  in  terms  of  a  numeric  name, 
an  alphanumeric  descriptor  ana  the  numeric  names  of  the 
three  profiles  (altitude,  delta -EPNL,  and  EPNL)  to  be 
used  for  each  of  eight  trip  length  categories.  Category 
1  has  been  arbitrarily  assigned  to  landings  and  the 
remaining  categories  to  takeoffs.  Category  # 1  ordinarily 
uses  altitude  profile  #1  because  this  is  the  landing  pro¬ 
file  generated  by  the  program  as  a  function  of  the 
glide  slope  and  runway  length.  Table  I  shows  the  card 
format  for  entering  this  data. 

2. )  ALTITUDE  PROFILES 

An  Altitude  Profile  is  defined  in  terms  of  a  numeric  name, 
the  number  of  X,  Y  coordinate  values  describing  the  profile, 
and  the  coordinate  values.  The  X  coordinate  is  the 
distance  along  the  flight  track  and  the  Y  coordinate,  the 
altitude.  The  origin  of  the  coordinate  system  is  arbitrary. 
However,  it  is  typically  on  the  end  of  the  runway  at  which 
the  aircraft  commences  its  takeoff.  The  landing  profile 
is  generated  by  the  program,  making  this  assumption.  For 
flight  paths  which  involve  operations  on  more  than  one 
runway,  a  common  profile  origin  point  must  be  chosen. 

This  origin  is  usually  related  to  one  of  the  runways  and 
the  profiles  involving  operations  from  the  other  runways 
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FIGURE  1.  BLOCK  DIAGRAM  OF 
PROGRAM  FLOW 
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TABLE  I 

AIRCRAFT  DESCRIPTOR  CARD  FORMAT 


Columns 

Contents 

1-2 

Kumerlc  Name 

3-14 

Alphanumeric  Descriptor 

17-18 

19-20 

20-21 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

25-26 

27-28 

29-30 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

33-34 

35-36 

37-38 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

41-42 

45-46 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

49-50 

51-52 

53-54 

Altitude  Profile  Name 

Delta  EPNL  Profile  Name 

EPNL  Profile  Name 

57-58 

59-60 

61-62 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

65-66 

67-68 

69-70 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

73-74 

75-76 

77-78 

Altitude  Profile  Name 

Delta-EPNL  Profile  Name 

EPNL  Profile  Name 

Mode 

Range 

Integer 

1-99 

Alpha 

Trip  Category 

Integer 

1-99 

01  (Landings) 

ft 

1-99 

It 

1-99 

Trip  Category 

If 

1-99 

0 2  (Takeoffs 

ft 

1-99 

0-500  Ml) 

ft 

1-99 

Trip  Category 

If 

1-99 

03  (Takeoffs 

ft 

1-99 

500-1000  ml) 

tf 

1-99 

Trip  Category 

ff 

1-99 

#4  (Takeoffs 

ff 

1-99 

1000-1500  ml) 

ft 

1-99 

Trip  Category 

V 

1-99 

05  (Takeoffs 

ft 

1-99 

1500-2500  ml) 

1-99 

Trip  Category 

ff 

1-99 

06  (Takeoffs 

ff 

1-99 

2500-3500  ml) 

ff 

1-99 

Trip  Category 

ft 

1-99 

07  (Takeoffs 

ff 

1-99 

3500-4500  ml) 

ff 

1-99 

Trip  Category 

ff 

1-99 

08  (Takeoffs 

ft 

1-99 

4500  +  ml) 

ff 

1-99 
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TABLE  II 

i 

ALTITUDE  PROFILE  AND  DELTA- 

EPNL  PROFILE 

\  * 

CARD  FORMAT 

|  Column 

Contents 

Mode  Range 

if  I  1‘3 

Numeric  Name 

Integer  1-99 

1  4-6 

Number  of  points 

.Integer  1-10 

i  1 

to  describe  profile 

■  7-8 

Blank 

|  9-14 

xi 

Real 

15-20 

Yi 

II 

21-26 

X2 

Tl 

27-32 

Y2 

IT 

33-38 

X3 

II 

39-44 

Y3 

1! 

45-50 

x4 

II 

51-56 

y4 

ll 

57-62 

X5 

II 

63-68 

Y5 

II 

69-74 

VO 

X 

II 

75-80 

y6 

1 

| 

To  enter 

more  than  six  coordinates, 

use  an  additional  card 

leaving  columns  1-8  blank.  X7  must 

appear  in  columns  9**14,  | 

*  etc . 

i* 

1 

1 

•7- 

9 

i 

must  be  adjusted  for  the  difference  in  origins.  Although 
the  program  generates  the  coordinates  of  the  landing  pro¬ 
file,  the  user  must  define  a  dummy  profile  whose  numeric 
name  is  "l".  It  must  be  the  first  altitude  profile  de¬ 
fined  in  the  program.  The  restrictions  on  the  coordinate 
values  are  that  no  number  may  be  negative  and  successive 
X  values  must  always  be  increasing.  Twelve  coordinate 
values  are  allowed.  Table  II  shows  the  card  format  for 
entering  this  data. 

3. )  DELTA-EPNL  PROFILES 

A  Delta-EPNL  Profile  is  defined  in  terms  of  a  numeric  name, 
the  number  of  X,  Y  values  describing  the  profile,  and  the 
coordinate  values.  The  X  coordinate  is  the  position  on 
the  flight  track  and  the  Y  coordinate  is  the  EPNL  correction 
to  be  added  to  the  EPNL  computed  from  the  EPNL  vs  Slant 
Distance  function.  As  with  the  altitude  profile,  the 
origin  of  the  coordinate  system  is  arbitrary,  however,  the 
origin  for  all  altitude  and  Delta-EPNL  profiles  must  be 
the  same.  The  only  restriction  on  the  coordinate  values 
is  that  successive  X  values  must  always  be  Increasing. 

There  is  no  restriction  on  the  Y  value.  The  card  format 
for  this  profile  is  identical  to  the  altitude  profile 
format  and  is  shown  in  Table  TI . 

4. )  EPNL  PROFILE 

An  EPNL  vs  Slant  Distance  Profile  is  defined  in  terms  of 
a  numeric  name,  and  two  lists  of  35  EPNL  values  each. 

These  EPNL  values  correspond  to  assumed  values  of  Slant 
Distance  starting  at  12.6  feet  (1()1*1)  and  increasing 
by  a  factor  of  lO^-l  up  to  31622  feet  (10^*5).  One 
list  is  to  be  used  when  the  aircraft  is  on  or  near  the 
ground  and  the  sound  propagation  path  is  from  ground 
location  to  ground  location.  The  second  list  is  to  be 
used  when  the  aircraft  is  airborne  and  the  propagation 
path  is  from  air-to-ground.  An  Integer  code  is  used  to 
specify  the  propagation  path  associated  with  the  data 
(1  ®  ground-to-ground,  2  »  air-to-ground).  The  program 
expects  the  data  to  be  on  eight  successive  cards;  however, 
the  cards  may  be  in  any  order  since  the  rank  of  the  data 
is  also  coded  on  each  card.  Table  III  shows  the  card 
format  for  entering  the  EPNL  data. 


TABLE  III 


EPNL  PROFILE  CARD  FORMAT 


Column 

Contents 

Mode 

Range 

1-2 

Numeric  Name 

Integer 

1-99 

3-4 

Propagation  Path  Code 

i? 

1-2 

5-7 

Lowest  EPNL  rank  (LI) 

t» 

1,11,21,31 

8-10 

Highest  EPNL  rank 

M 

10,20,30,35 

11-17 

epnll1 

Real 

18-24 

EPNLL1  +  ! 

tl 

25-31 

EPKLL1  +  i 

!1 

32-33 

EPNLL1  +  3 

l! 

39-45 

EPNLL1  +  4 

It 

46-52 

EPNLL1  +  5 

It 

53-59 

EPNLL1  +  6 

!» 

60-66 

EPKLL1  +  7 

II 

67-73 

EPNLL1  *  8 

It 

74-80 

EpNLL1  +  9 

II 

TABLE  IV 

VOLUME  OP  OPERATIONS  CARD  FORMAT 


Columns 

Contents 

Mode 

Range 

1-3 

Aircraft  Class  Numeric  Name 

Alpha 

1-99 

4-6 

Time  Period  Code 

Integer 

1-2 

7-14 

Volume  For  Trip  Length 

Category  § 1  (landings) 

Real 

>0 

15-22 

#2  (Takeoffs) 

If 

>0 

23-30 

*3  (  "  ) 

ff 

>0 

31-38 

#4  (  ) 

If 

>0 

39-46 

#5  (  "  ) 

II 

>0 

47-54 

#6  (  "  ) 

If 

>0 

55-62 

#7  (  "  ) 

ff 

>0 

63-70 

#8  (  "  ) 

If 

>0 

TABLE  V 

FLIGHT  PATH  IDENTIFICATION  CARD 

FORMAT 

Columns 

Contents 

Mode 

Range 

1-51 

Alphanumeric  Descriptor 

Alpha 

>0 

52-58 

Path  Position  To  Commence 
Computation 

Real 

>0 

59-65 

Path  Position  To  Terminate 
Computation 

Real 

>0 

66-70 

Interval  On  Path  Between 
Computations 

Real 

*0 

71-76 

Runway  Length 

Real 

>0 

71-80 

Glide  Slope  (in  degrees) 

Real 

>0 
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5.)  VOLUME  OP  OPERATIONS 


For  each  flight  path  the  average  volume  of  operations  on 
the  path  must  be  specified  by  trip  length  category  for 
each  aircraft  class  to  be  considered.  The  numeric  name 
of  the  aircraft  class,  the  time  period  code  (1  *  daytime, 

2  *  nighttime)  and  the  volume  of  operations  for  each  trip 
length  category  must  be  specified.  The  specified  volume 
of  operations  may  never  be  negative.  Zero  values  are,  of 
course,  valid.  Two  cards  per  aircraft  class  are  expected 
by  the  program.  The  first  must  be  for  the  daytime  period 
and  the  second  for  the  nighttime  period.  Table  IV  shows 
the  card  format  for  this  data. 

6.)  FLIGHT  PATH  IDENTIFICATION 

In  order  to  compute  the  contours  for  a  given  flight  path, 
an  alphanumeric  descriptor  and  five  parameters  describing 
the  flight  path  are  required.*  The  parameters  specify 
(1)  the  point  on  the  track  (distance  from  start  of  take¬ 
off  roll)  to  commence  computing  the  contour,  (2)  the 
point  on  the  track  to  terminate  computing  the  contour, 

(3)  the  interval  along  the  track  at  which  to  compute  the 
contour,  (4)  the  length  of  the  runway  on  which  the  air¬ 
craft  arrived  or  departed  and  (5)  the  glide  slope  used 
by  the  aircraft  approaching  said  runway. 

The  computations  may  be  made  by  moving  in  either  direction 
along  the  flight  track  (l.e.  the  specified  interval  may 
be  either  positive  or  negative.  However,  care  should  be 
taken  to  insure  that  the  stopping  position  is  of  lesser 
value  than  the  starting  position  if  the  interval  is 
negative  and  conversely,  if  positive. 

The  runway  length  and  glide  slope  need  not  be  specified 
if  no  landing  operations  are  to  be  considered  on  the 
current  path  segment.  The  card  format  for  the  flight 
path  identification  data  is  shown  in  Table  V. 


*  For  convenience  and  ease  of  interpretation  all  input  and 
output  data  are  with  reference  to  the  flight  track  (i.e. 
the  vertical  projection  of  the  flight  path  on  the  ground). 
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C .  Organizing  the  Input  Data  for  Processing 

Data  to  be  read  Is  organized  Into  blocks  corresponding  to 
each  flight  path.  Each  block  Is  initiated  by  a  Flight  Path 
Identification  Card  and  terminated  by  a  Blank  Card.  The 
required  data  elements  (aircraft  descriptor,  noise,  and  per¬ 
formance  data  along  with  the  volume  of  operations  data)  are 
contained  between  these  two  cards.  The  only  limitation  on 
the  order  in  which  the  elements  appear  is  that  aircraft 
decriptor,  noise,  and  performance  data  must  precede  the 
volume  of  operations  data.  Only  those  aircraft  whose  volume 
of  operations  has  been  specified  in  this  block  will  be 
considered  in  computing  the  contours  for  this  path. 

The  program  assumes  all  data  cards  to  specify  volume  of 
operations.  If  other  information  is  to  be  read  in  a  control 
card  specifying  the  type  of  data  and/or  the  action  to  be 
taken  must  preceed  the  information.  The  control  card  format 
and  list  of  commands  is  given  in  Table  VI.  An  example  of 
the  card  deck  setup  to  compute  contours  for  a  flight  path 
is  shown  in  Figs.  A-2. 

In  most  cases  the  aircraft  descriptor,  noise  and  performance 
data  will  be  the  same  for  all  flight  paths.  Therefore,  this 
data  need  be  entered  only  once  (for  example,  with  the  data 
for  the  first  flight  path) . 

Since  data  storage  space  is  limited,  however,  the  capability 
to  delete  information  no  longer  needed  is  desirable.  To 
accomplish  this  certain  cards  (listed  in  Table  VI)  will  allow 
particular  areas  In  storage  to  be  erased  or  fixed.  A  large 
amount  of  data  will  be  used  over  and  over  in  the  computations 
for  a  multi-flight-path  Job.  These  data  can  then  be  followed 
by  an  appropriate  '’fix"  command.  When  at  a  later  time,  it  is 
desired  to  clear  part  of  storage  to  make  room  for  more 
special  profiles,  only  the  information  read  in  after  the  last 
fix  command  will  be  erased.  This  way  only  the  ’’transient" 
data  (i.e.  data  used  for  one  or  two  paths)  will  be  erased, 
without  the  need  of  re-entering  all  "standard"  profiles. 

In  this  manner  many  flight  paths  with  many  special  profiles 
can  be  handled  even  though  only  limited  storage  space  is 
available.  The  volume  of  operation  storage  area  is  automatically 
cleared  upon  entering  the  data  for  a  new  flight  path. 

If  an  error  is  detected  during  the  data  input  or  computation 
phase,  an  error  message  is  printed.  These  messages  are 
tabulated  in  Table  VII.  Unless  an  input  error  is  irrecover¬ 
able  data  input  will  continue;  however,  no  contour  will  be 
computed  for  the  flight  path. 
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TABLE  VI 

CONTROL  CARL  FORMAT  AND  CODES  FOR  ENTERING  AIRCRAFT 
DESCRIPTOR,  NOISE,  AND  PERFORMANCE  DATA 


Column  Contents  Range 

1-3  Operation  To  Be  Performed  See  Below 

4-6  Number  of  Times  Operation  Is  To 

Be  Repeated 


103 

104 


111 


112 

113 

114 

115 
121 
122 


123 

124 

125 
999 


Operation  To  Be  Performed 
Read  an  Aircraft  Descriptor 
Read  an  Altitude  Profile 
Reed  a  Delta-EPNL  Profile 
Read  an  EPNL  Profile 

Fix  currently  defined  aircraft  descriptors  in 
permanent  table. 

Fix  currently  defined  altitude  profiles  in 
permanent  table. 

Fix  currently  defined  Delta-EPNL  profiles  in 
permanent  table. 

Fix  currently  defined  EPNL  functions  in 
permanent  table. 

Fix  all  current  aircraft  descriptor,  noise 
and  performance  data. 

Delete  all  aircraft  descriptors  not  in 
permanent  table. 

Delete  all  altitude  profiles  not  in  permanent 
table . 

Delete  all  Delta-EPNL  profiles  not  in  per¬ 
manent  table. 

Delete  all  EPNL  functions  not  in  permanent 
table 

Delete  all  aircraft  descriptors,  noise,  and 
performance  data  not  in  permanent  tobies. 

Stop  all  processing. 
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TABLE  VII 
ERROR  MESSAGES 


Number  Error  Message 

1  ALTITUDE  PROFILE  TABLE  FULL 
FOLLOWED  BY  #4 

2  ALTITUDE  PROFILE  UNDEFINED 
FOLLOWED  BY  #4 


3  AIRCRAFT  DESCRIPTOR  TABLE 
FULL  FOLLOWED  BY  #4 

4  CANNOT  ENTER  NAME  - 


5  CARD  MISSING  OR  OUT  OF 
ORDER 


6  CONDENSED  A/C  TABLE  FULL 


7  DELTA-EPNL  PROFILE  UNDEFIN¬ 
ED  FOLLOWED  BY  #4 


3  DELTA-EPNL  PROFILE  TABLE 

FULL  FOLLOWED  BY  #4 

9  DUPLICATE  ALTITUDE  PROFILE 

NAME  FOLLOWED  BY  #4 


10  DUPLICATE  AIRCRAFT  DESCRIP¬ 
TOR  NAME  FOLLOWED  BY  #4 


11  DUPLICATE  DELTA-EPNL  PRO¬ 
FILE  NAME  FOLLOWED  BY  #4 


12  DUPLICATE  PNLEFF  PROFILE 
NAME  FOLLOWED  BY  #4 


13  EPNL  FUNCTION  UNDEFINED 
FOLLOWED  BY  #4 


14  FUNCTION  EPNL-INVALID 
INPUT  PARAMETER 


Cause 

The  Altitude  Profile  just 
entered  cannot  be  stored. 

An  Altitude  Profile  named  on 
the  Aircraft  Descriptor  Card 
has  not  been  entered. 

The  Aircraft  Descriptor  Just 
entered  cannot  be  stored. 

The  Numeric  Name  of  the  errone¬ 
ous  profile  is  printed. 

The  data  input  routine  has 
detected  an  error  in  the  card 
sequence. 

The  number  of  unique  aircraft 
is  greater  than  50. 

A  Delta-EPNL  Profile,  named 
on  the  Aircraft  Descriptor 
Card,  has  not  been  entered. 

The  Del"a-EPNL  Profile  just 
entered  cannot  be  stored. 

An  Altitude  Profile  of  the  same 
numeric  name  has  been  previous¬ 
ly  defined. 

An  Aircraft  class  of  the  same 
numeric  name  has  been  previous¬ 
ly  defined. 

A  Delta-EPNL  Profile  of  the 
same  numeric  name  has  been 
previously  defined. 

An  EPNL  Profile  of  the  same 
numeric  name  has  been  previous¬ 
ly  defined. 

The  EPNL  Profile  named  on  the 
Aircraft  Descriptor  card,  has 
not  been  entered. 

The  EPNL  calculation  sub¬ 
routing  has  seen  an  error  in 
the  calling  parameters. 
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TABLE  VII  (Con’t) 


Number  Error  Message 

15  GLIDE  SLOPE  CANNOT  BE 
NEGATIVE 


16  INTERVAL  IS  OR  WRONG  SIGN 


17  INVALID  NUMBER  OF  ALTITUDE 

POINTS  FOLLOWED  BY  #4 

18  INVALID  NUMBER  OF  DELTA- 

EPNL  POINTS  FOLLOWED  BY  #4 

19  INVALID  OPERATION  CODE 

FOLLOWED  BY  CODE 


20  INVALID  PNLEFF  PROPAGATION 
PATH  FOLLOWED  BY  #4 


21  INVALID  PNLEFF  RANK  LIMITS 

FOLLOWED  BY  #4 

22  PNLEFF  PROFILE  NAMES  DO  NOT 

MATCH  FOLLOWED  BY  #4 

23  PNLEFF  PROFILE  TABLE  FULL 

FOLLOWED  BY  #4 


24  PROGRAM  CANNOT  COMPUTE  THIS 
CONTOUR 


25  RUNWAY  LENGTH  MUST  BE 
GREATER  THAN  ZERO 


26  RUNWAY  TOO  SHORT  FOR  GLIDE 
SLOPE 

27  STARTING  POINT  CANNOT  BE 
NEGATIVE 


23  STOPPING  POINT  CANNOT  BE 
NEGATIVE 


29  UNDEFINED  FOLLOWED  BY  THE 
AIRCRAFT  CLASS  NUMERIC 
NAME  AND  THE  VOLUME  OF 
OPERATIONS  DATA 


Cause 

The  Glide  Slope,  on  the  Flight 
Pain  Identification  Card,  is 
in  error. 

The  start,  stop,  and  interval 
parameters,  on  the  Flight  Path 
Identification  Card,  are 
inconsistent . 

The  Altitude  Profile  has  too 
many  or  too  few  points. 

The  Delta-EPNL  Profile  has  too 
many  or  too  few  points. 

The  Operation  Code  on  the  Con¬ 
trol  Card  cannot  be  recognized 
by  the  program. 

The  propagation  path  code  on 
the  EPNL  Profile  Card  is  not 
1  or  2. 

The  rank  parameter  on  the  EPNL 
Profile  Card  is  invalid. 

The  EPNL  Profile  Cards  are 
inconsistenly  identified. 

The  EPNL  vs .  Slant  Distance 
function  just  entered  cannot 
be  stored. 

An  error  was  detected  during 
the  input  of  this  Flight  Path 
Block . 

The  runway  parameter,  on  the 
Flight  Path  Identification  Card, 
is  in  error. 

The  glide  slope  intersection 
is  off  the  runway. 

The  start  parameter,  on  the 
Flight  Fath  Identification 
Card,  is  invalid. 

The  stop  parameter  on  the  Flight 
Path  Identification  Card,  is 
invalid . 

The  aircraft  class,  on  the 
Volume  of  Operations  Card,  has 
not  been  defined. 
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III.  PROCEDURES  USED  FOR  COMPUTING  CONTOURS 


As  discussed  in  Section  I,  the  purpose  of  the  computer 
program  is  to  generate  the  distance  relationships  between 
the  flight  track  and  contours  of  equal  NEF  value.  This  sec¬ 
tion  of  the  report  discusses  the  techniques  used  to  determine 
these  distance  relationships.  Specifically,  the  program 
generates  the  perpendicular  (sideline)  distance  from  the 
flight  track  to  each  contour  at  regular  intervals  along  the 
track.  The  method  chosen  to  compute  these  distances  is  a 
heuristic  one  based  on  knowledge  of  the  NEF  and  distance  be¬ 
havioral  relationships. 

The  dashed  curve  in  Fig.  2  shows  a  typical  relationship 
between  the  NEF  value  and  sideline  distance  at  a  particular 
point  on  the  flight  track.  An  important  attribute  of  this 
curve  is  its  smooth  decrease  in  NEF  value  with  increasing 
distance.  The  iterative  technique  shown  graphically  in  Fig. 

2  is  used  to  determine  the  sideline  distance  to  the  desired 
NEF  contour.  Two  trial  distances,  Y^  and  Yt2>  are  chosen. 
The  corresponding  NEF  values  are  then  computed  for  these 
distances  and  a  straight  line  is  drawn  through  these  points. 
The  point  of  intersection  of  this  line  with  the  desired  NEF 
value  yields  a  new  trial  distance,  Yta.  The  NEF  value  at 
this  new  distance  is  then  computed.  The  new  (t3)  and 
previous  (t2)  trials  are  used,  as  above,  to  compute  another 
trial  distance.  The  iterative  procedure  continues  until  the 
difference  between  two  successive  trial  distances  becomes 
sufficiently  small.  The  perpendicular  distance  reported  is 
the  last  trial  distance  calculated.  The  computations  proceed 
from  point  to  point  along  the  flight  track.  The  change  in 
shape  of  the  NEF  vs  sideline  distance  curve  is  gradual  for 
successive  points  on  the  track.  Thus,  the  last  sideline 
distance  calculated  for  the  contour  is  used  to  determine  the 
first  two  trial  distances  for  the  next  point  along  the  flight 
track. 

The  NEF  contour  computation  is  terminated  at  the  end  of  the 
flight  track  or  when  the  contour  closes.  The  contour  is 
assumed  to  have  closed  if,  during  the  iterative  procedure,  a 
trial  value  is  less  than  12.6  feet.  The  point  of  closure  on 
the  flight  track  is  determined  by  calculating  the  NEF  at  two 
nearby  points  on  the  flight  track.  A  logarithmic  interpola¬ 
tion  is  performed  to  obtain  the  point  of  closure. 

In  the  above  procedure,  it  was  required  to  calculate  the 
NEF  value  for  a  given  sideline  distance.  This  NEF  value  is 
computed  as  the  sum,  on  an  energy  basis,  of  the  NEF 


contributions  from  the  operations  of  each  aircraft  trip  length. 
For  purposes  of  computation,  each  aircraft  trip  length  category 
is  treated  as  an  individual  aircraft  type  with  known  noise  and 
performance  characteristics.  Figure  3  shows  the  geometrical 
relationship  between  the  aircraft  and  ground  positions  for 
two  aircraft  types.  The  EPNL  at  the  ground  location  Y  is 
determined  from  the  slant  distance  to  the  aircraft.  The  NEF 
contribution  of  the  aircraft  is  then  evaluated  according  to 
Equation  1,  Appendix  B  from  the  EPNL  and  volume  of  operations 
for  this  aircraft.  The  NEF  contributions  for  all  aircraft 
are  then  summed  according  to  Equation  3>  Appendix  B.  The 
summed  NEF  is  the  actual  Noise  Exposure  Forecast  Value  at 
the  point. 

If  it  is  desired  to  change  the  numeric  value  of  the  NEF 
contour  to  be  computed  the  array  NEFCON  (Fig.  C-l^C)  can  be 
set  equal  to  the  desired  values.  This  cannot  be  done  without 
limit,  however,  The  EPNL  information  is  supplied  for  slant 
distances  up  to  31622  feet.  A  contour  which  has  a  larger 
distance  than  31622  feet  to  the  flight  path  will  therefore 
be  in  error.  (The  value  for  31622  feet  will  be  computed). 
Simarily  a  contour  with  a  value  such  that  the  slant  distance 
is  less  than  12.6  feet  will  also  be  in  error. 
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APPENDIX  A 


SAMPLE  NOISE  EXPOSURE  FORECAST  CONTOUR  CALCULATIONS 


APPENDIX  A 

SAMPLE  NOISE  EXPOSURE  FORECAST  CONTOUR  CALCULATIONS 


The  NEF  program  needs  for  Its  execution  several  types  of 
data.  All  data  Input  is  on  80  column  punched  cards.  The 
various  formats  for  the  different  types  of  data  are  shown  In 
Figs.  A-1A  through  A-1G. 

Some  examples  of  input  data  are  shown  in  Fig.  A-2.  The 
results  of  this  simple  computation  are  shown  in  the  Figs.  A~3. 

The  program  can  handle  considerably  more  complex  problems., 
such  as  the  example  given  in  Fig.  A-^.  The  complexity  of  the 
volume  of  operations,  the  total  number  of  operations  and  the 
length  of  flight  track  for  which  the  contours  are  computed 
affect  only  the  execution  time  of  the  program.  The  first 
example  in  Appendix  A  (Figs.  A-3)  has  7  entries  in  two  classes 
for  a  total  of  65  aircraft.  The  contours  are  computed  up  to 
30,000  feet.  Execution  time  for  this  example  is  65  seconds 
using  an  IBM  System/360  Model  30.  The  second  example  (Figs.  A-4 ) 
has  69  entries  in  eight  classes  for  a  total  of  3^2  aircraft. 

In  addition  the  contours  are  computed  to  100.000  feet.  The 
execution  time  in  this  case;  On  the  same  machine,  is  368 
seconds . 
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FIGURE  A  -  1G.  PERFORMANCE  DATA  CONTROL  CARD 
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FIGURE  A  -  28.  EXAMPLE  OF  INPUT  DATA 
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FIGURE  A  -  3  A  EXAMPLE  OF  COMPUTER  OUTPUT  -  FLIGHT  PATH 

VOLUME  OF  OPERATIONS  LISTING  (EXAMPLE  I) 
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FAST  MYTHICAL  MUNICIPAL  AIRPORT 


FORECAST  1970 


POINT  ON 


FLIGHT  TRACK 

25 

0 

5096 

2000 

4521 

4000 

4115 

60  00 

3607 

8000 

3362 

lnnoc 

3601 

1 2000 

4  72C 

1 4000 

5005 

16000 

5885 

18000 

6896 

20000 

685t> 

22000 

6800 

24C0C 

6751 

26000 

6687 

28000 

6615 

30000 

6534 

CONTOUR  CLOSING 
POINT  ON  FLIGHT 

TRACK 


DISTANCE  TO  CONTOUR  I 


30 

35 

40 

3336 

2  796 

1906 

3386 

2418 

1590 

2977 

2061 

1307 

2608 

1743 

1053 

2  394 

1  5  75 

994 

2386 

1968 

1192 

3146 

2260 

1466 

3841 

2487 

1394 

4077 

2422 

1286 

4042 

2333 

1134 

3921 

2220 

920 

3840 

2081 

599 

3744 

1  9l  1 

0 

3633 

1  7C3 

3505 

1445 

3360 

1115 

23463 


FEET 


45 

50 

55 

1173 

646 

321 

9  20 

497 

236 

732 

384 

178 

590 

302 

137 

689 

365 

164 

746 

463 

201 

809 

362 

n 

636 

472 

0 

0 

7559 

13807 

11634 

FIGURE  A  -3  B  . 


EXAMPLE  OF  COMPUTER  OUTPUT 
DISTANCE  LISTING  (EXAMPLE  1) 
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EAST  MVTHlCAt  PUNICIFAl  AIPPC8T 


rnPECAsr  tooo 


PH I  NT  ON 


>■  r  TP  AC  K 

23 
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20  CO 
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*  30? 

37,70 

526  7 
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1 

3447 

!  4070 

1137? 
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14233 

13700 

166’7 

20030 

1  °64  3 

2  ’700 
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?  40Ca 

l«647 

260CC 

I01f  3 
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13437 

?  o  7  ■:  o 
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??OCP 
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?  «00C 
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163cC 
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40700 
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5  311? 
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1  P*7C 

541CC 
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560CC 
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0  303” 

1  °1r  3 

600.1,7 
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1  7777 

647CC 
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1747  1 
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1724/ 

707.7'“ 

1  7Cf  7 

72  )?o 

16887 
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1/  644 
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16427 

73707 

1620? 

P  70  70 

1  3367 

82700 

13724 

P40C'’ 

13477 

86777 

1*21  1 

P37CC 

14C4  7 

C007O 

14664 

07  713 

14  77S 

r4-7?C 
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t’AC? 

1OO30O 

171  co 

r.riNToup  cnsiNr 
POINT  ON  *1  I  Gf-T 
THAT* 


riSTANCE  TC  CONTOUP  IN 


?C 

33 

40 

6165 

4  760 

3607 

5713 

4C43 

3112 

44  75 

7422 

2  567 

41*6 

?  204 

2362 

4*36 

3051 

3367 

6622 

5303 

16P8 

sprq 

6161 

3720 

1  of  67 

6  748 

1670 

1C  86  3 

6307 

1555 

1  1  137 

6 .708 

3301 

1  17C0 

6104 

3 1  °2 

11268 

6  C66 

2053 

1  1  733 

5010 

266  5 

11711 

5731 

2316 

11227 

556) 

1802 

11114 

*745 

1  371 

lOfCfc 

3  1  07 

7A4 

1C968 

4  830 

0 

1C773 

4323 

10*77 

4175 

10414 

1781 

1'12’0 

372T 

1  00  «  ? 

2707 

3f  *1 

2161 

3676 

1  766 

34  C6 

3  1  6  1 

8  P*  3 
P616 

fl7  j  5 

7S«.» 

7644 

7263 

6  f  62 
6413 
53?4 

1 

3  7  3  7 

4  3  ’  2 

41  n 
7271 
2173 

P17?6  40713  ’3173 


FFfT 


45 

50 

53 

2  300 

2021 

17  7-1 

2277 

1350 

051 

1787 

1144 
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1610 

106  5 

801 

2322 

1442 

p  en 

2121 

1361 
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FIGURE  A-48.  EXAMPLE  OF  COMPUTER  OUTPUT  -  CONTOUR 
DISTANCE  LISTING  (EXAMPLE  2) 
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APPENDIX  B 

SUMMARY  OP  BASIC  NOISE  EXPOSURE 
FORECAST  CONCEPTS  AND  EQUATIONS 


APPENDIX  B 


SUMMARY  OF  BASIC  NOISE  EXPOSURE 
FORECAST  CONCEPTS  AND  EQUATIONS 


The  NEF  procedures  have  evolved  from  predecessor  procedures 
described  in  references  1  and  2.  A  basic  assumption  in  these 
predecessor  procedures  which  also  underlies  several  of  the 
procedures  used  abroad  for  estimating  noise  exposures  is  that 
estimates  of  the  total  noise  environment  resulting  from  aircraft 
operations  must  be  based  upon  a  measure  of  noise  that  can  be 
closely  related  to  subjective  reactions  of  people  to  the  noise 
plus  corrections  to  take  into  account  the  number  of  occurrences 
of  noise  intrusions  per  daytime  and  per  nighttime  period. 

The  NEF  procedures  outlined  in  this  document  are  based  on  the 
parallel  studies  of  references  3  and  4.  The  current  procedure 
combines  elements  of  the  procedures  from  the  predecessor  studies 
hence  certain  calculation  steps  or  constants  given  in  this 
document  differ  from  these  given  in  either  reference  3  or  4. 

In  calculation  of  NEF,  aircraft  noise  levels  are  expressed 
in  terms  of  the  effective  perceived  noise  level  (EPNL)  as  defined 
in  reference  5.  The  EPNL  represents  an  integration  of  the  tone- 
corrected  perceived  noise  level  (determined  at  half  second 
intervals)  over  the  upper  10  PNdB  of  a  flyover  time  history. 

The  EPNL  value,  therefore  includes  adjustments  for  the  effective 
duration  of  the  flyover  signal  and  for  the  presence  of  discrete 
frequency  components. 

In  estimating  the  total  noise  exposure  near  an  airport  or 
flight  path  one  is  faced  with  predicting  or  calculating  the 
noise  levels  resulting  from  a  wide  range  of  a  variety  of  air¬ 
craft.  To  simplify  the  calculations  and  provide  a  systematic 
method  for  calculation,  it  is  convenient  to  group  the  aircraft 
in  classes  based  upon  consideration  of  the  aircraft  noise 
characteristics  and  takeoff  and  landing  performance.  Each  class 
is  assigned  a  description  of  the  noise  in  terms  of  a  set  of 
EPNL  vs  distance  curves  and  a  set  of  takeoff  and  landing 
profiles. 

In  estimating  the  effective  perceived  noise  levels,  it  is 
assumed  that  the  noise  level  at  any  point  under  or  to  one  side 
of  the  aircraft  flight  path  (or  to  either  side  of  the  path 
during  takeoff  roll)  depends  upon  the  following: 

a)  Type  of  aircraft  (particularly  the  type  of  engine), 
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b)  Type  of  operation  (takeoff,  landing,  or  takeoff  roll), 

c)  Engine  power  setting, 

d)  Distance  to  aircraft  (dependent  upon  the  aircraft  takeoff, 
and  landing  profiles  as  well  as  the  location  of  the  flight 
track  with  respect  to  the  ground  observer). 

For  a  given  class  of  aircraft  at  a  particular  power  setting 
(i.e.  takeoff  power)  it  is  assumed  that  the  aircraft  noise 
characteristics  may  be  described  by  a  single  EPNL  vs  distance 
curve . 

The  total  noise  exposure  produced  by  aircraft  operations  at 
a  given  point  is  viewed  as  being  composed  of  the  effective  per¬ 
ceived  noise  levels  produced  by  different  aircraft  classes  fly¬ 
ing  along  different  flight  paths.  For  aircraft  class  i  on 
flight  path  J,  the  NEF  (ij)  can  be  expressed  as 

NEF  (ij)  *  EPNL  (ij)  +  10  log  -C 

(Eq.  1) 

where 

NEF  (ij)  ■  Noise  Exposure  Forecast  value  produced  by  aircraft 
class  (i)  along  flight  path  segment  (j). 

EPNL  (ij)  *  Effective  perceived  noise  level  produced  at 
the  given  point  by  aircraft  class  (i)  flying 
along  flight  path  segment  (J) 

K  *  Constant  normalizing  the  adjustment  in  NEF  values  due  to 
volume  of  operations.  Different  values  of  K  are  used  for 
daytime  and  nighttime  movements. 

C  =  Arbitrary  normalization  constant. 

K  (day)  is  chosen  so  that  for  20  movements  of  a  given  aircraft 
class  per  daytime  period,  the  adjustment  for  number  of  operations 
is  zero.  Hence, 

pii 

io  log  r/ar;; t  -  o  k  (day)  -  20 
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K  (night)  is  chosen  such  that  for  the  same  average  number  of 
operations  per  hour  during  daytime  or  nighttime  periods  the 
NEF  value  for  nighttime  operations  would  be  10  units  higher 
than  for  daytime  operation.  Hence, 


10  «*  10  log 


K  (day) 

K  (night) 


9_ 

15 


where  9  and  15  are  the  number  of  hours  in  the  nighttime  and 
daytime  periods  respectively,  and 

K  (night)  =  1.2 

The  value  assigned  to  C  is  75.  Choice  of  this  value  is  based 
upon  two  considerations. 


First,  it  is  desirable  that  the  number  assigned  to  the  NEF 
values  be  distinctly  different  in  magnitude  from  the  effective 
perceived  noise  level  so  that  there  Is  little  likelihood  of 
confusing  effective  perceived  noise  levels  with  NEF  values.  A 
second  aspect  is  the  desirability  of  selecting  a  normalization 
factor  that  will  roughly  indicate  the  size  of  the  NEF  value 
above  some  threshold  value,  indicating  the  emergence  of  the  noise 
exposure  from  levels  which  would  have  little  or  no  influence 
on  most  types  of  land  usage. 


With  the  above  choices  for  values  of  K  and  C,  Eq.  (1)  becomes 
NEF  (ij)  «  EPNL  (ij) 

+  10  log  [N  (day)  (ij)  +  16.67  N  (night)  (ij)]-88 

(Eq.  2) 

The  total  NEF  at  the  given  ground  position  may  be  determined 
by  summation  on  an  energy  basis  of  all  the  individual  NEF  (ij) 
values : 


NEF  =  10  log  l  l  antilog  } 

i  J 


(Eq.  3) 


B-3 


I 

I 

I 

I 

J 

I 

I 

! 

APPENDIX  C 

FLOW  DIAGRAMS  AND  PROGRAM  LISTINGS 


Input  Necessary 
Flight  And  Performance 
Data  To  Compute 
A  Contour 


SET  PARAMETERS  FOR 

CONTOUR  COMPUTATION 

1.  Set  Position  On  Path  Equal 
To  Specified  Starting 
Position. 

2.  Set  Number  Of  Contours 
To  Compute. 

3.  Set  First  Trial  Contour 
Distances . 

4.  Clear  Contour  Closing 
Points. 

5.  Clear  Line  Count  For 
Output. 


FIGURE  C-l  A .  FLOW  DIAGRAM  -  NEF  CONTROL  ROUTINE 


C  -  1 


Save  New  Trial  Dtitance 


FIGURE  C-IC.  FLOW  DIAGRAM  -  NEF  CONTROL  ROUTINE 


INITIALIZATION 


OUSF  C  -  2  A  ,  SUBROUTINE  rjr 


lnltl#H»#  for 

Examination  Of 
Romolnlng  Alice 


HOUR'  C-20,  SUBROUTINE  CARtMN 


Ineramant 
Numb#i  Of 
Unique  Aircraft 


t- 


I>.»3  !>«<»•* 

•Hi  *!Mi 


'**  »1'  '  ’  •  c  CARDIN 


Read  Aircraft 
Descriptor  Card 

Expunge  *  And  p|oce  Oota  |n 


Expunge  Command  DeleteJ  Dafo  Table  Entriej  Nof  Fixed 
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FIGURE  C  -  3D  .  FLOW  D  !  A 


no 


Cloiing  Y£C 

Point  y — 


C  -  16 


LASTPT 


N  E  L  A  T  P  T 


Fraction  Graotar 


FIGURE  C  -  5  A  .  FLOW  DIAGRAM  -  SUBROUTINE  A  N  E  F 


5  a  g 

>  2  o 

U-  ^  X 

S°l 

III 

hu 


00- 


2.1 2  I 

it  J  •“  o 
0#£°‘ 
=  i  e  z 

o  o  e  a. 
«unS 


c  -  ■£ 
,?0“  c 

X  •»  U-  0 

U  4)  M  D  UJ  •  • 

8  I  1^2 j 

'SlfSli 

3i<5-a6  ° 


ig  t  z 

3  0  o 


FIGURE  r  -  5  B  FLOW  DIAGRAM  -  SUBROUTINE  ANEF 


Get  Entry  Number 
In  Delta  -  EPNL 
Profile  Table 


function  SUBPROGRAM  X  NEW 


P I  ~  X  Value  of  Desired  Pt. 
P2  -  L  ist  of  X,  Y  points  on 
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NEF  CONTOUR  CONTROL  RGUTINF 
PURPOSE 

TO  COMPUTE  THE  COORDINATES  OF  NEF  CONTOURS  AS  A  FUNCTION 
OF  AIRCRAFT  FLIGHT  PATH  PARAMETERS  AND  AIRCRAFT  PERFORMANCE 

DESCRIPTION  OF  PARAMETERS 

X  STAR  T  -  LOCATION  ON  FLIGHT  PATH  TO  START 
CONTOUR  COMPUTATIONS  I  REAL) 

XEND  -  LOCATION  ON  FLIGHT  PATH  TO  END 
CONTOUR  COMPUTATIONS  (REAL) 

DX  -  INTERVAL  ALONG  FLIGHT  PATH  BETWEEN 
CONTOUR  COMPUTATIONS  (REAL) 

NAC  -  NUMBER  OF  UNIQUE  AIRCRAFT  TYPES  (INTEGER) 

CONTUP  -  VALUE  OF  NEF  CONTOUR  CURRENTLY 
BEING  COMPUTEO  (REAL! 

NTET(I)  -  TOTAL  NUMBER  OF  CURRENT  ENTRIES 
IN  DATA  TABLES  (INTEGER) 

I  =  1  -  AIRCRAFT  DESCRIPTOR  TABLE 

2  -  ALTITUDE  PROFILE  TABLE 

3  -  DEL  TA-E  PNL  PROFILE  TABLE 

4  -  E PNL  PROFILE  TABLES 
NFET(I)  -  NUMBER  OF  ENTRIES  IN  DATA 

TABLES  WHICH  ARE  FIXED  AND  MAY  NOT 
BE  DELETED  (INTEGER) 

I  =  I  -  AIRCRAFT  DESCRIPTOR  TABLE 

2  -  ALTITUDE  PROFILE  TABLE 

3  -  DELTA-E  PNL  PROFILE  TABLE 

4  -  E  PNL  PROFILE  TABLES 

NMAX(I)  -  MAXIMUM  ALLOWABLE  NUMBER  OF  ENTRIES 

IN  THE  VARIOUS  DATA  TA8LES  (INTEGER) 

I  =  L  -  AIRCRAFT  DESCRIPTOR  TABLE 

2  -  ALTITUDE  PROFILE  TABLE 

3  -  DELTA-EPNt  PROFILE  TABLE 

4  -  E PNL  PROFILE  TABLES 

DSCNAMU)  -  NUMBER  IC  NAME  OF  I-TH  ENTRY  IN 

THE  AIRCRAFT  DESCRIPTOR  TABLE  (INTEGER) 
DSCRPT(J,I)  -  12  CHARACTER  * LPHA NUMH R I C 

AIRCRAFT  DESCRIPTOR  FOR  THE  I-TH 
ENTRY  IN  THE  AIRCRAFT  CLASS 

DESCRIPTOR  TABLE.  J  IS  DIMENSIONED  FOR  3  TO 
ALLOW  FOR  THE  12  CHARACTERS.  (INTEGER) 
ALTUIk(J,I)  -  DIRECTORY  WHICH  CONTAINS  THE  NUMERIC  NAMES  OF 

ALTITUDE  PROFILES  USED  BY  THE  I-TH  ENTRY  IN 
THE  AIRCRAFT  CLASS  DESCRIPTOR  TABLE  FOR  THE 
J-TH  TRIP  LENGTH  CATEGORY.  (INTEGER) 
PGWDIR(Jfl)  -  DIRECTORY  WHICH  CONTAINS  THE  NUMERIC  NAMES  OF 

DEL TA-E  PNL  PROFILES  USED  BY  THE  I-TH  ENTRY  IN 
THE  AIRCRAFT  CLASS  DESCRIPTOR  TABLE  FOR  THE 
J-TH  TRIP  LENGTH  CATEGORY.  (INTEGER) 

FPNDIE  ( J,  l  )  -  DIRECTORY  WHICH  CONTAINS  THE  NUMERIC  NAMES  OF 


FIGURE  C-M A.  PROGRAM  LISTING  -  CONTROL  ROUTINE 
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EPNL  PROFILES  USED  3Y  THE  i-TH  tNIRV  IN  THE 
AIRCRAFT  CLASS  DESCRIPTOR  TABLE  FOP  THE  J-TH 
TP  IP  LENGTH  CATFGCRY.  (INTEGER) 

ALTNAM(-)  -  Nu.  EPIC  NAME  OF  THF  I-TH  ENTRY  IN  THE  ALTITUDE 
PROFILE  TABLE.  (INTEGER  I 

NALTI?)  -  NUMBER  OF  COOPOI  NATES  DEFINING  THE  I-TH  ENT  &  Y  IN 
THE  ALTITUDE  TABLE  (  INTEGER) 

ALTtKfJ.I)  -  THE  ALTITUDE  PROFILE  VLAUt  F OF  THE  I-TH  ENTRY 
IN  THF  TABLE  FOP  THE  J-TH  COOROI NAT E( 1=X , 2=4 ) 
rt NO  THE  K-TH  POINT  IN  THt  PROFILE  (PEAL) 

2(1)  -  THE  ALTITUjF  AT  THE  CURRENT  POINT  UN  THE  FLIGHT  PATH, 

EVALUATtJ  FPfM  THE  I-TH  ENTRY  IN  THE  ALTITUDE 
PROFILE  TABLE .  (REAL) 

POKNAHt!)  -  NUMERIC  AN 

POWNAMIi)  -  NUMERIC  NAME  OF  THE  I-TH  ENTRY  IN  THE  DECT  A-EPNL 
PROFILE  TABLE.  (INTEGER) 

NPLh(I)  -  NUMBER  Or  COORDINATES  DEFINING  THE  I-TH  ENTRY  IN 
THF  OELTA--FPNI  TABLE  (INTEGER) 

POWER  (  K»  J,  I  )  -  THE  OELTA--PNL  PROFILE  VALUE  FOR  THE  I-TH 

ENTRY  IN  THE  TABLE  FOR  THE  J-TH  COORDINATE 
{ 1= X,2=Y) AND  THE  K-TH  POINT  IN  THE  PROFILE 
(REAL) 

POWSPL(I)  -  THt  JE  L  TA-E  PNL  VALUE  AT  THE  CURRENT  POINT  ON  THF 
FL 1 3-  PATH,  EVALUATED  FROM  THF  I-TH  ENTRY  IN 
THE  DtcTA-EPNL  PROFILE  TABLE  (PEAL) 

PNLNAM(I)  -  NUMERIC  NAME  OF  THE  I-TH  ENTRY  IN  THE  EPNL 

profile  table,  (integer) 

PNLEFF ( K» J , I  3  -  THE  EPNL  PROFILE  VALUE  FOR  THE  I-TH  ENTRY  IN 

THE  TABLE  FOR  THE  J-TH  PROPAGATION  PATH 
(1--  GROUND-TO-GROUND,  2 =A [ R-T C-GR OU ND)  AND 
THt  K-TH  VALUE  IN  THE  PROFILE 

ACALC(I)  -  LIST  OF  ALTITUDE  PROFILES  REQUIRED  FOR  COMPUTING 
CONTOURS  ON  CURRENT  FLIGHT  PATH  (INTEGER) 

PCALC(I)  -  LIST  OF  DELTA-EPNL  PROFILES 

REQUIRED  FOR  COMPUTING  CONTOURS  ON  CJRRENT  FLIGHT 
PATH  (INTEGER) 

OPER(I)  -  NEF  ADJUSTMENT  FOR  VOLUME  OF-  OPERATIONS  ASSOCIATED 
WITH  THE  I-TH  UNIQUE  AIRCRAFT  TYPE  (REAL) 

APTP(I)  -  POINTER  TO  THE  ENTRY  IN  THE  ALTITUDE  PROFILE  TABLE 
USED  FOR  THE  I-TH  UNIQUE  AIRCRAFT  TYPE.  (INTEGER) 
PPTR(I)  -  POINTER  TO  THF.  ENTRY  IN  THE  DELTA-EPNL  PROFILF 
TAttLE  USED  FOR  THt  I-TH  UNIQUE  AIRCRAFT  TYPE 
(  INTEGER  I 

EPTR(I)  -  POINTER  TO  THE  FNTRY  IN  THE  EPNL  PROFILE  TABLE  USE') 
FOR  THE  I-TH  UNIQUE  AIRCRAFT  TYPE.  (INTEGER) 

I  DENT!  I )  -  ALPHANUMERIC  DESCRIPTOR  STRING  FOR  FLIGHT 

PATH  IDENTIFICATION.  (ALPHA) 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  CALLED 
CARDIN,  APO,  ANFF  ,  XNF  W ,  LASTPT 


(  ! J  M  f '  N  >  STAR  T,  XI  OX,  NAC  *  LONTU'R  ,  MAX  (4)  COM  l 

(  N  N  T  t.  T  (  4  )  ,  M  E  T  (  v  )  ,  ,AI.  T(12),  \POW<  12  )  ,  ACALCU2),  PCALCC12)  COM  2 

COM  <!»N  USCNA-K  A. )  ,  AlfNAMdP),  pi  !„i,AM(  1  2  )  .  PNL  NAM  (25  )  *  APT  K  (*><">  )  COM  3 

FIGURE  C-113.  PROGRAM  LISTING  -  CONTROL  ROUTINE 
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FOKTMAIN 


COM  4 
COM  5 
COM  6 


N  EF  CPN  (  1 )  s  2  5 
to EF  CON ( 2  )  »  30 
NEFCl>N(  3  )  s  35 
toEFCON( 4 )  *  4.0 
toEFCOM  5)  =  45 
NtFC0N<6>  =  50 
to  EFf.ONt  7  )  *  55 

to  M  A  X  (  1 )  =  20 
NMAX(2)  *  12 
NM\X( 3 )  «  12 
toMA  X  ( »t )  a  25 
00  l  1-1,4 
NTET( I )  =  0 

1  to  F  F  T  (  I  )  =  0 

2  CALL  CAPO  IN 

X  «  XSTAP  T 
NNFXT  *  7 
L  IN  CM  T  *  0 
00  3  N*  l ,  7 
XCL  OS  L  ( to  )  a  r. 

’  I  Y( N )  a  4000 

I  1  CALI  APPI  X) 
to  CON  a  totoFXT 

(•>l  40  N=l,NCON 
I  ON  TUP  *  toEK.ON(N) 

Y?  »  IY(N  ) 

Y 1  a  Y 2  *  0.707 
CALL  ANFF  (Y2,  NEF2) 

CALL  AN  CP  (Yl,  NtFH 

21  Y*  =  XNf*(Vt,  to E F 1  ,  Y2,  NFF2,  CCNTUR) 

IF  ( Y 0  -  3O00O.J  2j,  23,  22 

22  Y*1  a  3000  0, 

21  Cn.NT  ItoUt 

IF  I Yf  -  12.6)  24,  25,  25 
24  (  IN  1  (tour 

CALL  LASTPrix,  XCLlISFINI) 

V'  =  0. 

A 'IF  XT  a  NNFXT  -  1 

GO  T(  2:r 

2r>  H  (AAS(Y"i  -  Yl)  -  2.)  2 H ,  20,  ?6 

?F  IF  (AIiSIIYT  -  Yl)  /  YO)  -  0.001)  28,  20,  27 
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COMMON  PPTkf  50 ) ,  EPTPI50),  0PEP(50).  2(12),  POtoSPL  (12).  IfJFNTlm 

rnMM^  T  nsCRPT( 5,20) ,  ALTOIR( 0,20) ,  POWOI R{ 0 , 20 ) ,  EPNOIB(0,2O) 
COMMON  PNLCFF(  35,  2, 25) ,  AL T( 10,2 , 1 2 ) ,  PO WE R ( 10 ,2 , 1 2 ) 

UIMCNS  ION  tot'  FC ON <  7  )  ,  IY(7),  XCL0SE(7) 

FEAL  N  F  T  1 ,  NfF  2 
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FORT MAIN 


r 

c 

c 


c 

c 

r. 

c 


r. 


< 


< 


27  Y2  =  Y 1 
NEF2  =  NCFi 
VI  =  VO 

CALL  ANEF  (VI,  NFF  1 ) 
CO  TO  21 

28  l  Y(  N )  =  YO  ♦  0,5 
4*  CONTINUF 


L  INCN  T  =  L INC NT  *  1 
IF  (LINCNT)  A3,  42,  4  2 

42  L  1NCNT  *  -5C 

Y.P  I  T  E  (3,  3004)  IDENT 
WRITE  (3,  30C l )  NEFCON 

43  IX  «  X  +  0,5 

WRITE  (3,  3005)  IX,  (IY(N),  N»l  ,NCON) 
IF  (NNF.XT)  55,  55,  53 
53  X  =  X  +  OX 

IF  (X  -  XFND)  11,  11,  55 


55 


60 

61 

6? 

<!? 

64 


66 

67 

64 

7- 

71 


I  FLO,  =  ?. 

WPITF  (3,  3006) 
m  70  Ns  1,7 
IX  «  XCLUSE(N  )  ♦  C.5 
IF  (IX)  70,  70,  60 
GO  TO  (61,  62,  63, 
(3,  3031)  IX 


64,  65,  66,  67),  N 


WRITE 

go  to 

WRITE 
GO  TC 
WRITE 
GO  T( 
WPITF 
GO  Tf 

V  R  I  TF 
C.(j  TC' 

V  a  II  C 
GO  TO 

w«  i  tt 
I  EL  G 


6*7 
(  3. 
6« 

(  3, 
69 
(  3, 
69 

(  3, 

65 
(  3, 
6° 

(  3, 
■  1 


30  32 ) 
JC  3  3 ) 


IX 


I  X 


3C34)  IX 


3036 ) 
3'?6) 
3037  ) 


(UN  T  It.UF 
U1  TC  (71, 
VF  |1(  IS 


72  ) , 
3*  30  ) 


IX 
I  X 
I  X 

I  F  LG 


'2  G'l  TC  ? 

1  (CROAT  (•'.',  >?X,  'UlbTANCF  TO  CONTOUR  IN  FELT 

1  ox,  •ELIGMI  TRACK'  ,  71  !J  /  5X,  - - 

. . —  -  -  - ') 

4  EUR  WAT  (MS  15X,  13A4  /  9X,  ' - — - — 


/  11X,  '.MINT  ON'/ 
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FOUTrtAIN 


/  *■*  l. 

1 - 

h  0  f-’  *'AT 

(  10X, 

18, 

4X,  718) 

< 

<  « 

FORMAT 

(fwMAT 

(  ) 

7X, 

•CUNTOUP  CLOSING'  /  8X 

V  M 

l  /  fix, 

(  Hu  mat 

•  TRACK'  ) 
('+',  2 1  X, 

18) 

v  }?. 

FOP  VAT 

<  '  ♦ 

29X, 

1  3) 

v  <•> 

tor  vat 

<*♦'» 

37X, 

18) 

•i  i  it. 

rrm  oat 

(  '  ♦  •  . 

48X, 

18) 

*1  '•  ■>  'j 

hit  vat 

83X, 

18) 

\  ^  t. 

HIP  MAT 

(  » 

6 1  X, 

18) 

v  37 

FOR. VAT 

69  X, 

18) 

END 
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SUBROUTINE  CARDIN 
PURPOSE 

TO  READ  DATA  NECESSARY  TO  COMPUTE 
NEF  CONTOURS  FOR  A  FLIGHT  PATH 

USAGE 

CALL  CARDIN 
REMARKS 

THIS  ROUTINE  CONTROLS  THE  INPUT  OF  ALL  DATA  TO  THE  PROGRAM. 
THE  ROUTINE  READS  THE  BLOCK  OF  DATA  FOR  A  SINGLE  FLIGHT 
PATH  AND  PERFORMS  VARIOUS  CHECKS  ON  THE  DATA  TO  TRAP  LOGICAL 
ERRORS  WHICH  WOULD  RESULT  IN  ERRONEOUS  COMPUTATIONS. 

IF  ERRORS  ARE  FOUND  IN  THE  DATA,  THE  ROUTINE  AUTOMATICALLY 
READS  THE  BLOCK  OF  DATA  FOR  THE  NEXT  FLIGHT  PATH.  THE 
ROUTINE  DOES  NOT  RETURN  UNTIL  IT  HAS  SUCCESSFULLY  RE  AO  A 
COMPLETE  BLOCK  OF  DATA. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
SIN,  COS,  PDATA,  OPCOR 


SUBROUTINE  CARDIN 


COMMON  X  START,  XEND,  DX,  NAC  ,  CONTUR,  NMAXI4)  COM  l 

COMMON  NTETI  4 1  ,  NFETI4),  NALTI12),  NPOW112),  AC  ALC  (12),  PCALCU2)  COM  2 

COMMON  OSCNAMI 20) ,  ALTNAMI 12 ) »  POWNAMI 12),  PNLNAMI25),  APTRI50I  COM  3 

COMMON  PPTR ( 50 ) ,  EPTK(50),  0PERI50),  ZI12J,  PGWSPL(12),  I  0  6  NT ( 131  COM  4 

COMMON  DSCRPTI 3,20) ,  ALTDIRI  8,20) ,  POWDI R( 8 ,20 ) ,  EPNDIR(8,20)  COM  5 

COMMON  PNLEFFI  35,2,25)  ,  ALT( 10 ,2 ,12 ) ,  POWE R ( 10 , 2 , 12 )  COM  6 


INTEGER  ERRFLG,  ACALC,  PCALC,  EXPCRD,  TRAP,  ALTDIR,  POWDIR 
INTEGER  EPNDIR,  ALTNAM,  POWNAM,  PNLNAM,  CARu,  DSCNAM,  DSCRPT 
INTEGER  APTR,  PPTR,  EPTR 
REAL  NFLTS,  NITFLT 

DIMENSION  BUFF'  {  8 )  ,  TEMP(2,8),  NF  LT  S  (  2 , 8 , 20  ) 

NSEG  =  8 

-  INITIALIZATION 

1.  CLEAR  ERROR  FLAG 

2.  SET  LINE  COUNT  TO  -5 

3.  CLEAR  NUMBER  OF  CONDENSED  AIRCRAFT 

4.  CLEAR  NUMBER  OF  FLIGHTS 

5.  CLEAR  LIST  UP  REQUIRED  ALTITUDE  PROFILES 

6.  CLEAR  LIST  OF  REQUIRED  OELTA-EPNL  PROFILES 
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CARDIN 


2  NAC  =  0 
ERR  FL  G  =  I 
LINCNT  =  -5 
CO  6  1=1,2 

DO  6  J=  1 , 8 
DO  6  K= 1  *  20 
6  NFLTSI  I,J,K)  *  0. 

DO  10  1-1,12 
ACALCI I)  «  C 
10  PCALCI I )  =  0 

-  READ  CARD  WITH  FLIGHT  PATH  IDENTIFICATION  AND  PARAMETERS 

READ  (1*  10011  IDENT,  X  STAR  T,  XEND,  DX,  RWL,  GS 

-  WRITE  FLIGHT  PATH  IDENTIFICATION 

WRITE  (  3,  3001)  IDENT 
WRITE  ( 3,  3000) 

- ARE  ALL  PARAMETERS  ZERO 

IF  (XSTART  ♦  XENO  +  DX  +  RWL  +  GS)  15,  15,  20 

-  SET  ERROR  FLAG 

15  ERRFLG  *  0 
GO  TP  50 

-  WRITE  FLIGHT  PATH  PARAMETERS 

20  WRITE  (3,  3002)  XSTART,  RWL,  XEND,  GS,  DX 
WRITE  (  3,  3000) 

WRITE  (  3,  3003) 

-  STARTING  POSITION  OF  PATH  NEGATIVE 

TF  (XSTART)  22,  23,  23 

- WRITE  ERROR  AND  SET  FLAG 

22  WRITE  (3,  MCI) 

ERR  FIG  =  0 


-  ENDING  POSITION  ON  PATH  NEGATIVE 

23  IF  (XEND)  24,  25,  25 
- WRITE  ERR  OR  AND  StT  FLAG 

24  WRITE  (  "»♦  3<T?| 
fRRFLG  *  C 

C 

C  - ENDING  POSITION  GREATER  THAN  STARTING  POSITION 

FIGURE  C-128.  PROGRAM  LISTING  -  SUBROUTINE  CARDIN 
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25  IF  'XEND  -  XSTART)  26,  30,  20 
-  IS  INCREMENT  ALONG  PATH  NEGATIVE 

26  IF  (0X1  30,  27,  27 

- WRITE  ERROR  AND  SET  ERROR  FLAG 

27  WRITE  13,  39C3I 
FRPFLG  =  0 

GO  TO  3C 

-  IS  INCREMENT  ALONG  PATH  POSITIVE 

28  IF  (DXI  29,  29,  30 

- WRITE  ERROR  AND  SET  ERROR  FLAG 

29  WRITE  <  3,  3903) 

ERPFLG  *  0 

- CREATE  ALTITUDE  PROFILE  FOR  LANDINGS.  SET  COORDINATES  OF 

FIRST  POINT.  SET  NUMBER  OF  PROFILE  POINTS  TO  I 

30  NALT<  l)  =  l 

ALT  (  I,  I,  l  )  =  0. 

ALT  (  1,2,  1)  =  0. 

- ARE  PUNWAY  LENGTH  AND  GLIDE  SLOPE  BOTH  ZERO 

IF  CRWL  ♦  GS)  32,  50,  3? 

-  TEST  VALUE  OF  GLIDE  SLOPE 

3?  IF  (GS)  33,  50,  35 

- WRITE  FR.PPK  MESSAGE  AND  SET  ERROR  FLAG 

33  WRITE  (3,  39' 7) 

ERR  FI  G  =  P 
co  Tr  5r 

- COMPUTE  SECOND  A 'ill  THIRD  COORDINATE  PCI  NTS  LN  LANDING  PROFILE 

F  F  f)M  GLIDE  SLOP!  AND  MINE  AY  LENGTH 

36  GS  =  1.745m  -  >  #  .S 

TANGS  *  SIMTS)  /  Cl  SINS) 

ALT  12,1,1)  =  FWL  -  5',/IANf.S 
ALT (?, 2, 1  )  = 

ALT  (3,1,1)  =  R  wL  ♦  1  •  r  6 

ALT  (1,2,1  )  =  1  .E6  *  TANGS  ♦  \ 

-  IS  X  VALUL  OF  SlCCf.  i  IMINT  POSITIVE 
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CARDIN 


IF  ( AL T(  2 ,  1 ,  1  )  )  37,  37,  39 

- WRITE  ERROR  ANO  SET  ERROR  FLAG 

37  WRITE  (3,  3905) 

ERR  FIG  =  0 
GO  TC  50 

- SET  NUMBER  OF  PRUFILE  POINTS  TO  3 


39 

NALTi  1  )  = 

3 

50 

READ  (1, 

1002) 

IOPI, 

l 0P2 ,  BUFF 

501 

IF  I  IOPI) 

100, 

100, 

51 

51 

IF  (IOPI 

-  99) 

20C, 

2CG,  52 

52 

FXPCPD  = 

l 

IF  I  IOPI 

-  199) 

300, 

300,  53 

53 

IF  I  IOPI 

-  99  9) 

54, 

55,  55 

54 

WRITE  I  3, 

GO  TC  50 

3908) 

IOPI 

55 

STOP 

¥ 

i 

I 

l 


i 

I 


30c  iop  =  iopi  -  ion 

N  =  IOP 2 

IF  (N)  305,  305,  306 

305  N  =  1 
C 

306  DO  3 1C  1=1, N 
CALL  POATAI  IOP  ) 

IF  COP)  307,  337,  310 

307  EPRFLG  =  0 
310  CONTINUE 

LINCNT  =  0 
GO  TO  50 
C 
C 

pro  CARD  =  IOP 2 

IF  (CARO  -  1)  2n2t  206,  202 
20?  WRITE  (3,  3914)  IOPI,  I GP2 
ERR  FIG  =  0 
CO  TO  5f 
C 

■*0 6  NAME  =  IOPI 

CO  2^4  I  =  1 ,N SF G 
204  TEMPI  1,1)  =  BUFFI  I  ) 

C 

FEAT  (1,  IOC? )  IUP1,  1 0P2 ,  BUFF 

L 

IF  INAMf  -  IOPI)  2C»B ,  2C5,  2f8 
;r P  v  A  I T t  I  3 ,  35  14)  IOPI,  ICP2 

f  rp  h  r,  =  o 
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GO  TO  5C  i 
C 

705  IF  (  IOP  2  -  2)  200*  209,  200 
20c  DO  2IC  1=  l,NSEG 
210  TEMPI  2,1)  -  BUFFI  I  ) 

C 

L  INCNT  =  LINCNT  +  1 
C 

IF  (LINCNT)  2 1  B»  216,  216 
216  WRITE  13.  3001)  IDfcNT 
WRITE  13,  300  C ) 

WRITE  13,  3003) 

LINCNT  -  -5 
C 

218  N  =  NTETI l ) 

DO  220  I  ENT  =  1,N 

IF  (NAME-  DSCNAMI IfcNT) )  220,  224,  220 

220  CONTINUE 
C 

WRITE  13,  3909 )  NAME,  (J,  ITEMPU*J)»  1*1, 2),  J=1»NSEG) 
ERR  FL  G  *  0 
GO  TO  50 
C 

224  WRITE  (3,  3011)  NA  ME ,  (OSCRPTII,  IENT),  1=1,3), 

1  <J,  I  TEMPI  I ♦ J ) ,  1  =  1,2),  J=1 »NSEG) 

C 

TRAP  =  1 

DO  238  1=  1, N SEG 

NAME  =  ALTO  IP  I  I,  IENT) 

N  =  NTFTI2) 

DO  228  J*1,N 

IF  I  NAME  -  ALTNAMIJ))  228,  23C,  228 
228  CONTINUE 

WRITE  (3,  3910)  NAME 
TRAP  =  0 

230  NAME  =  POWOIRII.  IENT) 

N  =  N TETI 3) 

DO  232  J=l,N 

IF  (NAME  -  POWNAMl J )  )  232,  234,  232 
2  32  CONTINUE 

WRITE  I  3,  3Q  i  1 )  NAME 
TRAP  =  0 

2  34  NAME  =  FPNOIPI  I,  IENT) 

N  =  NTETI4) 

DO  236  J*1,N 

IF  (NAME  -  PNLNAMIJ))  2  3fc,  23  8,  23t 
236  CONTINUE 

WRITF  (3,  3<H2)  NAMt 
TRAP  =  0 
2  38  CONTINUE 
C 

IF  (TRAP)  2?9,  2  39,  244 
2  39  F&P  R  G  =  r 

r,o  to  s  r 

c 
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244  00  246  J= 1»  2 

DO  246  1=  If N SEG 

246  NFITSIJfl,  IENT>  =  NFLTS  (  J,  T  f  I  ENT)  ♦  TEHPIJ,!) 

GO  TO  50 
C 
C 

500  STOP 
C 

100  CONTINUE 

102  IF  (ERRFLG)  103,  103,  104 

103  WRITE  ( 3,  3913) 

GO  TO  2 

104  NN  «  NTETI l) 

DO  150  K*i,NN 

DO  150  J* IfNSEG  I 

IF  (NFLTSl  1,J,K)  ♦  NFLTSl  2  ,J , K)  )  150,  150,  105 

105  1A2  =  ALTDIRI  J,K )  . 

I  A3  *  POWDIR I  J ,K ) 

1A4  =  EPNDIRI J,K) 

c 

TOTALD  *  0.  | 

TOTALN  *  0.  ' 

DO  112  K  1*K»  NN 

DO  112  J 1=  l,NSEG  { 

DAYFLT  »  NFLTSl  lfJlfKl)  J 

NITFLT  *  NFLTSl  2,J1,K1) 

IF  ( DAYFLT  ♦  NITFLT)  110,  110,  106 
1C6  IF  (  IA2  -  ALTDIRI  Jl.KU  )  112,  107,  112 

107  IF  ( I  A3  -  POWDIRl  Jl.KU  )  112,  108,  112 

108  IF  (IA4  -  EPNDIRI  Jl.KU)  112,  110,  112 

C 

110  TOTALO  *  TOTALO  ♦  OAYFLT 

TOTALN  *  TOTALN  ♦  NUFLT 
NaTSI  l, Jl.KU  *  0. 

NFLTSI2.J  1,K1  I  *  0. 

112  CONTINUE 

IF  INAC  -  SOI  118,  118,  115 
US  NAC  «  1 

ERR  FIG  »  0 
WRITE  <3,  3915) 

118  CONTINUE 

NAC  *  NAC  ♦  ) 

C 

OPERINAC)  *  OPCOR I TOTALO,  TOTALN) 

C 

N  *  NTETI  2) 

CO  121  NENT*  1  »N 

IF  I  IA2  -  ALTNAMINENT) )  121,  126,  121 

121  CONTINUE 
STOP  l 

126  DO  125  L«l,N 

IF  (ACALCIL)I  123,  123,  122 

122  IF  (ACALCIL)  -  NENT)  123,  124,  125 
125  CONTINUE 

123  ACalCIL 1  *  NENT 
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124  APTR(NAC)  =  NENT 
C 


N  =  N  TET( 31 
DU  131  NENT*  1  »N 

IF  (IA3-  PGWNAMI  NENT ) )  131t  132.  131 

131  CONTINUE 
STOP  1 

132  DO  13C  L  =  1*N 

IF  (PCALCIL))  12ft.  12R,  127 

127  IF  (PCALC(L)  -  NENT)  13C,  129,  130 
13P  CONTINUF 

128  PCALC(L)  =  NENT 

129  PPTP (NAC )  =  NENT 
C 

N  =  NTET(4) 

CO  135  N  ENT*  1 ,  N 

IF  {  I A4  -  PNLNAMI  NENT  1  )  1  35,  1  36  ,  135 

135  CONTINUE 
STOP  1 

136  EPTR(NAC)  *  NENT 
C 

15C  CONTINUE 

IF  (ERPFLG)  155,  155,  152 

152  IF  (NAC)  2,  2,  153 

153  RETURN 

155  WRITE  (3,  3913) 

GO  TO  2 

4010  FORMATI'-V,  2015  /  (IX,  2015)) 

PETUPN 


C 

C 


3000  FORMAT  (9X,  - - 

l - 1  ) 

3001  FORMAT  (MS  15X,  13A4) 

3002  F0PMA1  (18X,  'START  AT',  F8.0,  15X,  'RUNWAY  LENGTH',  F8.0  / 

1  19X,  'STOP  AT',  F3.0,  17X,  'GLIDE  SLOPE',  F8.1  / 

2  1  8X ,  'INTERVAL',  F8.0) 

3003  FORMAT  ('O',  12X,  ' A  I RC RAF T  C LA  SSI F I C ATI  ON* ,  18X,  'VOLUME  OF  OPERA 

IT  IONS'  /  13X,  - - -  TRIP  LENGTH 

2  - •  /  14X,  'NAME  DESCRIPTION',  6X  ,  'CATEGORY', 

3  7X ,  '  DA  YT  IMF  NIGHTTIME'  /  1 4  X,  • -  - -  ---- 

3011  FORMAT  ('O',  6X,  110,  6X,  3A4 ,  110,  6X,  2FU.3  /  (35X,  110,  6X, 

1  2  F  1 1  •  3  )  ) 


3901  FORMAT  ( 1 9  X, 
3  y,r?  FORMAT  (19X, 
3"ni  FORMAT  ( 19X, 
39-7  FORMAT  ( 19X, 
FORMAT  (19X, 
39''5  FORMAT  (19X, 

3'i'9i  Format  ( •  o • , 

39  19  FORMAT  ('O', 

i  (3sx,  in, 

3110  FORMAT  (  1 9 X , 


•*  STARTING  POINT  CANNOT  BE  NEGATIVE') 

•*  STOPPING  POINT  CANNOT  RE  NEGATIVE') 

•*  INTERVAL  IS  OF  WRONG  SIGN') 

»*  GLIDE  SLOPE  CANNOT  BE  NEGATIVE') 

•*  RUNWAY  LENGTH  MUST  BE  GREATER  THAN  ZERO'  ) 

'*  RUNWAY  TOO  SHORT  FOR  GLIDE  SLOPE') 

18X,  •*  INVALID  OPFRATION  CODEC,  110,  ')♦) 

6 X ,  110,  5 X,  «*  UNUEF1NE0  *',  IlO,  6X ,  2F11.3  / 
6X ,  2FU.31) 

'*  altitude  profile  undefined  p»  in,  •  i • ) 
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3711  1-UP  MAT  (  i<7X,  •*  JFLTA-EPNL  PROF  lit  UNDEFINED  (*,  CIO,  *1*1 

3912  FORMAT  (19X,  •*  tPNL  FUNCTlUN  UNDEFINED  IS  11%  ')•) 

3913  FORMAT  I 4X,  '*+*  PROGRAM  CANNOT  COMPUTE  THIS  CONTOUR  *M'| 

3914  FORMAT  (‘OS  13X,  •*  CAPO  MISSING  UR  OUT  OF  QRuER,  A/C  CLASS*  ,  I 

1  *  ,  CAP  D  NO  1 «  I  3 J 

3915  FOPMAT  CO*,  13X,  •*  CONDENSED  AIRCRAFT  TAbLE  FULL') 

C 

10C1  FORMAT  (12AS  A3,  2F7.D,  F5.C,  F6.0,  F*t.n) 

1C 02  FORMAT  (213,  8F8.D 
END 
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SUBROUTINE  POATA 
PURPOSE 

TO  ENTER  AIRCRAFT  DESCRIPTORS  AND  PERFORMANCE  DATA 
USAGE 

CALL  PDATA  (OPCODE) 

DESCRIPTION  OF  PARAMETERS 

OPCODE  -  INTEGER  VALUE  TO  OENCTE  ACTION  TO  BE  TAKEN 

1  =  ENTER  NEW  AIRCRAFT  DESCRIPTOR 

2  =  ENTER  NEW  ALTITUDE  PROFILE 

3  =  ENTER  NEW  OELTA  EPNL  PROFILE 

4  =  ENTER  NEW  EPNL  PROFILE 

11  -  FIX  AIRCRAFT  DESCRIPTOR  TABLE 

12  *  FIX  ALTITUDE  TABLE 

13  =  FIX  DELTA  EPNL  TABLE 

14  =  FIX  EPNL  TABLE 

15  =  FIX  ALL  TABLES 

21  =  EXPUNGE  AIRCRAFT  DESCRIPTOR  TABLE 

22  =  EXPUNGE  ALTITUDE  TABLE 

23  =  EXPUNGE  OELTA  EPNL  TABLE 

24  =  E XPUNGE  EPNL  TABLE 

25  -=  EXPUNGE  ALL  TABLES 

REMARKS 

1.  EXPUNGE  COMMAND  DELETES  ONLY  THOSE  TABLE  ENTRIES  NOT 
FIXED. 

2.  ERROR  CONDITIONS 

A.  INVALID  OPCODE  -  NO  ACTION  TAKEN.  MESSAGE  IS 

PRINTED  ALONG  WITH  OPCODE  VALUE.  ROUTINE 
RETURNS  WITH  OPCODE  =  0. 

B.  INVALID  DATA  -  NO  ACTION  TAKEN.  MESSAGE  IS 

PRINTEO  ALONG  WITH  OFFENDING  OATA.  ROUTINE 
RETURNS  WITH  OPCODE  =  0. 

C.  AVAILABLE  TABLE  AREA  OVERFLOW  -  NO  ACTION  TAKEN. 

MESSAGE  IS  PRINTED.  ROUTINE  RETURNS  WITH 
OPCODE  =  0. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  CALLED 
NONE 


SUBROUTINE  PDATA  (OPCODE) 

COMMON  X  STAR  T »  XEND »  DX,  NAC ,  CONTUR,  NMAX(4) 

COMMON  NTET(  4 ) *  NF  E  T ( 4 ) »  NALT( 12 ) »  NP0W(12),  ACALf(l2),  PCALC(12> 
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COM  l 
COM  2 


C-37 


PDATA 


COK  3 
COF  4 
Cl>  3 

coy  h 


NSEG  =  8 


IF 

(GPCOOt 

- 

1  ) 

900,  11C,  10 

10 

IF 

(OPCODE 

- 

3) 

120,  13C,  11 

11 

IF 

(OPCODE 

- 

3) 

140,  90C ,  12 

12 

IF 

(OPCODE 

— 

11) 

900,  210,  13 

13 

IF 

(OPCODE 

- 

13) 

220,  230,  u 

14 

IF 

(OPCODE 

- 

151 

240*,  2  50,  15 

19 

IF 

(OPCODE 

_ 

21) 

900,  310,  16 

16 

IF 

(OPCODE 

- 

23 ) 

320,  330,  17 

17 

IF 

(OPCODE 

- 

25) 

340,  350,  18 

lfi 

GO 

TO  TO" 

«OC  WRITF  (3,  3900)  OPCODE 


110  READ  II,  I0C1)  NAME  *  DTEMP,  (  ATEMP(  I  )  ,  PTEMPI  I  )  »ETEMP(  I )  •  1*1, NSEG) 

M  *  N  T  E  T  (  l  ) 

N  =  M  ♦  1 

IF  ( N-NMA XI  II)  111,  111  ,  910 

111  CO  112  1*1, M 

IF  (DSCNAM(I)  -  NAME)  112,  911,  112 

112  CONTIMIF 
CSCNAMN)  *  NAME 
i:n  1 1 a  i*i,3 

1  14  DSCPPTI  I,N)  *  DTEMPI  I  ) 

PO  116  I*  1,N  SEG 
ALTCIM  I,N)  *  ATEMPI  I  ) 

POWOIP  I  I , N )  *  PT> MP(  I  I 
lit  FPNiHM  I,N)  *  E TEMPI  I  ) 

MET!  1  )  *  N 
RETURN 

i !  C  Vf'lTE  (3,  34  ID 

w.  P  I  T  I  I  it  39'  l )  NAML 

goto 

d  i  uw  i  t r  r»,  Tom 

V.PITC  (3,  3CM)  NAME 
cr»  t«i  999 


1  2*  NN  *  • 

(•CAD  (1,  1003)  NAME,  NP,  (  I  AL  TTE  M I  J  ,  I  )  ,  1-1,2),  J*l»6) 

FIGURE  C-13B.  PROGRAM  LISTING  -  SUBROUTINE  PDATA 


COMMON  0SCNAMI20),  ALTNAm(12),  POWNAMI 12),  PNL  NAM ( 25 ) ,  APT  P  ( 50  ) 
COMMON  PPTR ( 50) ,  EPTtUSD,  0P6R150),  ZI12),  POUSPLI  12),  I0ENTI13) 
COMMON  0SCRPTC3.2  3),  A  L  ID  I  H  (  ft ,  2D  ,  POWDI  P  1 8  ,20  I  ,  EPNIH  R(3  •  20) 
COMMON  PNLEEFI  .3*,,  2,25)  ,  ALT!  ID?  ,12)  ,  POfcE  R  1 10  ,2,12) 

1NTEGEP  OPCODE,  0 SC. NAM ,  USCRPT,  ALTNAM,  POWNAM,  PNLNAM 
INTEGER  PRP »  ALTOIk,  PCfcDIR,  EPNDIP,  DTFMP,  ATEPP,  PTEMP,  ETEMP 
DIMENSION  DTEMPI 3),  ATEMPI8),  PTEMPI 8 ) ,  ETFMP(p) 

LIMENSION  ALTTEMI  10,2)  ,  PUWTEMH  0,2  )  ,  E  PNTEM  {  35 , 2 ) ,  TEMPIIOI 
EQUIVALENCE  (  ALTTF  MI  1  ,  l  )  ,  POWTEMIl,!),  EPNTEMI  1,1)) 


C-38 


POATA 


121  NN  =  NN  ♦  6 

If  (NP  -  NN J  123,  123,  122 

122  P  HAD  U,  1004)  ( (  ALTTEMI  J  ,1 )  ,  1*1,2),  J*7,10) 
GO  TO  121 

123  M  =  N  T  E  T  ( 2) 

N  =  M  ♦  1 

IF  (NP)  922,  922,  124 

124  IF  INP  -  10)  125,  125,  922 

125  IF  (N  -  NMAXI  2  > )  126,  126,  923 

126  00  127  1*1, M 

IF  ( ALTNAMI I )  -  NAME)  127,  924,  127 

127  CONTINUE 

DO  128  1*1,2 

00  128  J= 1, NP 

128  ALT ( J , I ,N )  *  ALTTE  MI J , I  ) 

NALT(N)  ■  NP 

ALTNAMi N  )  =  NAME 
NTETI2)  =  N 
RETURN 


922 

WRITE 

<3, 

3922) 

WRITE 

(  3, 

390  1) 

NAME 

GO  TO 

999 

923 

WRITE 

(  3, 

392  3) 

WRITE 

<  3, 

3901) 

name 

GOTO 

999 

924 

WRITE 

(3, 

3924) 

WRITE 

(3, 

390  1 ) 

NAME 

GO  TO  999 
C 
C 

130  NN  «  0 

READ  (  1,  1903)  NAME,  NP  ,  (  (  POWTEM ( J  ,  I  )  t  1=1,2),  J=l,6) 

131  NN  *  NN  ♦  6 

IF  (NP  -  «N)  133,  133,  132 

132  READ  (  1.  1004)  ((  POWTEM <  J  ,  l )  ,  1=1,21,  J=7,10) 

GO  TO  131 

133  M  *  NTETI3) 

N  =  M  ♦  l 


IF  (NP) 

932,  932, 

134 

134 

IF  (NP  - 

10)  135, 

1 35 1 

,  932 

135 

IF  (N  - 

NMAX(  3) ) 

136, 

136, 

933 

136 

DO  137 

I 

3  i, M 

IF  (POWNAMI  I)  -  NAME  ) 

137, 

934 

137 

CONTINUE 

DO  138 

I 

=  1,2 

CO  138 

J 

3  1  ,NP 

138 

POWER  ( J, 

I 

,N )  3  POWTEMI 

J,I  ) 

NPOW(N) 

« 

NP 

PUWNAM ( N 

) 

=  NAME 

NTETI 3) 

= 

N 

RETURN 

9  32 

WRITE  (j 

1 , 

3932) 
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WRITE 

(3, 

39011 

NAME 

GO  TO 

999 

933 

WRITE 

(3, 

3933) 

WRITE 

I  3, 

3901) 

NAME 

GO  TO 

999 

934 

WRITE 

(3, 

3934) 

WRITE 

(3, 

3901) 

NAME 

GO  TO 

999 

C 

c 

140  ERR  FL  G  =  1 

DO  155  L*l,8 

READ  (lf  1002)  NAME,  PRP»  Ll»  L2 ,  TEMP 
IF  U  -  1)  141,  141,  142 

141  NAM  El  =  NAME 
C 

142  IF  ( NAME1  -  NAME)  942,  143,  942 

942  WRITE  (3,  3942) 

ERRFLG  *  0 

143  IF  (PRPI  943,  943,  144 

144  IF  IPRP  -  2)  145,  145,  943 

943  WRITE  (3,  3943) 

ERR  FLO  ■  0 

145  IF  (L  'i )  945,  545,  146 

146  IF  (L 2  -  35)  147,  147,  945 

147  IF  (L2  -  LI  -  101  148,  148,  945 

148  IF  (L2  -  LI)  945,  945,  150 

945  WRITE  t  3,  3945) 

ERRFLG  *  0 

15C  IF  (ERRFLG)  155,  155,  1  51 

151  11=0 

00  152  1*L1,L2 

11=11+1 

152  EPNTEM(  I,PRP)  =  TEMPU1) 

155  CONTINUE 

IF  (EFRFLG)  949,  949,  156 
C 

156  M  =  NTET(  4) 

N  *  M  ♦  1 

IF  (N  -  NMAX(  4 )  )  157,  157,  947 

157  CO  158  1*1, M 

IF  (PNLNAM(I)  -  NAMED  158,  948,  158 

158  CONTINUE 
C 

160  00  162  1*1,35 

DO  162  J=  l, 2 

162  PNLEFF(  I,  J,N)  *  EPNTEM(I,J) 

PNLNAM(N)  *  NAME  1 
NTET(  4  )  *  N 
RETURN 

947  WRITE  (3,  3947) 

WPITfc  (  3,  3901)  NAME  1 
GO  TO  999 

948  WRITE  (  3,  3948) 

WRITE  <  3,  3901)  NAME  l 

FIGURE  C-14D.  PROGRAM  LISTING  -  SUBROUTINE  PD4&A 
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GO  TO  999 

949  WRITE  <3,  3901)  NAME  1 
GO  TO  999 

999  OPCODE  *  0 
RETURN 
C 
C 

210  NFET(  1)  *•  NFETl  1) 

RETURN 

C 

220  NFET(  2)  =  NTE  T(  21 
RETURN 
C 

230  NFET( 3)  =  NTE T ( 3) 

RETURN 

C 

240  N  FF T( 4 )  =  NTF  T( 4) 

RETURN 

C 

250  00  252  1=1,4 

252  NFETl  I)  =  NTE T(  l  ) 

RETURN 

C 

C 

3 1C  NTET(l)  =  NFETl  II 
RETURN 
C 

320  NTET( 2)  =  NFETl 2) 

RETURN 

C 

33"  NTETC31  =  NFFTl 3) 

RETURN 

C 

340  NTEH4)  =  NFE  T(  4 ) 

RETURN 

C 

350  DO  352  1*  1,4 

352  NTET(  I)  =  NFETl 1 1 
RETURN 
C 
C 

1001  FORMAT  112,  IX,  3A  4 ,  IX,  9(312,2X1) 

1002  FORMAT  (212,  213,  10F7.3) 

1C"3  FORMAT  (213,  2X,  12F6.0) 

1004  FORMAT  <8X,  12F6.0) 

3900  FORMAT  ('O',  IPX,  •*  INVALID  OPER ATI  ON  C ODE  (',  14,  •  I  •  I 

3901  FORMAT  (21X,  'CANNOT  ENTER  NAME  =  ',  14) 

3910  FORMAT  ('O',  18X,  •*  AIRCRAFT  DESCRIPTOR  TABLE  FULL') 

3911  FORMAT  PC',  18X,  •*  OUPLICATE  AIRCRAFT  DESCRIPTOR  NAME' » 

3922  FORMAT  PC',  1BX,  •*  INVALID  NUMBER  OF  ALTITUDE  POINTS') 

3923  FORMAT  PC',  18X,  »*  ALTITUDE  PROFILE  TABLE  FULL') 

3924  FORMAT  ('O',  18X,  DUPLICATE  ALTITUDE  PROFILE  NAME') 

393?  FORMAT  PC',  18X,  •*  INVALIO  NUMBER  OF  OELT  A-EPNL  POINTS') 

3)33  FORMAT  ('O',  18X,  •*  OELTA-EPNL  PROFILE  TABLE  FULL') 

3934  FORMAT  ('O',  I8X,  •*  DUPLICATE  OELTA-EPNL  PROFILE  NAME') 
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PDATA 


3942 

FORMAT 

( *0*, 

3943 

FORMAT 

I  'O', 

3945 

FORMAT 

1*0*, 

3947 

FORMAT 

(  •  0  • , 

3948 

FORMAT 

(  'O', 

ENO 


18X,  •*  PNLEFF  PROFILE  NAMES  DO  NUT  MATCH*) 
18X,  **  INVLAID  PNLEFF  PROPAGATION  PATH*) 
18X,  •*  INVALID  PNLEFF  RANK  LIMITS*) 

18X,  PNLEFF  PROFILE  TABLE  FULL*) 

18X,  •*  DUPLICATE  PNLEFF  PROFILE  NAME') 
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SUBROUTINE  LAST°T 
PURPOSE 

TO  COMPUTE  THE  POINT  ON  THE  FLIGHT  TRACK  WHERE  THE 
CONTOUR  CLOSES 

USAGE 

CALL  LASTPT  (X,  XCIOSE) 

DESCRIPTION  OF  PARAMETERS 

X  -  THE  CALLING  VALUE  OF  THIS  PARAMETER  IS  A  POINT 

ON  THE  FLIGHT  TRACK  NEAR  WHICH  THE  CONTOJR  CLOSES  . 

XCLOSE  -  THE  RETURNED  VALUE  OF  THIS  PARAMETER  IS  THE  POINT 
ON  THE  FLIGHT  PATH  WHERE  THE  CONTOUR  CLOSES 
I  IE,  THE  NEF  VALUE  AT  THIS  POINT  IS  EQUAL  TO  THE 
NEF  VALUE  OF  THE  CONTOUR  CURRENTLY  BEING 
COMPUTED! • 

REMARKS 

1.  A  VALUE  OF  ZERO  WILL  BE  RETURNED  IF  THE  ROUTINE 
CANNOT  COMPUTE  A  CLOSING  POINT. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
ANEF,  PARAM,  XNEW 

METHOD 

THE  ROtriNE  PICKS  TWO  POINTS  ON  THE  FLIGHT  TRACK  AND 
COMPUTES  THE  NEF  VALUES  AT  THESE  POINTS.  INTERPOLATION 
BETWEEN  THESE  POINTS  USING  THE  VALUE  OF  THE  CURRENT  CONTOUR 
BEING  COMPUTED  YIELDS  THE  CONTOUR  CLOSING  VALUE  OF  X  TO 
BE  RETURNED. 

ONE  OF  THE  POINTS  CHOSEN  IS  THE  CALLING  VALUE  OF  X.  IF 
THE  ALTITUOE  OF  ANY  AIRCRAFT  IS  LESS  THAN  12.5  FEET  AT 
THIS  POINT  THE  ROUTINE  TERMINATES  AND  RETURNS  WITH  X 
EQUAL  TO  ZERO. 

THE  SECOND  POINT  IS  VARI4BLE  AND  IS  CHOSEN  IN  THE  FOLLOWING 
MANNER.  A  POINT  ON  THE  TRACK  ONE  COMPUTATION  INTERVAL 
CDX!  LESS  THAN  THE  CURRENT  VALUE  OF  X  IS  CHOSEN.  IF  ANY 
ALTITUOFS  APE  LESS  THAN  12.5  FEET  AT  THIS  POINT,  ANOTHER 
POINT  IS  CHOSEN  HALF  AGAIN  THE  DISTANCE  BETWEEN  THESE  TWO 
POINTS.  THE  ALTITUDE  TEST  IS  PERFORMED  AGAIN  AND  IF 
ANY  ALTITUDE  T5  LESS  THAN  12.5  FEFT  THE  DISTANCE  IS 
HALVED  ONCE  AGAIN.  THIS  SEQUENCE  IS  REPEATED  JNTIl 
A  USABLE  VALUE  OF  X  IS  OBTAINED.  IF,  AFTER  10  TRIALS, 

NO  USABLE  VALUE  OF  X  ! S  OBTAINED  THE  ROUTINE  TERMINATES 
ANO  RETURNS  WITH  X  EQUAL  TO  ZERO. 
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SUBROUTINE  USTPT  (X,  XCLOSE) 

C 

COMMON  X  STAR  T ,  XENO,  OX,  NAC ,  CONTUR,  NMAXI4)  COM  1 

COMMON  NTETI  4) ,  NFE  T(  4) ,  NALTU2),  NPOWI12),  ACALCI12),  PCALCC12)  COM  2 
COMMON  0  SCNAM(  20 ) ,  ALTNAMU2),  POWNAM( 12) ,  PNLNAMI25),  APTR(50)  COM  3 
COMMON  PPTR  (  50) ,  EPTR(  50)  ,  OPERI50),  Z<12),  POWSPLI l 2)  ,  IDENT<13)  COM  4 
COMMON  DSCRPTI 3,20),  ALTDI RI 8,20) ,  POWDI R(8  ,20 ) ,  EPNDIR(8,20)  COM  5 

COMMON  PNLEFF ( 35, 2 » 2 5)  ,  ALT ( 1 0,2 , 12 ) ,  POWER! 10 ,2 ,12)  COM  6 


C 

INTEGER  APU 


C 

C  -  CHECK  AIRCRAFT  ALTITUOES  AT  CALLING  VALUE  OF  X,  IF  ANY 

C  ARE  LESS  THAN  12.5  FEET  TERMINATE  COMPUTATIONS  FOR 

C  THIS  CONTOUR  - 

C 

C 


DO  5  1*1,  NAC 
J  =  APTR (  I) 

IF  (Z(J)  -  12.5)  30,  5,  5 
5  CONTINUE 


C 

c  -  IF  ALTITUDES  ARE  OK  COMPUTE  NEF  AT  THIS  POINT  - 

C 

CALL  ANEF  (0.,  NEF  2 ) 

C 

c  -  N0W  COMPUTE  NEF  DIRECTLY  UNOER  PATH  AT  SOME  POINT 

C  WITH  A  LESSER  VALUE  OF  X - 

C 

XLESS  *  X  -  DX 
XINT  ■  DX 


C 

CO  10  IT  *  1,10 
C 


C  - CALCULATE  PERFORMANCE  PARAMETERS - 

C 

CALL  APO  (XLESS) 

C 

C  - CHECK  FOR  ANY  ALTITUOES  LESS  THAN  12.5  FEET 

C 

DO  7  1*1, NAC 


J  =  AP  TP  (  I  ) 

IF  IZ(J)  -  12.5)  8,  7,  7 
7  CONTINUE 


C 

C  -  COMPUTE  NEF  AT  THIS  POINT  - 

C 

CALL  ANEF  <0..  NEF  1 ) 

C 

C  -  INTERPOLATE  TO  FIND  VALUE  OF  X  WHERE  CONTOUR  CLOSES 


C 

XCLOSE  *  XNEW  (XLESS,  NEF  I  ,  X,  NEF2,  CONTUR) 

IF  (ABSIX  -  XCLOSt)  -  ABS(OXI)  40,  40,  30 

C 

c 

f.  - TRY  ANOTHER  VLAUE  OF  X - 

FIGURE  C-I4B.  PROGRAM  LISTING  -  SUBROUTINE  LASTPT 
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LASTPT 


C 

8  XINT  =  XINT/2 . 

XLESS  *  XLESS  ♦  XINT 
ID  CONTINUE 


- ROUTINE  CANNOT  COMPUTE  CLOSING  VALUE  OF  X - 

30  X CLOSE  *  0. 


- RESTORE  ALT  I TUOE  S  AND  DELTA-EPNL 'S 

40  CALL  APDIXJ 
PETURN 
C 

END 


FIGURE  C-UC.  PROGRAM  LISTING  -  SUBROUTINE  LASTPT 
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SUBROUTINE  ANEF 


PURPOSE 

TO  COMPUTE  THE  NEE  VALUE  AT  A  POINT  P  ,  A  GIVEN 
DISTANCE  PERPENDICULAR  TO  THE  FLIGHT  TRACK. 

USAGE 

CALL  ANEF  <Y,  NEFVAU 
DESCRIPTION  OF  PARAMETERS 

V  -  PERPENDICULAR  DISTANCE  FROM  POINT  TO  FLIGHT  TRACK 
NEFVAL  -  NEF  VALUE  AT  POINT  P 

REMARKS 

1.  IT  IS  ASSUMED  THAT 

A.  NAC  IS  GREATER  THAN  ZERO 

B.  ALT  I TUOE  S  AND  EPNL  CORRECTIONS  HAVE  PREVIOUSLY 
BEEN  TABULATED. 

C.  NO  ALT  I TUOE  IS  EQUAL  TO  ZERO  IF  THE  CALLING 
VALUE  OF  Y  IS  EQUAL  TO  ZERO. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
EPNL,  SQRT,  ALOGIC 


SUBROUTINE  ANEF  (Y,  NEFVAL) 


COMMON  X  STAR  T ,  XEND,  DX,  NAC,  CONTUR,  NMAXI4)  CUM  I 

COMMON  NTETI4),  NFET14),  NALTII2I,  NPOWI12),  ACALCU2I,  PCALCI 12)  COM  2 
COMMON  QSCNAMI20),  ALTNAM(  12 1  »  POWNAMU2),  PNLNAMI25 )  ,  APTRI50)  CUM  3 
COMMON  PPTP ( 50 ) ,  EP  TR ( 50) ,  0PERI50I,  ZI12),  PQWSPLI 12) ,  I  DENT ( 131  COM  A 
COMMON  DSCRPTI 3,20),  A  LTD  I R I  8, 20) ,  POWDI RI 8 ,?0 )  ,  EPNDIR(8,20)  COM  5 

COMMON  PNLEFFI  35,  2,25) ,  ALT! 10,2 ,12)  ,  POWE  RI  10  ,2 , 12  >  COM  6 


INTEGER  APTR  ,  PP TP,  EPTR,  PTR 

PEAL  NEFLIN,  NEF,  NEFVAL 

NEFL  IN  =  0. 

-  COMPUTE  NCF  FOR  EACH  AIRCRAFT  AND  COLLECT  CUMULATIVE  SUM 

DO  20  1=1, NAC 

-  PICK  UP  POINTER  TO  ALTITUDE  FOR  THIS  AIRCRAFT  ANO 

CALCULATE  SLANT  DISTANCE - 

PTR  =  APTR  I  I  ) 

DIST  =  SQkT(Y*Y  ♦  L ( P TR  )*Z ( PTR) ) 

C  - CALCULATE  PROPAGATION  FACTOR - 
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A  NEF 

C  1.  If  LESS  THAN  OR  EQUAL  TO  ZERO  USE  GROUND  TO 

C  GROUND  PROPAGATION. 

C  2.  IF  BETWEEN  ZERO  AND  ONE  INTERPOLATE  BETWEEN 

C  GRUUNU  TO  GPQUNO  AND  AIR  TO  GROUND  PROPAGATION. 

C  3.  IE  GREATER  THAN  OR  EQUAL  TO  ONE  USE  AIR  TO 

C  GROUND  PROPAGATION 

C 

C  - NOTE 

C  PROPAGATION  FACTOR  IS  CALCULATED  ON  THE  BASIS  Of  THE 

C  ANGLE  Of  ELEVATION  Of  THE  AIRCRAFT  WITH  RESPECT  TU 

C  THE  OBSERVER  AND  THE  GROUND  PLANE.  PF  =  0  AT 

C  ARCS  INI  AND  PF  =  1  AT  ARCS  I  N<  0 . 1 25  > , 

C  4  DEG  la  MIN  AND  7  DEG  11  MIN  RESPECTIVELY. 

C 

PF  =  <Z(?TK)  /  0 1  ST  -  0.075)  /  0.  05 
C 

C  -  PICK  UP  POINTER  TO  PNLEFF  LIST  FCR  THIS  AIRCRAFT  - 

C 

P  TP  =  EPTRI  H 
C 

IF  IP  F )  5,  5,  3 
3  IF  (PF  -  I.)  fc,  7,  1 
C 

C  - GET  FPNL - 

C 

5  P  =  EPNLIPTR,  I,  0 l ST ) 

GO  TO  10 

C 

6  P  =  (1.  -  PF)*£PNLtPTP, l.DIST)  +  PF  *E  PNL  I  PT  R ,  2 , 01  ST  > 

GO  TO  K 
C 

7  P  =  EPNLIPTR, 2,0 l ST) 

C 

c  - PICK  UP  POINTER  TO  E  PNL  CORRECTION  FOR  THIS  AIRCRAFT 

C  AND  COMPUTE  TOTAL  EPNL  FOR  THIS  AIRCRAFT  - 

C 

10  PTR  =  PPTRII) 

P  =  P  +  POWSPLI  PTR  I 

C 

C  -  TRUNCATE  EPNL  IF  LESS  THAN  80 - 

C 

NEF  =  o. 

IF  If  -  75.)  15,  19,  13 

13  IF  IP  -  8C . )  14,  14,  15 

14  P  =  (80.  ♦  OPERim  *  IP  -  75.)  *  0.2 
C 

C  -  COMPUTE  NEE  FUR  THIS  AIRCRAFT  - 

C 

1 5  NEF  =  P  *  OPFRI  I  ) 

C 

C  -  CONVERT  NEF  FRUM  LOG  TO  LINEAR  AND  ADD  TO  SUM  - 

C 

10  r  EFL IN  =  NEELIN  ♦  10. ** I NEF /I 0. ) 

2?  CONTINUE 

C 

c  -  CONVERT  SUM  FR  CM  LINEAR  TO  LOG - 

C 

NFFVAL  *  n.  *  AL0G10I  NEF  LI  N)  -  75. 

C 

RETURN 

C 

END 

FIGURE  C-15B.  PROGRAM  LISTING  -  SUBROUTINE  A  NEF 
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C 

C  . . . . . . 

c 

c  subroutine  apd 

c 

PURPOSE 

TO  COMPUTE  THE  NECESSARY  ALTITUDES  AND  DELT A-EPNLS 
AT  A  POINT  ALONG  THE  FLIGHT  TRACK. 

USAGE 

CALL  APD  IX) 

DESCRIPTION  OF  PARAMETERS 

X  -  DISTANCE  ALONG  FLIGHT  TRACK  FROM  START  OF 
TAKEOFF  ROLL. 

REMARKS 

I.  ROUTINE  ASSUMES  THAN  TWO  PROFILE  POINTER  TABLES 
(ACALC  AND  PCALC)  HAVE  BEEN  BUILT  CONTAINING  THE 
PROFILE  ENTRIES  FROM  WHICH  TO  CALCUf.ATE  THF 
PARAMETER  VALUES. 

?.  ROUTINE  ASSUMES  THAT  SPECIFIED  PROFILES  HAVE 
BEEN  DEE  INFO. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
CURVE 

METHOD 

ROUTINE  PICKS  UP  PROFILE  ENTRY  FROM  AFOREMENTIONED 
PREVIOUSLY  BUILT  TABLES  AND  COMPUTES  THE  PARAMETER  VALUE 
FROM  THE  SPECIFIED  PROFILE  FOR  THE  SIVEN  VALUE  OF  X. 
PARAMETER  VALUE  IS  THEN  PLACEO  IN  TABLE  WHICH  IS  PARALLEL  TO 
THE  PROFILE  POINTER  TABLES. 


SUBROUTINE  APD  IX) 


COMMON  X  START,  XEND,  OX,  NAC  ,  CQNTUR,  NMAX  14)  COM  1 

COMMON  NTETI4),  NFETI4),  NALTI  12  )  ,  NPCIW(12),  ACALC  (12),  PCALC  (12)  COM  2 
COMMON  DSCNAM( 2D  )  ,  ALTNAM( 12 ) ,  POWNAMU2),  PNLNAM(25),  APT  P ( 50  )  COM  3 
COMMON  PPTRfBP),  EPTMSP),  0PER(50>,  1(12),  PUWSPL(I2),  IDENU13)  COM  4 
COMMON  DSCRPT( 3,20) v  ALTDI R ( 8 ,20) ,  PCWDI R( 8 ,20 ) ,  EPNDIR(8,20)  COM  5 

COMMON  PNltFFI 35, 2,25)  ,  ALT( 10 , 2 , 1 2 )  ,  PUWE R < 1 0 , 2 1 1 2 )  COM  6 


DIMENSION  DUM Y (  2,  I n ) 

INTEGER  ACALC,  PCALC 

- COMPOTE  ALTITUDES - 

M  =  NTETJ  2) 

DO  10  1*1, rf 

C  - PICK  UP  PROFILE  NUMBER 


FIGURE  C-16A.  PROGRAM  LISTING  -  SUBROUTINE  APD 
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J  -  ACALC (  I  > 

TEST  FUR  ZERO  lit.  END  OF  STRING )  - 

IF  (J)  ?u,  2",  1 5 

PICK  UP  NlJHftfcR  OF  POINTS  DEFINING  THIS  PROFILE 
15  N  =  NALT(  J  ) 

i-UAo  this  profile  into  dummy  array _ --- 

DO  10  K=  1 , 2 
DO  IP  L=  1 » N 

1 8  CUM Y l K t L I  *  ALTIL fR, J) 


- CALCULATE  PARAMETER  VALUE 

lq  )  *  CURVE!  X »DUM  Y , N ) 


COMPUTE  EPNL  ADJUSTMENTS 


M  =  NTET! 3) 

20  DO  29  1=1  ,M 
J  =  PCAl  Cl  I) 

IF  (J)  30,  30,  2  5 
25  N  =  NP0W! J ) 

DO  ?8  K=  1 , 2 
DO  28  L=1,N 

28  CUMYI  K,L  )  =  POWER(L,K,J| 

29  POWSPL(J)  =  CURVE!  X,DUMY,N) 

30  RETURN 


FIGURE  C-16B,  PROGRAM  LISTING  -  SUBROUTINE  APD 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTION  £  PNL 
PURPOSE 

TO  COMPUTE  AN  EFFECTIVE  PERCEIVED  NOI SF  LEVEL  AS 
FUNCTION  OF  DISTANCE  FR CM  AIRCRAFT  TO  OBSERVER. 


USAGE 

VAR  IABLE 


=  EPNULIST,  PRPGTN,  DIST) 


DESCRIPTION  OF  PARAMETERS 

LIST  -  INTEGER  TO  SPECIFY  WHICH  OF  18  PNL  VS.  DISTANCE 
FUNCTIONS  TO  USE  FROM  PNLEFF  TABLE  ( INTEGER  I 
PRPGTN  -  PROPAGATION  PATH  1  -  GROUND  TC  GROUND 
(  INTEGER)  2  -  AIR  TO  GROUND 

EPNL  -  EFFECTIVE  PERCEIVED  NOISE  LEVEL  (REAL) 

DIST  -  DISTANCE  FROM  AIRCRAFT  TO  OBSERVER  {REAL) 


REMARKS 
1 


FUNCTION  IS  DEFINED  FOR  ONLY  TWO  PROPAGATION  PATHS 
AND  FOR  DISTANCES  FROM  12.5  TO  31500  FEET.  IF  DIST 
IS  OUTSIDE  OF  DEFINED  RANGE  OR  PRPGTN  IS  NOT  EQUAL 
TO  1  OR  2,  AN  ERROR  MESSAGE  IS  PRINTEO  ALONG  WITH  THE 
VALUES  OF  THE  THREE  CALLING  PARAMETERS. 

EPNL  VALUES  RETURNEO  FOR  INVALID  DATA 


IF  DISTANCE  IS  LESS  THAN  12. 5  FEET  THE  EPNL 
AT  12.5  FEET  WILL  BE  RETURNED. 

IF  DISTANCE  IS  GREATER  THAN  31500  FEET 
AT  31  500  FEET  WILL  BE  RETURNED. 

IF  PRPGTN  IS  NOT  EQUAL  TO  1  OR  2 


THE  EPNL 


A  VALUE  OF 


SUBROUTINES 

ALOGIO 


ZERO  WILL  BE  RETURNEO 


AND  FUNCTION  SUBPROGRAMS  CALLED 


METHOO 

ROUTINE  FIRST  LOOKS  UP  TWO  PNLEFF  VALUES  IN  TABLE  FOR 
DISTANCES  BOUNDING  CALLING  DISTANCE.  INTERPOLATION 
BETWEEN  THESE  VALUES  IS  PERFORMED  TC  YIELD  DESIRED  EPNL, 


REAL  FUNCTION  EPNL ( L I  ST  , PRPGTN. DI ST) 

COMMON  XSTART,  XENO.  DX,  NAG.  CONTUR  ,  NMAXI4) 

COMMON  NTETI4),  NP£  T ( 4 ) »  NALTI12),  NP0WII2),  ACALC (12).  PCALC 1 12 ) 
COMMON  DSCNAMI 20 ) .  ALTNAMI 12 ) .  POWNAMI 12 ) ,  PNLNAMI 25  > .  APT  R 1 50 ) 
COMMON  PPTPI50),  EPTR  I  50  »  OPERI50),  ZU2),  P0WSPLU2),  I OENT  1 13 1 
COMMON  DSCRPTI 3,20),  A L TD I R ( 8,20)  ,  PQWDI R I  8 ,20 ) ,  EPND1R(8,20) 
COMMON  PNLEFF!  35,2,25) ,  ALT! 10,2  ,12  >  #  PUWE R 1 10 , 2 , 12 ) 


INTEGER  PRPGTN 


FIGURE  C-17A,  PROGRAM  LISTING  -  FUNCTION  EPNL 


COM  l 
COM  2 
COM  3 
COM  4 
COM  5 
COM  6 


C-50 


EPNL 


^al  10gd,  logoi 

tables  areMaseo  l--tHGARITHy[C  SCALC  uppn  which  pnleff 
IOG0  =  10.  *  AlOGlOfDIST)  -  I  c. 


CHECK  FOR  INVALID  PARAMETERS 


IF  (PRPGTNI  14,  14,  2 
2  IF  (PRPGTN  -  2}  ilf  llf  l4 

tl  IF  ( l  OGD  *•  l.|  15  it  |  j 

12  IF  n  oon  -  35.  )  ‘l3t  'l6,  ‘l6 

-  PARAMETERS  ARE  VALID  _ 

rauMCAre  fraction  part  of  log  distance 

13  I  LOGO  =  LOGO 
LOGO!  =  ILOGD 


COOK  UP  BOUNDING  PNLEFF  VALUES 


G1  =  PNLEFF! ILOGD, PRPGTN, LIST) 

C2  *  ^leff(ilogd,i,pRpg;n;l,st, 


INTERPOLATE 


Return  C1  *  ,C2*CU*<l-OGD-LOGDII 


parameters  are  invauo 


1A  EPNL  *  0. 

GO  TO  17 

c 

15  G0NT0S17NLEFM1,PPPGTN,UST> 

C 

c  16  tPNL  =  P^CEFF(  35, PRPGTN, LIST) 

17  muTpV3’  30nn  UST’  PRPGTN’  dist 

*  »  DIST  *»,  F9.0I 


FIGURE  C-178.  PROGRAM  LISTING  -  FUNCTION  EPNL 
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FUNCTION  XNE'W 


? UK POSE 

COMPUTE  A  VALUc  OF  X  CN  THE  CURVE  Y  =  M  *  LGG(X)  ♦  rt 
GIVEN  The  VALUE  OF  Y  FOP  WHICH  X  IS  TO  BE  COMPUTED  AND 
TWO  X  , Y  CnUPUlNATES  ON  THE  CURVE 


USAGE 

VARIABLE  =  XiMtWtXlt  Yl,  X2  ,  Y  2,  YNEW) 

DESCRIPTION  OF  PARAMETERS 

XI  -  VALUE  OF  X  AT  FIRST  COORDINATE  POINT 

Yl  -  VALUE  OF  Y  AT  FIRST  COORDINATE  POINT 

X2  -  VALUE  OF  X  AT  SECOND  COORDINATE  POINT 

Y2  -  VALUE  OF  Y  AT  SECOND  COORDINATE  POINT 

VNEW  -  VALUC  OF  Y  FOR  WHICH  X  IS  TO  BE  COMPUTED 
XNEW  -  VALUE  OF  X  AT  Y  =  YNEW 

REMARKS 

1.  IF  Yl  AND  'Y2  ARE  THE  SAME,  THEN  XNEW  IS  SET 
EQUAL  TO  X2. 

2.  IF  EITHER  XI  OR  X2  ARE  EQUAL  TO  OR  LESS  THAN  ZERO, 

then  xnew  is  set  wual  to  zero. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

MfcTHOD 

THE  EQUATION  USED  IS  GIVEN  IN  STATEMENT  30. 


REAL  FUNCTION  XNE w  I  XI,  Yl  ,  X2 ,  Y2 ,  YNEW) 
- CHFCK  FOR  INVALID  PARAMETERS - 


I  E 

(Yl  - 

Y  2 ) 

2, 

1C 

2 

IF 

(XI) 

?0, 

20, 

3 

3 

IF 

(X2) 

20, 

20, 

30 

-  IF  Yl  AND  Y 2  ARE  THE  SAME 

10  XNEW  =  X2 
RETURN 


r 

C  - IF  XI  AND  X 2  APE  LESS  THAN  OR  EQUAL  TO  ZERO 

C 

2C  XNEW  =  0, 

RETURN 

C 

C  -  IF  PARAMETERS  ARF  OK - 


30  XNEW  -  X2  *  (  X2/Xl)**( ( YNEW-Y2) /( Y2-Y1) ) 
RETURN 
C 

END 


FIGURE  C-18A.  PROGRAM  LISTING  -  FUNCTION  XNEW 


C-52 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


FUNCTION  OPCOR 


PURPOSE 

COMPUTE  THE  NEF  CORRECT  ION  FOR  A  GIVEN  NUMBER  OF 
DAYTIME! 0700  -  2200)  AND  NIGHTTIME  (2200  -  07001 
AIRCRAFT  OPERATIONS. 


USAGE 

VARIABLE  =  OPCORINDAY,  NNITE) 


DESCRIPTION  OF  PARAMETERS 

NDAY  -  NUMBER  OF  DAYTIME  OPERATION  (REAL! 
NNITE  -  NUMBER  OF  NIGHTTIME  OPERATIONS  (REAL) 
OPCOR  -  VALUE  OF  NEF  CORRECTION  (REAL) 


REMARKS 

1.  IF  THE  TOTAL  NUMBER  OF  OPERATIONS  (DAY  AND  NIGHT) 

IS  LESS  THAN  OR  EQUAL  TO  ZERO,  A  VALUE  OF  -100  WILL 
BE  RETURNED  AS  THE  NEF  CORRECTION. 


SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
AL0G1P 


METHOD 

THE  DEFINITION  OF  THE  NEF  CORRECTION  FOR  A  GIVEN  NUMBER 
OF  DAYTIME  AND  NIGHTTIME  OPERATIONS  IS  GIVEN  IN  STATEMENT  1. 


REAL  FUNCTION  OPCOR!NDAY,  NNITE) 

REAL  NDAY,  NNITE 

- CHECK  FOR  IN VA LI D  PARAMETERS - 

A  =  NDAY/20.  ♦  NNITE/1.2 
IF  (A)  5,  5,  l 

-  PARAMETERS  ARE  VALID.  CALCULATE  CORRECTION  - 

1  OPCOR  =  10.  *  ALOGIO(A) 

RETURN 

-  PARAMETERS  ARE  INVALID.  SET  CORRECTION  EQUAL  TO  -100 

5  OPCOR  =  -100. 

RETURN 

END 


FIGURE  C-19.  PROGRAM  LISTING 


FUNCTION  OPCOR 
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R (-  A L  FUNCT  ICN  CURV*1 


purpo  se 

rnvporFS  a  v  value  f^r  a  given  value  of  x,  w-m-rf 

V  =  FIX)  IS  DEFINED  FY  X  ,Y  POINTS  CCNNFCTFD  DY  '  r?Mi.tlT 
LINE  SEGMENTS.  (SEC  Rtmwj 

USAGE 

VNTIAliLF  =  C  IP  Vf  ( XL ,  Dim,  NPTS) 

Ut  SCR  I  PI  ION  Cf  p  AP.  ALTERS 

XL  -  VALUE  OF  X  F OH  WHICH  Y  IS  Tr  PE  Cf^PUTtO 

Pli^Y  (  l  » J)  -  T  A  °LE  Or  X,Y  POINTS  It  DEFINE  EU.CTIuN. 

I  =  1  FOR  X i  ?  FCR  Y 

NPTS  -  NIJMKFP.  OF  X,Y  POINTS  DEFINING  FUNCTION 

TEE' AR  K S 

1.  GUHY  ARRAY  PIJST  PE  CRCLRED  SUCH  THAT  VALUFS  OF  X 
ARE  CONSTANTLY  P'CRFASIVS 

?.  IF  NPTS  IS  GCUAL  TO  (JR  LFSS  I  HA  j  /FRO,  A  VALUr 
CF  ZERO  WILL  PF  R  ETLRN'FD 

A.  IF  THE  VALLE  OF  XL  DCFS  NOT  LIE  WITEIN  Thp  PFFINED 
RANGE  PF  THF  FUNCTION,  THFN  THE  VALLE  OF  Y  FOR  THE 
DEFINED  VALUE  CF  X  WHICH  IS  CLOSEST  TO  XL  WILL  8F 
RETLENFD.  (IE  TH c  FUNCTION  LOCKS  LIKE  T^IS.) 


X 


SUHR  OUT  I  NFS  A\X  FUNCTIONS  REQUIRED: 
NONE 


FIGURE  C-20A.  PROGRAM  LISTING  -  FUNCTION  CURVE 
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^r:Al  Filler  ICN  CURVF  (yi.  r,/wv  MTtl 
OJHMSIfiN  OIMVIMC,  ’  ’  NPTS) 

,  "  '^pTS  -  1)  4,  b,  | 

’  ,XL  "  ™*vn,j))  s,  s ,  ^ 

'  1  l=^f  NPI S  ' 

1  'Wu\y:  nu ”  U'J"  '■■  *.  ’ 

'»(J  TO  S 

«cpwii, ln!;vw;i,,?;irv,7,,,,*,XL  ■  ^“’‘n/ioiKYd.ji 

topn 

*i  CU3  Vf  =  OU^Y  (  2  i ) 

RE  TUM 
(>  CUR  Vf-  =  0 

h  m  rj 


FIGURE  C-206.  PROGRAM  LISTING  -  FUNCTION 


CURVE 
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APPENDIX  D 


COMPUTATION  OP  CONTOURS  FOR  PARALLEL  RUNWAYS 


APPENDIX  D 

COMPUTATION  OP  CONTOURS  FOR  PARALLEL  RUNWAYS 

The  NEP  program  computes  the  loci  of  points  of  equal  NEF 
value  resulting  from  operations  along  a  single  flight  path.. 
Although  it  is  possible  to  compute  the  NEP  contours  for  many 
different  flight  paths  in  sequence,  the  program  dees  not  take 
into  account  the  possible  influence  of  operations  on  other 
nearby  flight  paths. 

Larger  airports  in  particular,  are  likely  to  have  runways 
parallel  to  each  other.  Consequently  the  operations  on  one 
runway  may  influence  the  location  of.  the  contours  about  the 
other.  There  are  two  ways  tc  solve  this  problem.  One  is  to 
compute  the  NEP  contours  fer  :he  two  flight  paths  directly, 
using  a  modification  of  the  "normal"  NEF  program.  Although 
conceptually  very  simple,  this  approach  requires  considerable 
more  main  storage,  because  the  two  paths  may  well  have  different 
profiles  associated  with  them.  As  all  efforts  have  been  made 
to  keep  storage  requirements  to  a  minimum,  so  that  users  of 
small  computers  may  use  the  program,  expanding  storage 
requirements  for  this  purpose  seems  inconsistent. 


The  path  chosen  was  to  use  the  output  of  the  NEF  program  as 
input  lor  a  merging  program,  which  constructs  a  summed  contour 
from  the  digital  form  of  the  NEP  contour  output.  This  program 
is  completely  separate  from  the  NEF  program  although  it  needs 
the  output  of  the  NEF  routine  as  input.  In  this  report  the 
merging  program  has  been  named  MERGENEF.  Dependent  upon  the 
needs  of  the  user,  however,  the  program  can  be  Incorporated  in 
the  main  NEP  program.  Input  and  output  media  for  MERGENEF  as 
well  as  the  particular  way  in  which  it  may  be  called  will  be 
dependent  upon  the  operating  system  and  the  preference  of  the 
programmer. 

The  MERGENEF  program  will  accept  the  perpendicular  contour 
distances  to  the  flight  tracks  of  two  parallel,  arbitrarily 
displaced  tracks,  and  will  compute  the  proper  NEF  contour 
distances  for  the  operations  on  the  two  paths  together.  In 
order  to  perform  this  task  the  following  input  is  required: 

a.  Title  Card.  This  card  is  for  identification  and  can 
contain  up  to  80  characters  of  identifying  information. 


D-l 


TABLE  VIII 


MERGED  CONTOUR  CONTROL  CARD  #1  FORMAT 


Columns 

Contents 

Mode 

Range 

1-7 

Distance  along  flight  track 

Integer 

>0 

to  commence  merging 

ft 

8-14 

Increment  along  flight  track 

ti 

between  successive  computations 

tt 

15-18 

NEF  value  of  1st  merged  contour 

Real 

tt 

19-22 

"  tt  »  2nd  " 

It 

tt 

23-26 

f!  f!  ff  3pd  t* 

tl 

ft 

27-30 

"  "  "  4  th  !'  " 

It 

tt 

31-34 

"  "  "  5th  "  " 

tt 

tl 

35-38 

"  »  "  6th  "  " 

tl 

ft 

39-42 

H  ii  n  7  th  "  *• 

ft 

tl 

43-46 

«  n  it  gth  "  « 

tt 

ff 

47-50 

ii  H  it  9  th  "  *' 

tl 

tt 

51-54 

n  •,  10tL  i.  ii 

ft 

tt 

TABLE  IX 

MERGED  CONTOUR  CONTROL  CARD 

#2 

Columns 

Contents 

Mode 

Range 

1-7 

Distance  along  the  flight  track 
to  terminate  computation 

Integer 

>0 

8-14 

Offset  in  the  X-direction 

tl 

tt 

15-21 

Offset  in  the  Y-direction 

tt 

•t 

D-2 


TABLE  X 


CONTOUR  CONTROL  CARD  FORMAT 


Columns  Contents 


1-7 

Distance  along  flight  track  for 
the  first  data  point 

8-14 

Increment  along  flight  track 
between  successive  data  points 

15-18 

NEF  value  of  the  1st 

NEP 

contour 

19-22 

"  "  "  "  2nd 

tt 

V 

23-26 

ti  ii  n  it 

tt 

Tt 

27-30 

»  "  "  "  4th 

m 

IT 

31-34 

n  "  "  "  5th 

tt 

t! 

35-38 

1.  ..  ..  .1  6th 

tt 

ft 

39-42 

"  ”  "  ”  7th 

tt 

ft 

43-46 

"  "  »  "  8th 

tt 

ft 

47-50 

II  II  II  II  gth 

it 

ft 

51-54 

"  "  "  "  10  th 

tt 

It 

Mode  Range 

Integer  >0 

ii  it 


Real 

tt 

it 

tt 

ft 

tf 

tt 

ft 

tt 


ft 

tt 

tt 

tt 

ft 

ft 

ft 

tt 

ft 

tt 
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TABLE  XI 


DATA  CARD  FORMAT 


Columns 

1-7 

8-14 

Distance 

to 

Contents 

/.-distance 

1st 

NEF 

contour 

Mode 

Integer 

It 

Range 

>0 

tl 

13-21 

If 

2nd 

tl 

ft 

It 

•I 

22-28 

11 

It 

3rd 

tt 

tt 

It 

»? 

29-35 

ft 

It 

4th 

If 

11 

tl 

tt 

36-42 

It 

tt 

5th 

tt 

It 

It 

IT 

i»3_49 

IT 

tt 

6th 

tt 

11 

ft 

11 

50-56 

17 

If 

7th 

tt 

tt 

II 

It 

57-63 

It 

It 

8th 

It 

It 

ft 

tt 

64-70 

?l 

tt 

9  th 

tt 

tt 

Vt 

71-77 

ft 

It 

10th 

ft 

ft 

tt 

*t 

Columns 

1-7 

Closing 

TABLE 

CONTOUR  CLOSING 

Contents 

distance  of  1st 

XII 

CARD  FORMAT 

NEF  contour 

Mode 

Integer 

Range 

>0 

8-14 

If 

IT 

tt 

2nd 

Tt  11 

•I 

ft 

15-21 

tl 

It 

tt 

3rd 

ft  II 

:t 

?! 

22-28 

If 

If 

tl 

4th 

11  If 

it 

If 

29-35 

It 

11 

It 

5th 

tl  II 

it 

ft 

36-42 

tl 

If 

11 

6th 

tl  tt 

IT 

It 

43-49 

tl 

IV 

It 

7th 

tt  tt 

tt 

If 

50-56 

M 

II 

tt 

8th 

If  tt 

11 

ft 

57-63 

It 

tt 

ft 

9t'n 

It  II 

tt 

It 

64-70 

tt 

It 

»f 

10th 

ft  It 

It 
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b.  Merged  contour  control  card  HI.  This  card  specifies  the 
location  along  the  flight  path  where  contour  calculations 
are  to  commence  and  the  increment  along  the  flight  track 
for  which  the  contour  distances  are  to  be  computed.  The 
set  of  NEF  contour  values  for  which  the  contour  distances 
are  to  be  computed  must  also  be  specified  on  the  card. 

c.  Merged  contour  control  card  H2  .  This  card  specifies  the 
point  at  which  the  merging  process  is  to  be  terminated  as 
well  as  the  offset  of  the  runways  (flight  paths). 

For  each  of  the  two  contours  to  be  merged  the  data  to  be 
processed  must  appear  in  the  following  sequence: 

a.  A  title  card.  This  card  identifies  one  of  the  contours  *  * 
be  merged  (up  to  80  characters  are  allowed). 

b.  A  contour  control  card.  This  card  details  the  starting 
location  and  the  interval  along  the  track  of  successive 
data  cards  and  the  NEF  values  of  the  successi\'e  da^a  poirts 
on  each  card. 

c.  Data  cards  containing  the  distance  along  the  flight  trad 
and  the  distances  to  the  contours  for  the  given  trac’'. 

d.  A  blank  card. 

e.  A  card  with  the  closing  point  of  the  various  NEF  contours. 
The  columns  pertaining  to  contoui*s  which  do  not  close  In 
the  interval  under  consideration  are  left  blank. 

f.  The  same  sequence  as  a  through  e  but  now  for  the  second 
track. 

The  various  card  formats  are  described  in  Tables  VIII 
through  XII. 

All  distances,  which  are  read  into  or  generated  by  the 
program  must  be  uniquely  defined.  In  order  to  make  input  to 
MERGENEF  compatible  with  the  output  of  the  NEF  program  the 
data  describing  the  contours  for  each  runway  are  accepted  in 
the  same  form  as  the  output  generated  by  the  NEF  program.  For 
purposes  of  output,  however,  a  unique  origin  must  be  specified. 
MERGENEF  assumes  that  the  point  of  takeoff  roll  of  the  first 
runway  for  which  data  are  entered  is  the  new  origin.  As  the 
contours  obtained  will  generally  be  plotted,  the  output  of 
MERGENEF  has  been  organized  with  this  use  for  the  data  in  mind. 


To  find  the  correct  offset  distances,  do  the  following: 

1 .  Choose  one  of  the  two  runways  to  be  the  datum  line. 

2.  Standing  at  the  beginning  of  this  runway  and  looking  along  the 
flighttrack,  is  the  other  runway  on  the  left  or  on  the  right? 

3.  If  the  other  runway  is  on  the  left,  YO  is  negative. 

If  the  other  runway  is  on  the  right,  YO  is  positive. 

4.  Still  looking  down  tho  flighttrack,  is  the  beginning  of  the  other 
runway  ahead  of  you  or  behind  you? 

5.  If  the  othe-  runway  is  ahead,  XO  is  positive. 

If  the  other  runway  is  behind,  XO  is  negative. 

If  the  other  runway  is  right  next  to  you,  XO  is  z3ro. 


EXAMPLE:  Runways  25-L  and  25-R  East  Mythical  Municipal  Airport  1970. 


1. 

Choose  1st.  runway  to  be  25-L. 

1. 

Choose  1st.  runway  to  be  25-R. 

2. 

25-R  is  to  the  right. 

2. 

25-L  is  to  the  left. 

3. 

YO  =  +  7000' 

3. 

Y0  =  -  7000' 

4. 

Beginning  of  25-R  is  ahead. 

4. 

Beginning  of  25-L  is  behind. 

5. 

XO  =  +  2000' . 

5. 

X0  =  -  2000'. 

For  results  compare  Fig.  D-2. 

For  results  compare  Fig.  D-3. 

FIGURE  D-l .  ALGORITHM  TO  FIND  CORRECT  OFFSET  DISTANCES. 


The  output  for  each  point  along  the  flight  track,  whi;'h  Is 
considered,  by  the  program  as  the  datum  line,  consists  of  sets 
of  numbers.  The  first  row  gives  the  contour  distances  to  the 
right  of  the  right  hand  runway.  The  next  two  rows  give  the 
contour  distances  between  the  two  runways,  those  closest  to 
the  right  hand  runway  first,  and  the  final  row  gives  the 
contour  distances  to  the  left  of  the  lefthand  runway.  All 
references  to  left  and  right  are  for  a  person  looking  down  the 
flight  track.  All  distances,  which  are  to  the  left  of  the 
datum  are  given  a  negative  sign. 

The  offset  of  the  runways  is  also  referenced  t0‘  the  same 
origin,  and  a  displacement  in  the  X-dir-ection  (along  the 
track)  given  by  the  appropriately  signed  distance.  The  offset 
in  the  y-direction  follows  the  same  convention  as  is  used  in 
specifying  the  output.  If  the  second  runway  is  to  the  right 
of  the  first,  the  sign  of  the  y-offset  is  positive,  on  the 
other  hand,  if  the  second  runway  is  to  the  left  of  the  first 
the  y-offset  is  negative.  A  simple  algorithm  to  find  the 
correct  offset  values  is  given  in  Fig.  D-l. 

The  examples,  which  we  have  carried  through  in  Appendix  A, 
can  also  serve  here  to  illustrate  the  procedure.  We  assume 
that  each  of  the  two  examples  gives  the  operations  on  only  one 
of  two  parallel  runways.  The  desired  contours  are  then  obtained 
by  using  the  MERGENEF  program.  The  results  are  shown  in 
Figs.  D-2  and  D-4  for  the  1970  and  1990  forecasts  respectively. 
Fig.  D-3  shows  an  alternate  computation  of  the  1970  contour 
to  demonstrate  how  either  one  of  the  runways  can  be  chosen  for 
the  datum. 

A  few  words  of  caution  are  in  order.  The  MERGENEF  program 
is  in  essence  an  interpolation  scheme  which  generates  a  summation 
curve  from  two  overlapping  curves.  It  is  clear,  that  such  a 
program  is  most  accurate  when  a  large  number  of  points  are 
known  between  which  the  desired  curve  must  be  fitted.  The 
routine  becomes  therefore  inaccurate  when  only  a  few  contours 
are  still  existence  at  a  particular  X-value.  In  particular  if 
only  one  or  two  contours  are  left  (at  fairly  large  X-values) 
the  program  cannot  compute  the  contours  accurately  and  the 
resultant  output  shows  a  diverging  rather  than  a  closing 
contour.  The  reason  is  that  at  those  points  the  program 
performs  on  extrapolation,  rather  than  an  interpolation. 

Similarly  the  contour  with  the  smallest  NEF  value  is  always 
extrapolated  (it  is  always  further  away  than  the  contours  it 
was  generated  from)  and  therefore  it  is  better  to  supply  the 
program  with  one  lower  value  NEF  contour  than  is  of  interest 
(i.e.  if  the  25  NEF  contour  is  desired  MERGENEF  will  compute 
its  location  more  accurately  If  the  20  NEF  contour  Is  also 
known  to  the  program) . 
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As  it  is  necessary  to  punch  the  NEF  contour  information  into 
cards  for  input  to  MERGENEF  it  is  more  convenient  land  it 
eliminates  errors)  to  have  the  computer  punch  out  the  necessary 
information  while  the  F^F  contours  are  computed. 

One  way  to  achieve  this  is  shown  in  Fig.  D-5.  The  modx Med 
NEF  routine  MODNEF  will  punch  out  the  cards  needed  by  MERGENEF 
in  the  right  sequence.  It  is  up  to  the  user,  however,  to 
supply  the  identifying  and  control  cards  needed  to  run  the 
merging  program.  This  can  be  added  after  all  cards  are  punched 
out  and  sorted  into  parallel  runway  sets.  If  cards  are  not 
the  medium  chosen  to  punch  the  MODNEF  output  (but  for  example 
a  tape  or  disk),  the  user  will  have  to  make  sure  that  the 
needed  control  information  is  inserted  in  the  MODNER  output 
stream.  Such  extensive  modifications  are  not  within  the  scope 
of  this  text. 

A  flowchart  of  the  MERGENEFF  program  is  given  in  Fig.  D-6, 
a  program  listing  in  Fig.  D-7. 


D-8 


DOUBLE  PUNWAYSf T  25l-25»  197C 

START*  C  INTERVAL-  2000  CONTOURS  25.  30.  35.  40.  45.  50.  55. 

STOP  -  16000  X,Y  OFFSET  f!F  2ND  RUNWAY-  2000  ,  7000 


EAST  MYTHICAL  MUNICIPAL  A  |  RPfjfl  T  -  FORECAST  1970  25-L 


STAPT- 

0  INTERVAL* 

2000 

CONTOURS 

25. 

30. 

35.  40.  45.  50. 

55 

EAST  MYTHICAL 

MUMCIPAL  AIRPORT  - 

FORECAST 

1970  25- 

R 

STAF  T- 

0  INTERVAL- 

2000 

CONTOJRS 

25. 

30. 

35.  40.  45.  50. 

55 

DISTANCE  Tfl  MERGED 

CONTOUR 

(IN  FFCT) 

FROM 

FIRST  RUNWAY 

CONTOUR  VALUES 

25. 

30. 

35. 

40. 

45. 

50. 

55. 

X 

VALUE-  2)00. 

TO  RIGHT 

17117 

10P45 

0700 

P907 

P  173 

764  6 

7321 

BETWEEN 

4057 

5091 

5A29 

6346 

4771 

RUNWAYS 

2575 

1594 

92? 

510 

07 

TO  IFFT 

—4  54 ft 

-3401 

-2423 

-1492 

-920 

-407 

-236 

X 

VAL  UF  =  4OC0. 

TO  FIGHT 

l  1  4  40 

103C4 

0421 

8S91 

7020 

7497 

72  16 

Bf  T  WFEN 

4531 

44  14 

6372 

64*9 

6912 

RUNWAYS 

211  3 

1307 

736 

390 

36 

TO  LEFT 

-AC  80 

-2*-H7 

-2065 

-  1 3  3  8 

-732 

-384 

-178 

X 

VALUE-  4C00 • 

TO  FIGHT 

1  IC71 

90^3 

9063 

»3P7 

7732 

73P4 

7178 

8L  TWEEN 

3  725 

49J  A 

5696 

6764 

660° 

4  96  3 

RUNWAYS 

2°12 

1  769 

1057 

604 

289 

13 

TO  LEFT 

-3631 

-2617 

-1744 

-1054 

-590 

-302 

-137 

X 

VALUE*  9000. 

TO  RIGHT 

MT  24 

96  I  4 

R745 

8053 

7590 

7  3C2 

7137 

3ETWFEN 

4277 

4742 

5944 

6395 

6710 

6984 

RUNWAYS 

74  17 

1492 

1012 

686 

366 

41 

TO  (EFT 

-3  Jf  | 

-24 r  | 

-1577 

-004 

-t  P9 

-3(5 

-164 

X 

VALUE-  more. 

TO  RIGHT 

1 0  ?  6  5 

930S 

P475 

7094 

76*9 

7365 

7 1  04 

HI TWEFN 

4  4  30 

5415 

5  »BA 

6313 

4633 

6958 

RUNWAYS 

2-301 

14  71 

1213 

754 

447 

171 

Ttl  LEFT 

-3 1 1  3 

-2  BO0 

-107) 

-1192 

-746 

-463 

-2C1 

X 

VALUE-  12000. 

TO  RIGHT 

10(63 

ocp. 

P  470 

8195 

7746 

746? 

7201 

BETWEEN 

4076 

5779 

6244 

4551 

RUNwA VS 

2261 

1447 

802 

34? 

6800 

T(  IFFT 

-4  >22 

-  3  1  4  ■> 

-226) 

-1466 

-509 

-  342 

c 

X 

VAL'JC*  14096. 

TO  PK.HT 

1  1  7*3 

1016  7 

5764 

3467 

T030 

7342 

c 

HE  TWEEN 

4640 

44  10 

6176 

PUNWAVS 

2(.2-3 

1408 

6  31 

664* 

TU  LCFT 

-507*, 

-  i«67 

-  ?  40  3 

-139A 

-4*7 

0 

0 

X 

VALUE*  160CO. 

TO  FIGHT 

1211° 

1 0  A®  4 

941  9 

9401 

7687 

0 

0 

BETWEEN 

4171 

5576 

6360 

RUT  WAYS 

767  1 

1277 

3  6  A 

TU  LEFT 

-•4  =>  l  1 

-4040 

-242  A 

-1288 

-472 

0 

*\ 

FIGURE  D-2.  EXAMPLE  OF  COMPUTER  OUTPUT  - 
DISTANCE  LISTING,  EXAMPLE  1 
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MERGED  CONTOUR 


70UBt  6  PlINWAVSE  T  25L-24R  1070  ALTFPNATF  COMPUTATION 

START*  P  !NUMV/>1  =  2300  CONTOlMS  25.  JO.  J5.  40.  45.  5C.  55. 

STOP  *  1600C  X.V  ()F  P  Sf  T  OF  2ND  MUNWAV*  -2000.  -TOCC 

iAST  mvthical  municipai  mnrur  -  fipfcast  1070  25-p 

STAPT*  0  INTERVAL*  2000  CC.'NTPJStS  25.  JO.  J5.  40.  45.  50.  55. 

:  AST  MV  TM  f  T  AL  ,''JM('l0Al  At»PPPT  -  TIPfCAST  1070  25-L 

START*  r  tNTf  KVAl  =  2000  c."'i  T3JRS  25.  JO.  35.  40.  45.  5C.  55. 


I  ISTANO  TO  *rpf.lL>  L  ■  INTI  l»W  UN  FEfT»  HfJV  F  I K  i  T  PUNRAV 

3CMPU8  V  Al  I Jf  S  ?4.  10.  16.  40.  45.  50.  55. 

X  VALUF*  0. 

TP  F  tl.MT  5  117  Jf»46 

A?  TV.FPN 
MUNfcA VS 

Tf  I  r  P  T  -11  54t-  -  1040  1 

X  VALUE*  2A0 0. 

TC  P  I  l.HT  4  *  4  0  Mc4 

IF  TV.FPN 
MONV.AVS 

Ti  IfPT  -llOPO  -0«pv 

X  value*  4°0 n . 

TP  P I  f'.t  <  T  4071  ;cf  I  70f.l  l  J  07  7  37  304  1  7rt 

Of.TV.FFN  -4.0P7  -4771  -5047  -4.305  -5  7 1 0  -50H4 

Rl  INV.  A  V  S  -1274  -21-1  -1J03  -7  35  -  Jn0  -  Jf. 

Tf  LIFT  -  |  P  6  11  -04l?  -f  74*.  —  *1 C  54  -7431  -73  '“2  -TIJ7 


2  7‘>a 

1  <!C  7 

1  1  73 

644 

321 

-4424 

-4405 

-4"70 

-44«3 

-6012 

-7047 

-  1  OOA 

-1173 

-443 

-226 

-•>42  1 

-A412 

-  7023 

-7407 

-7234 

1  5°1 

020 

407 

736 

-A  *^6 

-  S6'l? 

-<■  ’63 

-6600 

-5  >6  3 

• .?  A  f*  •» 

—  1  6  ®  4 

-027 

-510 

-H  7 

-  F  30° 

-7732 

-  Tl»«. 

-7  1  7« 

x  v a t  ur  *  4  K<’. 

TP  t  ■  i>*<T  ‘*7 

BP  TwFPN 
P.INV.A  VS 

TC  LFPT  -I" 


26  |  4 

1  7A» 

1  04* 

-44-7 

-  •»  a  r  t 

-  4  f ;  H  7 

-2  777 

-1  yt>' 

- 1  C  54 

-0401 

-41  ft 

-  7404 

6  C  6 

33? 

1  37 

6  11  J 

-66  3  3 

-605A 

-•04 

-?r  > 

-1  3 

7?>A0 

-73?  5 

-  71?  4 

X  VA1UE*  AP?P. 

TP  p  |  i',nT  JJf.6  7144  l  ‘74  spo  3»>5  144 

BE  Tv.Fr*.  -4C'<“  -3075  -‘■P'"'  -*  245  -5452  -4*7- 

PUNPAVS  —24?  5  -1*,A.  -1011  -655  -344  -41 


Ti  lU  T 

- 1 0<  1? 

-  n  p  < 

»  V Al UF  =  IOOO-. 

TC  R  Ji.mT 

i»?  3 

/ur  i 

(it  Twe  r *> 
PI.N4A  vr. 

Tf  IFFT 

-11  ft  7 

-1^1^ 

■  V Al UF  *  17000. 

TI.  P  If.ll T 

<,7M 

3  1  6  > 

BE  Tv.FF’1 
PlINWA  VS 

T  (  IFFT 

-i/C  7/. 

-  1  0  -4  7 

X  VAL'JF*  140,'". 

TI  PIC.Hf 

4 1 1 0 

3  A  P  5 

111  T  nF  F \ 
JllNPA  VS 

T  1  rp  r 

-17-11 

-1 !  14‘ 

X  VAI  HE  *  16"00. 

Ti  Pir.Pl 

7  •>  1  4 

4704 

1[  T * F F ? 

RU NV.A  VS 

T  IfPT 

-1  4  »>  7  r 

-  1 1  173 

“'*7  1 

-F  l°? 

-  7746 

-7443 

-  7  2  r  1 

|  47- 

1  1  05 

7<»<. 

46  3 

201 

4  7  14 

-  5  4  4 ,7 

-  1  37 

-4447 

70?  3 

-172C 

-755 

-  44  f 

-171 

1 

-  “4t.t 

-va.’O 

-7  3(.  7 

n 

??i.4 

1447 

34? 

r\ 

4  i  70 

-  54 4  { 

-4  un 

2  14  1 

-|4..o 

-“03 

-  347 

fin 

.  0  Vi  J 

.’6  A  7 

6 

*5 

7  40'. 

1401 

417 

0 

n 

4326 

7  J?i 

047H 

-4  722 

-  1  4  7  I 

-  J  7  A  F 

-66  3  l 
-4  3  0 
-  74  72 

0 

0 

24-  ’ 

1301 

4  76 

0 

0 

42  1  » 

2  A  1  4 
<•34  1 

-6 -'06 
-1  304 

-“147 

-  174 

0 

c 

A 

FIGURE  0-3, 


EXAMPLE  OF  COMPUTER  OUTPUT  -  MERGED  CONTOUR 
DISTANCE  LISTING,  ALTERNATE  COMPUTATION  OF 
EXAMPLE  I 
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DOUBLE  RUNWAYSET  25L-25R  1990 

START*  0  INTERVAL*  2000  CONTOURS  25.  30.  35.  40.  45.  50.  55. 

STOP  *  100000  X,Y  OFFSET  OF  2ND  RUNWAY*  2000,  7000 

EAST  MYTHICAL  MUNICIPAL  AIRPORT  -  FORECAST  1990  25-1 

START*  0  INTEPVAL*  2000  CONTOURS  25.  30.  35.  40.  45.  50.  55. 

EAST  MYTHICAL  MUNICIPAL  AIRPORT  -  FQRFC AST  1990  25-R 

START*  0  INTEPVAL*  2^00  CONTOURS  25.  33.  35.  40.  45.  50.  55. 


DISTANCE 

TO  MERGED  CONTOUR  I  IN 

FEFTI 

FROM  FIRST  RUNWAY 

CONTOUR  VALUES 

25. 

30. 

35. 

40. 

45. 

60. 

65. 

X 

VALUE*  20C0. 

TC  RIGHT 

1  <>  7  7  0 

1  3  l  7« 

1177*, 

10699 

9810 

9021 

7333 

RFTWFEN 

4117 

4970 

5690 

RUNWAYS 

2352 

l  550 

910 

TO  LEFT 

-6‘»15 

-5253 

-4062 

-7117 

-2279 

-1550 

-051 

X 

VALllF*  4000. 

TO  RIGHT 

13665 

I  2228 

11051 

10113 

9277 

8550 

7951 

BETWEEN 

3672 

4716 

5458 

5091 

RUNWAYS 

2771 

1804 

1  141 

612 

TO  LEFT 

-5714 

-4457 

-7431 

-2570 

-17P8 

-1144 

-654 

X 

VALUE*  6000. 

TO  RIGHT 

I2LF2 

I  1442 

10425 

956H 

879  7 

5144 

7654 

BETWEEN 

4352 

5239 

5  560 

6387 

RUNWAYS 

24  17 

1427 

1076 

7  38 

TO  LEFT 

-5325 

-4167 

-32  0° 

-2763 

-1619 

-1064 

-803 

X 

VALUE*  ROCO. 

TO  PIGMT 

1?  loo 

It  156 

10204 

9362 

8619 

d  065 

7803 

nr  TWEEN 

4567 

5375 

5924 

6261 

RUNWAVS 

346« 

2324 

1446 

793 

Til  LEFT 

-5?77 

-4600 

-3052 

-3367 

-2722 

-1442 

-580 

X 

VALUE*  1COCO. 

TO  FIGHT 

12269 

11507 

10641 

10367 

9322 

8442 

7H«P 

BTTWf.EN 

4504 

5529 

6201 

RUNWAYS 

2346 

1362 

703 

TO  LEFT 

-7  369 

-6622 

-5507 

-3675 

-2821 

-1363 

-711 

X 

VALUF*  120C0. 

TC  RIGHT 

14*73 

1  3o24 

124CI 

10697 

9»75 

3370 

7712 

dl  TWFEN 

4432 

5611 

6  291 

RUNWAYS 

74J8 

1179 

362 

TG  LFFT 

-9<  4 7 

-SHOO 

-6  l  5  l 

-  3729 

-2720 

-1196 

-394 

X 

VALUE*  1 4000 . 

TO  RIGHT 

IFF  61 

15017 

13221 

10598 

9266 

9197 

7396 

BETWEFH 

4539 

5P14 

RUNWAYS 

2335 

853 

6634 

TO  LEFT 

-l 1 c72 

- 10RF  7 

-6352 

-  3fc°4 

-2000 

-897 

-91 

X 

VALUE*  16000. 

TO  RIGHT 

19265 

17106 

1  .77  31 

10867 

9117 

T9PI 

7091 

BETWEEN 

4707 

6  144 

RUNWAYS 

2343 

406 

6«20 

TO  LEFI 

-I4iri 

-10472 

-63B7 

-3751 

-1893 

-  4  84 

C 

X 

VALUE*  iRoro. 

TO  RIGHT 

22685 

185R2 

17HA  J 

10764 

8895 

7484 

0 

BETWEEN 

4977 

RUNWAYS 

I6P8 

6  6C0 

TO  LEFT 

- 1 7  I R  l 

-11626 

-6709 

-3600 

-1594 

0 

0 

FIGURE  0-4.  EXAMPLE.  OF  COMPUTER  OUTPUT  -  MERGED  CONTOUR 
DISTANCE  LISTING,  EXAMPLE  2 


D-n 


DISTANCE  TO  MERGEO  CONTOUR  IIN  FEET)  FROM  FIRST  RUNWAY 


CONTOUR  VALUES 

25. 

30. 

35. 

40. 

45. 

50. 

55. 

X  VALUE*  20000. 

TO  RIGHT 

262  IP 

19570 

14000 

10617 

8506 

0 

0 

BETWEEN 

5334 

RUNWAYS 

1251 

TO  LEFT 

-2C2P3 

-12450 

-6926 

-3426 

-1205 

0 

0 

X  VALUE*  22000. 

TC  RIGHT 

28429 

19955 

139  30 

10427 

8203 

0 

0 

BETWEEN 

5765 

RUNWAYS 

738 

TO  LEFT 

-21410 

-12929 

-6828 

-3204 

-697 

0 

0 

X  VALUE*  24000. 

TO  RIGHT 

28373 

19924 

13835 

10203 

7694 

0 

0 

BETWEEN 

6271 

RUNWAYS 

0 

TO  LEFT 

-21393 

-13013 

-6715 

-2928 

-49 

0 

0 

X  VALUE*  26000. 

TO  RIGHT 

28848 

20175 

13743 

9925 

7037 

0 

0 

BETWEEN 

RUNWAYS 

6993 

TO  LEFT 

-22117 

-13051 

-6587 

-2598 

0 

0 

0 

X  VALUE*  28000. 

TO  RIGHT 

29790 

20271 

13604 

95  85 

0 

0 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-22861 

-13193 

-6444 

-2184 

0 

0 

0 

X  VALUE*  30000. 

TO  RIGHT 

30038 

2023C 

13428 

9168 

0 

0 

0 

BETWEEN 

3868 

RUNWAYS 

2490 

TO  LEFT 

-23015 

-13161 

-6259 

-1663 

0 

0 

0 

X  VALUE*  32000. 

TO  RIGHT 

30072 

20153 

1  3232 

86  44 

0 

0 

0 

BETWEEN 

4797 

RUNWAYS 

1325 

TO  LEFT 

-23C73 

-13072 

-6040 

-1000 

0 

0 

0 

X  VALUE*  34000. 

TO  RIGHT 

30183 

20O7P 

13012 

7979 

0 

0 

0 

BETWEEN 

5794 

RUNWAYS 

78 

TO  LEFT 

-23183 

-12992 

-5794 

-406 

0 

0 

0 

X  VALUE*  360C0. 

TO  RIGHT 

30246 

19Q79 

12759 

7405 

0 

0 

0 

BETWEEN 

RUNWAYS 

6°35 

TO  LEFT 

-2321° 

-12891 

-5517 

0 

0 

0 

0 

FIGURE  D  -4. 


EXAMPLE  OF  COMPUTER  OUTPUT  -  MERG  E  D  CO  NTOU  R 
DISTANCE  LISTING,  EXAMPLE  2  (CONTINUED) 


DISTANCE  TO  MERGED  CONTOUR  UN  FEET!  FROM  FIRST  RUNWAY 


CONTOUR  VALUES 

25. 

30. 

35. 

40. 

45. 

50. 

55. 

X  VAt UE=  380C0 . 

VC  n I GHT 

BE • WEEN 
RUNWAYS 

30242 

19862 

12476 

0 

0 

0 

0 

TO  LEFT 

-23219 

-12766 

-5203 

0 

0 

0 

0 

X  VALUE*  40000. 

TD  RIGHT 

BETWEEN 

RUNWAYS 

30251 

19736 

12154 

0 

0 

0 

0 

TO  LEFT 

-23226 

-12633 

-4B39 

0 

0 

0 

0 

X  VALUE*  42000. 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30251 

1 c  5F  7 

117P6 

0 

0 

0 

0 

TO  LEFT 

-23221 

-124P8 

-4424 

0 

0 

0 

0 

X  VALUE*  44000. 

TD  RIGHT 

BETWEEN 

RUNWAYS 

30242 

19447 

11359 

0 

0 

0 

0 

TO  LEFT 

-22211 

- 1 2  330 

-3956 

0 

0 

0 

0 

X  VALUE*  46C0C. 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30235 

19285 

lC85o 

c 

0 

0 

0 

TO  LEFT 

-23203 

-12161 

-3397 

0 

0 

0 

0 

X  VALUE*  48000. 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30232 

19112 

10273 

0 

0 

0 

0 

TO  LEFT 

-232C0 

-11979 

-2637 

0 

0 

0 

0 

X  VALUE*  50000. 

TO  RIGHT 

30234 

13927 

9533 

0 

0 

0 

0 

BETWEEN 

4130 

RUNWAYS 

242 

TO  LEFT 

-23203 

-11785 

-1199 

0 

0 

0 

0 

X  VALUE*  52000. 

TO  RIGHT 

30244 

18730 

7927 

0 

0 

0 

0 

BETWEEN 

6805 

RUNWAYS 

163 

TO  LEFT 

-23213 

-11578 

-843 

0 

0 

0 

0 

X  VALUE*  54000. 

TO  RIGHT 

30263 

18519 

7680 

0 

0 

0 

0 

BETWFEN 

68  3  > 

RUNWAYS 

123 

TO  LEFT 

-  '2;' j2 

-11355 

-585 

n 

0 

0 

0 

FIGURE  D-4.  EXAMPLE  OF  COMPUTER  OUTPUT  -  MERG  ED  CONTOUR 
DISTANCE  LISTING,  EXAMPLE  2  (CONTINUED) 


D  - 1  3 


DISTANCE  TO  MERGED  CONTOUR  (IN 

FEET! 

from 

FIRST  RUNWAY 

CONTOUR  VALUES 

X  VAL UE  =  56000. 

25. 

30. 

35. 

40. 

45. 

5C . 

55. 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30292 

182^2 

7445 

6879 

7ft 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  5  P000 . 

-23263 

-11  115 

-32ft 

0 

0 

0 

0 

TO  RIGHT 

BF  TWEEN 
RUNWAYS 

30336 

1  8049 

72  30 
6925 

24 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  60000. 

-23309 

-10859 

-88 

0 

0 

0 

0 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30397 

17789 

7050 

6979 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  62000. 

-23374 

-  0  5  F  3 

0 

0 

0 

0 

0 

TO  RIGHT 

BF  T  WE  EN 
RUNWAYS 

304B0 

1750ft 

0 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  64000. 

-23462 

-10285 

0 

0 

0 

0 

0 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30593 

17205 

0 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  66000. 

-235FO 

-9962 

3 

0 

0 

0 

0 

TP  RIGHT 

BETWEEN 

RUNWAYS 

30743 

16876 

0 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  68000. 

-23740 

-9610 

0 

0 

0 

0 

0 

TO  RIGHT 

BETWEEN 

RUNWAYS 

30944 

16517 

0 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  70000. 

-23951 

-9223 

0 

0 

0 

0 

0 

TO  RIGHT 

BETWEEN 

RUNWAYS 

31210 

16121 

0 

0 

0 

0 

0 

TO  LEFT 

X  VALUE*  72000. 

-24234 

-B7e3 

0 

0 

0 

0 

0 

TO  RIGHT 

BETWEEN 

RUNWAYS 

31  574 

15682 

0 

0 

0 

0 

0 

TO  LEFT 

-29622 

-8  306 

0 

0 

0 

0 

0 

FIGURE  D-4. 


EXAMPLE  OF  COMPUTER  OUTPUT 
DISTANCE  LISTING,  EXAMPLE  2 


-  MERGED  CONTOUR 
(CONTINUED) 


D-14 


DISTANCE 

TU  MERGED 

CONTOUR 

(IN 

FEET! 

FRO*  F IRST 

RUNWAY 

CONTOUR  VALUFS 

2  5. 

3f  • 

35.  4 

0. 

45. 

50. 

55 

i . 

X  VALUE®  74000. 
TO  RIGHT 

32061 

15183 

0 

0 

0 

0 

o 

BETWEEN 

RUNWAYS 

TO  LEFT 

-26130 

-7760 

J 

0 

0 

0 

0 

Y.  VALUE*  76000. 
TO  RIGHT 

32  757 

14626 

0 

c 

0 

0 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-26899 

-7  008 

0 

0 

0 

C 

0 

X  VALUE®  78000. 
TO  RIGHT 

33F  78 

1  3°30 

o 

0 

0 

0 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-27133 

-6242 

0 

0 

0 

0 

0 

X  VALUE®  aoooo. 
TO  RIGHT 

36871 

13024 

0 

0 

0 

0 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-29412 

-4996 

0 

0 

0 

c 

0 

X  VALUE®  82000. 
TO  RIGHT 

45  890 

1  1460 

0 

0 

0 

c 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-46585 

1373 

672 

-43P7 

0 

0 

0 

0 

0 

X  VALUE®  84000  • 
TO  RIGHT 

105704 

10465 

0 

0 

0 

0 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-99396 

203 

-3245 

0 

0 

0 

0 

0 

X  VALUE*  86000. 
TO  RIGHT 

108C15 

10044 

0 

0 

c 

0 

0 

BETWEEN 

RUNWAYS 

TO  LEFT 

-102114 

6700 

181 

-2*02 

0 

0 

0 

0 

0 

X  VALUE®  88000. 
TO  RIGHT 

l 1Q4C° 

ofol? 

0 

0 

0 

0 

0 

BtTWEFN 

RUNWAYS 

TO  LEFT 

-104830 

6822 

156 

-2346 

0 

0 

0 

0 

0 

X  VALUE®  9C000. 
TO  RIGHT 

1 1 2  8  Q9 

9172 

0 

0 

0 

0 

0 

BETWEEN 

RUNWAYS 

TG  LEFT 

-107453 

6847 

129 
-1  6  T6 

0 

0 

0 

0 

0 

FIGURE  D-4.  EXAMPLE  OF  COMPUTER  OUTPUT  -  MERGED  CONTOUR 
DISTANCE  LISTING,  EXAMPLE  2  (CONTINUED) 


D  - 1  5 


v  I STANC E 

TO  y{P 

~»F f  CkiNTnu»  TIN 

FEET) 

FRQN  Ei*ST 

runway 

CON TO UK  VALUES 

25. 

30. 

35.  ’»0. 

45. 

50.  55. 

X  VAl  IE=  ATI 

« 

TO  RIGHT 

115500 

P  T?r 

3  0 

C 

0 

0 

9f TWEEN 

6*74 

RUNWAYS 

co 

TT  LEFT 

-11015B 

-  1352 

T  C 

0 

n 

0 

X  VALUF=  C4Q00 

• 

TO  FIGHT 

11P243 

a  ?5<: 

3  J 

0 

C 

r* 

BETWEEN 

6:>C4 

RUNWAYS 

t  5 

TO  LETT 

-  1 1  3 1  1  5 

-  ac4 

3  0 

n 

0 

0 

X  VAL 'JE=  ‘VftfTO 

• 

Tf  RIGHT 

;  ?  •  i 

•»-»p  « 

0  0 

0 

0 

n 

BETWEEN 

RUNWAYS 

2~> 

TO  l EFT 

-1 1  f  236 

-3=2 

?  0 

0 

0 

0 

X  VAL UE  =  9-iOOO 

• 

Tr  EIGHT 

1?4?F2 

"Mf  ■* 

> 

G 

0 

0 

BE  TWEE’. 

RUNWAYS 

6474 

TO  l EF T 

-I1R5P? 

/*■ 

°  0 

0 

r 

0 

FIGURE  0-4.  EXAMPLE  OF  COMPUTER  OUTPUT  -  MERGED  CONTOUR 
DISTANCE  LISTING,  EXAMPLE  2  (CONTINUED) 


D-l  6 


r>  nr«or*nr^ 


c  - 

c 

c 

C  THIS  IS  A  SiPHF  MODIFif ATI')N  OF  T«E  NEF  lOVTSr*  PPUTINE 

C  I K  THIS  Fr»*  Th£  ROUTINE  PUNCHES  OUT  THE  CA*DS 

C  PECJIM}  *0*  W  »£»S£\E=  C7NTDJ?  pEP&EVG  *DJTl*E 

r  Aii  ow*  functions  of  thf  hcjtine  are  u\:haigeo 
C 

c  - 

( 

r  SuhPTUTIVES  AX')  FuKCllUN  $J:»?»rGRAPS  caileo 

C  CAPLIN,  fPO,  ANE-,  XNEW.  IASTPT 

t 


COH  1 
t  ■*“  2 

CCM  3 
4 

COM  S 
c^h  6 

r.I“CNSIDN  Vf=r3NC7»,  I  v  l  7 ) ,  xClOSFfM 

r  - 


XiTA=T,  XF  VI),  DX,  VAC,  CtKlT  U*1 ,  N»*AXf4| 

Cf-VIN  ».TET|P),  NFFTI4I,  \4lH121.  N?0rftl»l,  ACAICI12I,  PCAICII2I 
Cl  “H  )?i  f)SC\A,M?''l ,  ALTN4XI121,  »PxXAN|l>),  PNIVAVI25I,  APTPI50) 
t'-M)N  PPlcfEri,  CPTPI52I,  LDEMS3»,  2II2I,  P»»*-SPL  I  121  .  IDFNTII3) 
C:*«*(“fiN  OSCC^TI  3,?0»  ,  AITDI-I3.2DI,  PfiUDI R  t  3,20  I  ,  £PV:>mP,?Dl 
r.-VMn*i  PMFFF  135,2,25)  ,  UUn.M?),  P0»ERI 10,2, 12» 


THF  MUY  irClS  IS  INTP^OJCID  FO*  USE  mITH  THE  PUNCHING  CF  THE  LAST  CAPO. 
DIMENSION  IXLISC7! 


PEAL  x.FFl 

,  NFf? 

N^FC.INI  1  1 

=  25 

me  rear,  i?) 

=  3- 

f.f  rcnvi  31 

=  73; 

vff  r.nvi  4» 

s  U  o 

Nr  FOTNI5I 

=  45 

NEFCOVf 6 1 

=  5r 

NEFCUM  7  ) 

=  55 

NVAX (  l  I  = 

?5 

N  H  A  X  (  2 1  = 

51 

NVAXI?)  = 

51 

NMAX(4l  = 

30 

001  1=1, 

4 

NTCT (  I  1  = 

0 

WEFT  III  = 

n 

call  ca?oin 

THE  IDENTIFYING  I NFOPPAT I DN  l)=  THE  OATA  IS  PUNCHED  HE-lt. 

WPITE<2,2010I  IDtNT 
2010  F  CRM  AT ( 13AA) 

C 

FIGURE  D  -5  A .  PROGRAM  LISTING  -  MODNEF  CONTROL  ROUTINE 


D  - 1  7 


“OONEF 


r. 

C  JST  AMJ  JCX  A KE  !NraC<Xir  =  D  TO  CREATE  THE  3ESI8ED  PilNCHf !)  OUTPUT, 

f 

JST  =  XSTAPT 
JO*  =  OX 

ITE  <2 , 2rOl 1  JST,  JOX.  JNFFC  INfl 1,1*1,71 
20"1  FrPMATl?I7,i?i41 
C 

C  - 

X  =  XSTAHl 
NNEXT  =  7 
LIKCST  =  0 
Of'  3  N=l,7 

XCICStINl  =  O. 

3  I Y 1 N I  =  4r*0 
C 

11  CALL  APDtXl 
MCCN  *  NNEXT 
f 

DC  43  N=1,NCCN 
CCNTU*  =  nef:cmni 

Y?  =  !  Y  C  M I 
Y1  =  V?  *  0.7*7 
CALL  ANEF  IY2,  NFF?I 
CALL  ANFF  1Y1,  NEFll 

C 

21  YC  =  XKFWIY1,  F-EF1,  Y2,  NEf2»  CCNTU31 
IF  <YC  -  3000C.I  23,  23,  22 

22  YO  =  33000. 

23  CONTINUE 

IF  IY.'  -  12. tl  2*.  25,  25 
C 

24  CONTINUE 

CALL  LASTPTiX,  XCLOSEINll 
YC  =  0. 

NNEXT  *  NNFXT  -  1 
GO  TO  28 
C 

25  IF  (AdSIYC  -  Y 1 1  -  2.1  2s* ,  28,  26 

26  IF  IA9SC1Y0  -  Yil  /  YOI  -  3.7311  23,  28,  ?7 

27  Y2  =  Y1 
NFF?  =  NEFl 
Y 1  *  YP 

CALL  ANFF  ( Y l ,  NEFil 
GO  TO  21 
C 

?«  IY(N)  =  YC  +  C.5 
40  CONTINUE 
C 
C 

l INCNT  =  LINCNT  ♦  1 
IF  ( L  INCNT I  43,  42,  42 
42  l  INCNT  =  -50 

WPITE  (3,  3034)  I  DENT 
WPITE  (3,  3001)  NEFCON 


1 

I 

i 

1 

i 

f 

i 


FIGURE  D  —  5  E, .  PROGRAM  LISTING  -  MODNEF  CONTROL  ROUTINE 


D-  18 


HrONFE 


r 

43  IX  =  X  ♦  0.5 

WF  I T  t  (3,  "*005  I  IX,  (IVIN),  N=l,NC0NI 

c  - - - 

r 

f  TK)  DISTANCE  Z' I  C^*';  THF  THAC<  AND  The  C3NTDJR  DISTANCES  are  PUNCHED  mere. 

r 

W“|Tr(2,2ro2)  IX,  llYtn,  1*1,  NCON I 
200?  rr<RMAIIlll7) 
f 


f 

r. 


IF  IN'FXT)  55,  55,  53 
C 

51  X  =  X  ♦  DX 

( 

II  ( X  -  *FND)  II,  11,  55 

l. 

5  5  1  F  l  C.  =  2 

«pirr  n,  wm 

On  70  *1*1,7 

i  x  *  xi  lose  ini  ♦  " . 5 
ir  iix)  7 r,  7r,  60 

50  f.C  Tfl  161,  52,  63,  64,  65,  56,  67),  N 
61  *  h  I T  r  (3,  3011)  IX 
r.r  TO  5  0 

*  2  W« ITf  (3,  301? )  I  X 
01  T'l  60 

6  3  wr-'ITF  (1,  3 0 v  3  I  IX 

f.o  rn  6  0 

64  wr<  !Tf  (3,  1C14  )  I  x 

cr  to  60 

6  5  V.PITE  (3,  30  35)  IX 
r.fl  TO  6C 

66  WRITE  (3,  3015  »  I  X 
on  TO  6C 

(jt  W'-ITE  (3,  1C?  7 )  IX 

6  0  1  F  i  r,  =  l 

n  conti  nut. 

r  _ 


c 

(  THE  rttANK  CAP!'  AT  THE  FND  Of-  THE  CONTOUR  DISTANCF  DATASET  IS  PUNCHED  HERE 

C 

KR IT E( 2, 2C 00) 

2000  fcrmati*  *) 

C 

C  THE  CARD  CONTAINING  THE  CLDSlNG  POINTS  DE  THE  CONTOURS  IS  PUNCHED  HE  RE 

C 

Cfl  2004  1  =1  ,7 

2004  ixcisin  *  xrifSEm  ♦  *>.5 

wEITF<2,2C021  (ixnsil),  I  *  1,7) 

r 

r.  - - - - - 

on  TO  r*l,  V) ,  I T LG 
71  WRITE  (3,  303C  ) 


FIGURE  D-5C.  PROGRAM  LISTING  -  MODNEF  CONTROL  ROUTINE 


D-19 


7?  r.P  TO  2 


3COI 

FORMAT 

l‘C*. 

3?x, 

•disiam: 

l 

°x ,  * 

FLIGHT 

TRACK*,  713 

c 

3DOA 

I 

3CPS 

fofmai 

i  *  i  • . 

1SX, 

I3AA  / 

FORMAT 

( icx. 

lb. 

4X,  7IBI 

3006 

format 

i  *0*  I 

3030 

FORMAT 

i  •♦• , 

7  X  , 

•CONTOUR 

1 

/  «X, 

•  TP  AC 

K*  ) 

3C3  l 

Fl PMAT 

C  *♦*, 

?1X, 

IB) 

303? 

FORMAT 

1  •♦*, 

?cx. 

1°) 

3C?7 

FfOMAT 

c  . 

T7X, 

m 

3C?<. 

ri RMAT 

ft^x. 

i«) 

?°3  5 

format 

<•♦*. 

hW, 

ip) 

303  ^ 

FtKMAl 

(•♦*. 

<  IX, 

IR) 

30  3  7 

format 

1 '♦*, 

FOX, 

m 

le  to  :»intour  in  feet*  /  n<,  *point  on*/ 

/  'it,  * -  -  -  - 

oXt  « - 

- .  I 


ClDSING*  /  ax,  ‘POINT  OK-  FLIGHT* 


f  NO 


FIGURE  D-5D.  PROGRAM  LISTING  -  MODNEF  CONTROL  ROUTINE 


D-20 


W.arg* 

Contours 


FIGURE  D 


6B.  MERGENEF  ROUT 


D  -  2 


D  -22 


GU 


or>r»  ooorvooonnnooooonooooonoooooooooo 


MfRGENEE 


DIMENSION  Tl(  10 » , T2I  10),TP(22), TM.(  22 )  , T*S (  21 1 
01  MENS  I DN  CMU0),C1(10),C2!  10)  ,CP(22)  ,CN(22)  ,CB<21) 
DIMENSION  RX1CLSI 10) ,HX2CLS( 13) 

INTEGER  XiCLS(lO),  X2CL  SI 10 ) .  I RR AY  1 200 » 10 ) 

INTEGER  DXM,  CXI,  0X2,  XO,  YO,  XEND 
INTEGER  XBE3 ,  X9EG1 ,  XBEG2 
DIMENSION  LBLMI 20)  ,  LBLK20),  LBL2I20) 

DIMENSION  MYP(IO),  MYM(IO),  MY  A I  1C),  MY9(10) 


NEE  CONTOUR  MERGING  ROUTINE 


PURPOSE.  TO  ACCEPT  CONTOURS  FROM  TWO  PARALLEL  RUNWAYS 
AND  GENERATE  Tut  MERGED  CONTOURS. 


INPUT. 

THE  DATA  FOP  THE  MERGER 

1.  LA  TEL  F CP  THE  MEPCfP 

2.  X-START,  X-INTERVAL,  AND  ThE  CONTOUR  VALUES  FOR  THE  MERGER 

3.  X-END,  X-OFFSET,  v-OFFSET  OF  THE  2ND  =>UNWAY  FFOM  THE  1ST 

THE  DATA  FCR  THE  FIRST  RUNWAY 

1.  LABEL  FOR  THE  FIRST  RUNWAY 

2.  X- START,  X-INTERVAL,  CONTOUR  VALUES 

3.  X-VALUE,  Y-VALUES  FOR  EACH  CONTOUR 

4.  THE  CONTOURS  ARt  TERMINATED  WITH  A  BLANK.  CARO 

5.  THE  SET  OF  DATA  FOR  THE  FfPST  RUNwAY  IS  TERMINATED 
WITH  THE  X-VALUES  FOR  THE  CLOSURE  OF  EACH  CONTOUR 

THE  DATA  FOR  THE  SECOND  RUNWAY  IS  SIMILAR  TD  THAT  OF  THE  FIRST 


READ  THE  DATA  FOR  THE  MfcRGCk 

300  READ  (1,10031  LBLM 

READ  ( 1 , 1003 i  XBEG,  DXM,  (CM|I),  1=1,10) 

IF  ( OXM )  205,205,301 

301  CONTINUE 

THE  PROGRAM  STOPS  WHEN  THE  X-INTERVAL  FOR  MERGER  IS  3 

DO  l  I  =  1 , 1 C 
IF  ( CM ( I))  2,2,1 
1  CONTINUE 
GOTO  3 
21=1-1 
3  NCM  =  I 

READ  (1,1(0!)  XEND,  XO,  YC 
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C  READ  THE  DATA  RIP  THE  FIRST  RUNWAY 

READ  (1,1003)  L8L1 

REAP  (1,100P)  XBEG1,  0X1,  (Cl(l),  1=1,10) 
00  4  I  =  1,10 
IE  (CHI))  5,?, 4 

4  CONTINUE 
GOTO  6 

5  1=1-1 

6  NCI  =  I 
J  =  l 

7  READ  (1,1002)  IX,  (IPFAY(J,I),  1=1, NCI) 

IF  ( 1  PR  AY ( J ,  1 ) )  302,3C2,P 

30?  IF  (IX)  9,9,8 
8  J  =  J  *  1 
GOTO  7 

o  NXl  =  J  -  l 

READ  (  1,1002)  ( XI  CL  SI  I ) ,  I  =  1 , NC 1  I 
DO  7 CO  NR  = 1 , NC 1 
700  PXICLS(NR)  =  XICLS(NR) 

C 

C  READ  THE  DATA  FJF  THE  SECOND  RUNWAY 

READ  (  1 ,  1CC3  )  L9L2 

REAP  (1,1C00)  XHEG2,  0X2,  (C2(I),  1=1,10) 
DO  10  I  =  1,10 
IF  ( C  2 ( I  ))  11,11,10 

10  CONTINUE 
GOTO  12 

11  1=1-1 

12  NC2  =  I 

J  =  NXl  ♦  1 

13  FEAO  (1,1002)  IX,  ( I RF  AY ( J , I ) ,  1  =  1, NC2) 

IF  ( I  PR  AY ( J ,  1 ) )  303,303,14 

303  IF  (IX)  15,19,14 

14  J  =  J  ♦  1 
GOTO  13 

15  NX2  =  J  -  NX!  -  l 

READ  (  1,1002)  ( X2CL  S (  I  I ,  1  =  1, MC2) 

DO  701  NP=  1 » NC 2 
701  R  X2CL  S ( NF )  =  X2CLS(NR) 


WP  I  te 

(3,3010 

WP  I  TE 

THE  INPUT 

DATA 

FOR  IDENTIFICATION 

PURPOSES 

WP  I  TE 

( 3,3011 ) 

IBLM 

WP  I  TE 

( 3,3012 i 

X8EG, 

CXM,  ( C  Y  (  I  )  ,  1=1, 

NCM) 

WRITE 

( 3,3012) 

XEND, 

XO,  YO 

WP  I  TE 

(3,3020) 

WP  I  TE 

(3,301 1 1 

LBL  1 

WP  I  TE 

( 3,30121 

XBEGl 

,  DXl,  (Cl( I),  1=1 

,n;  i ) 

WPI  TE 

(3,3020 1 

WRI  TE 

(3,3011) 

LBL  2 

WRITE 

( 3,3012) 

XBEG2 

,  0X2,  <C2( I ) ,  1=1 

,  NC  2 ) 

LINES 

=  60 

FIND  THE  FIRST  X-VALUE  WHEFE  BOTH  SETS  OF  CONTOURS  ARE  DEFINED 
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XT  =  XBEG1 

IF  (XT  -  XBF  '»2  -  XOI  H,17,17 

16  XT  =  X3EC2  ♦  XO 

17  X  =  X8EG 
R  X  =  X 

I F ( X  »  10,702,  IP 
702  R X  =  1. 

18  IF  (X  -  XT)  10,20,20 
10  X  =  X  +  DX.m 

RX  =  X 
GOTO  18 

CHECK  FOR  CLOSURE  OF  THE  CONTOUR 

20  CONTINUE 

J  =  IX  -  XBFGl)/nxi  +  1 
XP  =  (J  -  1)  *  DX1  ♦  X8EG1 
DO  26  I  1 , NC l 
IF  (IRRAYIJ, I ))  22,22,28 

28  CONTINUE 

IF  (XlCLSim  25,25,21 

21  IF(X-  XICLSI  I  )  )  23,22,22 

22  T1 (I)  =  0. 

GOTO  26 

23  IF  (XP  +  0X1  -  XICLSI  I )  )  25,24,24 

24  TKIl  =  ((  XICLSI  I  )-X)/(XlCLSm  -  X  P  I  )  *1  RR  AY  ( J  ,  I  ) 

GOTO  26 

25  T1II)  =  (  (X-XP)/0X1 )*IIPPAY(  J*i, I )-lRRAY( J, I ) )  +  IRRAY( J,  I) 

26  CGNTINUE 

INTERPOLATE  TO  F I  NO  THE  Y-VALUES  FOR  THE  TWO  CONTOURS  SETS  AT  X 

J  =  IX  -  XO  -  X3EG2J/0X2  *  1 
XP  =  (J-l)  *  0X2  +  XO  *■  XH  LG? 

K  =  J  ♦  MX1 
DC  36  I = l , NC  2 
IF  I  I RR AY  IK, I ))  32, 32, 38 
38  CONTINUE 

IF  ( X2C  L  S ( I ) )  35, '35, 31 

31  IE  ( X  -  X2CLSI  I)  -  XO'  33,32,3? 

32  T2(  I  )  =  0. 

GOTO  36 

33  IF  (XP  +  0X2  -  X2CLSII)  -  XO )  35,34,34 

34  T2II)  =  (IX2CLSI  I) — XI/CXPCLSC I  )-XP)  )*IRRAY(K,I  ) 

GOTO  36 

35  T2 (  I)  =  ( ( X  -  XP)/DX25  * 

C  I  I  KRAY ( K  + 1 , I  )  -  I RR AY ( K , I  ) )  ♦  I RR  AY ( K , I  I 

36  CONTINUE 

Tl,  T 2  ATE  THE  SET  OF  INTERPOLATED  Y-VALUES  FOR  THE  CONTOUR 
VALUES  OF  THE  FIRST  AND  SECOND  RUNWAYS  AT  X 

INTERPOLATE  UNDER  FLIGHT  PATH 
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C  CXI  IS  THE  CONTOUR  VALUE  ONI  FIRST  RUNWAY  IUNDER  PATH) 

C 

350  NJ  =  NCI-1 

Dfl  351  1*1,  NJ 

IF  (XICLS1L1 )  351 ,351,352 

351  CONTINUE 

C  X 1  =  (0.  -  T l (NC l-l ) )/ ( T1 (NC L )  -  Tl(NCl-l))* 

C  (C11NC1)  -  CKNU-m  ♦  C1INC1-1) 

GOTO  360 

352  JS  =  L+l 

DP  353  J  =  J  S  *  NC  l 
I F ( x 1CL  S ( J ) ) 355,  355,366 

356  IT  (X  -  XICIS(J))  353,354,354 

353  CONTI  Nil  F 

1356  CXI  =  Tl<J-2)/ITltJ-2)-TUJ-in*(LHJ-l)-Cl(J-?J  )  ♦  CLIJ-2) 

GOTO  360 

355  I F ( J-2 )  1356,1355,1356 

1355  WF  I  T E ( 3 , 304 8 )  X , C l ( J ) 

WRITE13* 30 11)  mil 
LINES  =  LINES  ♦  6 
GOTO  390 

354  CXI  =  ALOGIPXICLS  IJ-U/RX)  /  UPC,  ( R  Xl  CL  S(  J  - 1 ) /RX  ICl  S  (  J  ))  * 

C  (Cl ( J )  -  Cl(J-l))  ♦  C1(J*1I 

CX2  IS  THE  CONTOUR  VALUF  ON  SCCOND  RUNWAY 

360  NJ  *  NC  2-1 

00  361  L  =  1 , N J 
IF  1X2CISIL))  351,361,362 

361  CONTINUE 

CX2=  (C.  -  T  2  I  NC  2— 1 1 ) / 1 T  2 1 NC  2 )  -  T2INC2-1))* 

C  IC7IMC2)  -  C2(NC.2-1>)  ♦  C2HC2-1) 

GOTO  357 

362  JS  =  L«-l 

DO  363  J= J  S  »  NC  7. 

I F | X2CL  S ( J  » )  365,365,366 
366  IF  (X  -  X2CL  S ( J )  I  363,364,364 

363  CONTINUE 

1366  CX2  =  T2I J-2) /( T2I J-2)-T2( J-l) ) *IC2( J-l )-C2( J-2))6C2( J-2) 

GOTO  357 

365  IFIJ-2)  1 36 5, 1345, 1366 
1365  WRITE(3,304P)  X,  C2IJ) 

KR I TF( 3, 301 1 )  L8L2 
LINES  =  UNIS  ♦  6 
GOTO  ?.«'« 

364  CX2  =  AL0G1RX7CLSI  J-D/AHSIP.X-XO)  )  /  AL  flG  I RX2CL  S  ( J- 1  I  /*  X2CL  S I  J  l  )  * 
1  (C2 < J)-C2C J-l )  )*C2< J-l) 

357  Y  =  0. 

K  =  4 
GOTO  520 

359  CCX1  =  lC.*ALf  G10I1C.**1CX1/10.  )  1C,**1  CC2/10.  )  ) 

Y  =  YO 
K  =  4 
GOTO  5 10 

36H  F  C  X  2  =  16.*ALOG)OUO.**(CX?/lO. )  6  10. CC1/19. I) 
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GENERATE  THE  CONTOUR  VALUE  VS.  V  FUNCTION  AT  THE  GIVEN  « 

ALONG  THE  LINE  TO  THE  RIGHT  OF  THE  TWO  CONTOURS  AND  TO  THE  LEFT. 


CP, TP  ARE  THF  ASSOCIATED  LISTS  OF  CONTOUR  AND  Y-VALUtS  TO  THE  RIGHT  OE 
THE  TWO  RUNWAYS. 


CNtTM  ARE  THE  ASSOCIATED  LISTS  UR  CONTOUR  AND  Y-VALUES  TO  THE  LEFT  OF 
THF  TWO  RUNWAYS. 


THE  CASE  WHERE  THE  FIRST  RUNWAY  IS  THE  RIGHTMOST 
OP  WHFPE  TWO  PATHS  ARE  SUPERPOSED 

IT  <YO)  401,401*411 

401  NP  =  0 

DO  404  I  =  I, NCI 

Y  =  Tl(I) 

IF  ITKIII  404,404,40  2 

402  K  =  I 
GOTO  520 

403  NP  =  NP  +  I 

CP(NP)  =  10.*ALOG10(10.**(C1 (  11/10.  )  ♦  13.  **(CC2/ ID.  I  ) 

TP(NP)  =  Y 

404  CONTINUE 

D°  1410  I  =  1,  NC  2 

Y  =  T2( I )  +  YC 

I F  t  Y I  1410,1410,1405 

1405  i<  =  5 
GOTO  510 

1406  NP  =  NP  +1 

f  P(NPJ=10.*ALf  G10(  1P.*»(C2(  II/10.I+10.**(CC1/10.  »  I 
TP (NP  I  =  Y 

1410  CONTINUE 
NO  =  1 
GOTO  600 

1411  NM  =  0 

DO  403  1  =  1, Nr  2 

Y  =  -T2 ( I )  ♦  YO 

IF  (T2(II)  407,407,405 

405  K  =  1 
GOTO  510 

406  NM  =  NM  1 

CN(NM>  =  1C.*AUG10UO.**(C2<  I  I  /  10 .  )  ♦  1  0.  **  (  CC  l  /  1 0  .  U 
TM(NM)  =  -T2  (  I) 

407  CONTINUE 

DO  1450  1=1, NCI 

Y  =  -T1III 

IF(Y-YO)  1415, 14S0, 145C 

1415  K  =  5 
GOTO  520 

1416  NM  =  NM  ♦  1 
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■  1 

cmnmi  =  ic.*^inGioi  if  .♦*(Ci  i  n/io.  i 

♦  10.**«CC2/10.  1)  B 

1 

TH(NM)  =  Y  -  YC 

1  . 

1450 

CONTINUE 

MO  =  1 

GOTO  650 

I 

1451 

NP  =  NP  ♦  1 
fPCNP)  =  crxi 

IP(MP)  =  1. 

1 

NM  -  MM  ♦  1 

I 

CN1NM)  =  CCX2 

] 

TM(NM)  =  -1. 

f 

r 

YP  =  C. 

1 

YM  =  Y1 

GOTO  420 

C 

I 

C 

The  CASE  PHFPE  The  SECOND  PUNuA*  is 

TKf  RIGHTMOST  » 

c 

j 

411 

NP  =  0 

1 

1 

DC  414  1  =  1 ,NC2 

Y  =  T2( 11  ♦  YC 

IF  (T?m»  414,414,41? 

412 

K  =  2 

i 

1 

GOTO  510 

413 

MP  =  MP  ♦  1 

cpimpj  -  io.*Ainc,io<  ic.**<c?m/io.) 

♦  1 D. **( CC 1/10. 1 1 

TP1NP1  =  T? ( I 1 

414 

CONTINUE 

CC  1 46.0  I  =  l.NCl 

y  =  Tim 

1KY-Y0I  146C, 1466,1425 

1425 

K=f- 

GMO  520 

1426 

MP  =  MP  ♦  1 

CP1NPJ  =  lO.MALOGlOf  ic-.**icim/io.) 
TPINP)  =  v  -  Y0 

♦  n.**(CC2/10) »  1 

1460 

CCMTIMUE 

NO  =  2 

GOTO  600 

1461 

MM  =  0 

DC  417  1  =  1, NCI 

Y  =  -TKI) 

IF  (TUI))  417,417,415 

415 

K  =  2 

GCTO  520 

*  T 

416 

NM  =  MM  ♦  1 

CN(NM)  =  10.*ALnG10(ir-.**(Cl(l)/10.) 
TMCNM)  =  Y 

♦  13.**CCC 2/10.1 1 

417 

CONTINUE 

Or  1470  I  =  1  ,  NC 2 

Y  =  -T2  (  II  ♦  YC- 

1 

i 

I F ( Y  )  1435,1470,1470 

1435 

K  =6 

GrTO  510 

1 

1436 

NM  =  NM  *•  1 

I 

CM  NM )  =  10.  *  ALOGlO(lC.**(C2m /10 

)  *  l 0.**(CC1/10. ) 1 
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TM(NM)  =  Y 

1470  CONTINUE 
NC  =  2 
GOTO  650 

1471  NP  =  NP  ♦  1 
CP1NPJ  =  CCX2 
TP1NP)  =  1. 

NH  ■  NM  ♦  1 
CNCNM)  *  CCXi 
TM(NM)  =  -l. 

YP  =  YO 

YM  =  0. 

C 

C  CB,  TB  APE  THE  ASSOCIATEO  l ISIS  OF  CONTOJR  VALUES  ANJ  Y-VALUES 

C  FOR  THE  REGION  PFTwCEN  THE  TWO  F UNwAYS  AT  THE  GIVEN  X-DISTANCE 

C 

420  Y  =  0. 

OY  =  YO/20. 

DO  423  1=1,21 

K  =  3 
GOTO  510 

421  K  =  3 
GOTO  520 

422  CB( 1 1  =  10.*AL0G10( 1C.**(CC1/1G. »  ♦  10.**l CC2/10. 11 
TBU)  *  Y 
Y  =  Y  ♦  OY 

423  CONTINUE 

INTERPOLATE  BETWEEN  CONTOURS  TO  GET  THE  Y-VALUES  F3*  THE  MERGED  CONTOURS. 

MYA,  MY B  APE  TWO  POSSIBLE  Y-VALUES  F QP.  EACH  CONTOUR  BETWEEN 

43 C  DC  442  1=1 , NCM 

MYA ( I )  =  0 

MYBU)  =  C 

IF  (YO)  431 , 442  , 431 
431  DO  441  J= 1 , 20 

IF  (CM(I)  -  CB(JI)  433,434,434 

433  IF  (CMI  I »  -  CPU  +  m  441,435,435 

434  IF  ( CM ( I )  -  CPU+ll )  ‘435,435,441 

435  IE  (CB(JM)  -  Cfl  ( J ) )  437,436,437 

436  TY  =  <TB(J+1)+TB(JI  >/2. 

GOTO  438 

437  TY  =  (CM(I)  -  C 13  ( J )  )  /  (CB(JH)  -  CB(JM  * 

C  (TBU+1)  -  TP(  JJ  )  +  TR(JI 

438  IF  ( MYB ( I ) )  440,439,440 

439  MYR( I )  =  TY 
GOTO  441 

440  MYA (1)  =  Tv 

441  CONTINUE 

442  CONTINUE 

MYP  IS  THE  LIST  OF  Y-VALUES  FOR  THE  CONTU JR$  TO  THE  RIGHT  OF  THE  RUNWAYS. 
DO  455  1=1, NCM 
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vfcOGFNEF 


451 

452 

453 

454 
459 

45o 
45  7 


45  6 

C 


455 


MYP< I  )  =  ? 

J  =  2 
IF  INF  -n 
IF  <CM(  l  »  ■ 
no  453  J  = 

IF  ICM< I  I  - 
CONTINUE 
GOTO  455 
IF  (J  -  NPJ 
TT  »  TPIJ-1) 

IF  t  TT  I  467  ,456  »457 
TT  =  TP( J)/APS( TP(J)  ) 

MYpin  =  ic.**((Oti)  -  cpu-i n/icpui 

ALOGl^IAHSI TP( J  )/TT > ))*TT  ♦  YP 
GOTO  4*5 


453,455,451 

CPI l I »  45*, 452, 45? 

1  ,NP 

(PIJ1)  454,453,453 


45= ,45*, 45* 


CPI J-  111* 


VYPUJ  =  (C.KII)  -  CF  (J-l  )  l/ICPCJ  >  -  CPCJ-IM* 
ITP(J)  -  TPJJ-ll)  ♦  TP  I J- 1 1 ♦  YP 
CONTINUE 


461 

462 

463 


C 

C 

C 

c 

c 

c 


4 1  4 
46  9 

466 

467 


468 

C 

46  5 


211 


C 

c 


MYM  IS  THC  LIST  OF  Y— VALUE S  FOR  THF  CONTOURS  TO  THE  LEFT  OF  T HF  RUNWAYS. 


DO  465  I  =  l.OCO 


MYM ( 11=0 

J  =  2 

IF  (  NM  -1)  465 , 465 , 46  1 

IF  (  CM!  ■  )  -  C  L:  1 1  1  )  46P,462,462 

DC  463  J  =  1,NM 

IF  (CMC  I  )  -  CN(Jl)  464,463,463 

CONTINUE 

GOTO  464 

IF  Ij  —  N'J  ‘  /  , 4  6  o •  -  6 

TT  =  TM(J-l) 

IF  (TT1  467,466,467 
TT  =  TMlJJ/ABSlTMIJ)  ) 

”Y.M(I)  =  IP.**UCH(N  -  CNIJ-l)  I/ICNIJ) 
ALOGIOIABSI  TM(  J)/TTMI*TT  ♦  YM 


CNIJ-l)  )* 


COTO  465 

NYMII)  =  ICM(I)  -  CNIJ-l )) /(CNIJ )  -  CNIJ-l))* 
(TM(J)  -  T M  C  J— 1)1  +  TM  (  J- II-*-  YM 


CONTINUE 


WRITE  THE  OUTPUT 
LINES  =  LINES  ♦  6 

54  LINES  APE  PUT  ON  FACH  dAGE  OF  OUTPUT 

IF  (LINES  -  48)  212,212,211 

LINES  =  0 
WRITE  (3,30211 
WRITE  ( 3,3020) 

WRITE  (3,3022)  (CM(I),  1  =  1,  NC*) 

WRITE  (3,3020) 

THE  CONTOUR  VALUES  TER  THF  M E R G F R  ARE  WRITTFN  AT  THE  TOP 
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THt  X-VALUF  FCR  THE  MERGED  CCNTOUR  Y-VAIUES  fS  WRITTEN 
THE  V-VALUES  ARE  WRITTEN 


212  WRITE  13*  30?  3 )  X 

WRITE  (3,10241  ( MY? (II*  1=1, NCM) 

WRITE  ( 3 ,30? 5 ) 

DO  2?1  1*1, NCM 

K  =  1 

M  =  MYA  (  I  ) 

IF  (M)  240,230,240 

221  CONTINUE 
WRITE  (3,302ft) 

DP  222  1*1, FT v 

K  =2 
H  =  MYB(  I  I 
IF  (MJ  240,730,240 

222  CONTINUE 

WRITE  (3,30291  (MYM(I),  1=1, NCR) 

WRITE  (3,3020) 

THE  X-VAIUE  IS  INCREMENTED 
AND  END  CHECKS  ARE  P£KF  DP MED 

399  X  =  X  +  DXM 

IF  (X  -  XEND)  201,201,300 

201  XP  =  DX 1  * ( NX  1 - 1 )  *  XBFG1 
IF(X-XP)  202,300,300 

202  XP  =  DX2*  ( NX2  -  1)  ♦  XO  ♦  XBEG2 
IF  (  X  -  XP)  203,300,300 

2C3  RX  =  X 
GOTO  20 

IF  THIS  MFRGEC  CONTOUR  IS  OUNE,  A  NEW  SET  OF  DATA  IS  READ  IN. 
205  STOP 

230  GOTO  1231, 232, 233, 234, 235, 236, 237), I 

231  WRITE  (3,3031) 

GOTO  238 

232  WRITE  (3,3032)  • 

GOTO  238 

233  WRITE  (3,3033) 

GOTO  238 

234  WRITE  (3,3034) 

GOTO  238 

235  WRITE  (3,3035) 

GOTO  238 

236  WRITE  (3,3036) 

GOTO  238 

237  WRITE  (3,3037) 

238  GOTO  (221, 222), K 

240  GOTO  (241,242,243,244,245,246,247) , I 

241  WRITE  (3,3041)  M 
GOTO  248 

242  WRITE  (3,3042)  M 
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GOT 0  ?4R 

243  WPITE  (3,30^3)  M 
GCT'1  248 

244  fcPITf  (3,30441  M 
GOTO  248 

245  WPITE  (  3,3045)  M 
GOTO  24B 

246  y  F  I  Tfc  (3,30^6)  'I 
GOTO  24 B 

247  WP1TF  (3,3047)  V 
24P  GOTO  (  221,22/  )  ,K 


TO  Gt'T  THE  CONTOUR  VALl  F  FOR  T  HC  FIRST  RUNWAY  AT  .4  GIVEN  Y-VAUJE 

510  CCl  =  0. 

T  Y  =  Y 

IF  (TV)  511,512,51 2 

511  TY  =  -TY 

512  IF  (Tl(l)l  530,530,513 

513  J  =  1 

IF  ( TY  -  T 1 ( 1 ) )  514,514,517 

514  NJ  =  NCI  -  1 

DO  516  J  =  1  ,  N J 

IF  (TY  TIIJI)  515,515,516 

515  IF  (TY  -  TllJ+lll  516,517,517 

516  CONTINUE 

CCl  =  ( TY/T1 (NCI ) I*CC1 (NCI )  -  CXI)  ♦  CXI 
GOTO  530 

517  IF(TY)  1517,151b, 151? 

1516  TY  =  1. 

1517  IF  (Tl(Jtl)i  518,518,515 

518  CCl  =  ALOG ( T 1  (  J  )  /  TY )  /  A  LOG  (  T  1  (  J  ?  >  *  (  C  X  1  -  ClIJI)  ♦  CKJI 

GCT 0  530 

519  CCl  =  Atn:;(Tl  (J)/TY)/AlCC(Tl(J)/ri(j4i*)*(Cl(J  +  n-Cl(Jl  I+CIIJ) 
530  GOTO ( 406 ,413,421 ,368, 1406,1436)  ,K 

TO  GFT  THE  CONTOUR  VALUE  FOR  THE  SECOND  RUNWAY  AT  A  GIVEN  Y 


520 

CC2 

=  0. 

TY 

=  Y  - 

Y  0 

I  F 

(  TY) 

82  1  ,522, 

52? 

521 

TY 

=  -TY 

522 

IF 

( T? ( 1 ) )  540,540,523 

523 

J 

=  I 

IF 

(  TV  - 

T2(  1  )  ) 

524, 5?4,5?7 

524 

NJ 

=  NC2 

-  1 

DO 

526  J 

=  1 ,  N  J 

IF 

(TY  - 

T2( J  )  ) 

525,525,526 

525 

IF 

<TY  - 

T2  (  J  +  1  1) 

526, 527, 5. 

526 

( ONTIMUE 

C  C  2  =  (TY/T2(NC2) )*(C2(NC2)  -  C  X  2 ) 
GOTO  540 

520  IF (TY)  1527, 152b, 1527 
1526  TY  =  1. 


CX2 
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1527 

528 

529 
540 
6C0 


601 


602 

603 

604 

605 

606 

60  8 
609 

607 


650 


651 


652 

653 


654 


If  (T2(JU))  528,528,529 

CC?  =  AL0G(T2( J)/TV)  /  AL0GU2U)  >*(CX2  -  C2U))  ♦  C  2<  J  I 
GOTO  540 

CC2  =  AL0GiT2(J)/TY)/AL0G(T2U)/T2(J*ll)*(C2UM)“C2U))*C2U) 

GOTO  (403,416,422,358,1416, 1426) ,K 
NC  =  0 

NM1  =  NP  -  1 

00  6C7  I  =  1 , MM  1 

IP1  =  1  *•  l 

00  607  J=I?1,MP 

lf(CPUI  -  CP  (  I )  )  60 1  , 60  2 , 6  OP 

CT  =  CP (  I  ) 

CPU  >  -  CPU) 

CPIJ)  =  CT 
CT  =  TP (  I  ) 

TP ( I )  =  TP ( J ) 

T  P ( J )  =  CT 
GOTO  60? 


I F ( A B S ( T P < I ) )  -  AHS ( TP( J | ) )  603,604,606 

C P ( I  )  =  CP( I )  +0.01 

C P ( J )  =  CP(J)  -  0.01 

GCTO  601 

DO  605  K=J,NM1 

C  P ( K )  =  C  P  ( K  + 1 ) 

TP  I K )  =  TP(KU  ) 

NC  =  NC  ♦  1 
GCTO  607 

CPU)  =  CPU)  -  0.01 
CP(J>  =  CPU)  ♦  0.01 
GOTO  607 
IF(ABSITPU)) 

TP(I )  =  TP(I ) 

TP ( J )  =  TP( J )  -  1. 

CONTINUE 
NP  =  NP  -  NC 
GOTO  (1411,1461),  NO 
NC  =  0 
NM1  =  MM  -  1 
DO  657  1=1, N«1 
I P 1  =  l  +  1 
DD  657  J  = I P 1 ,NM 
I  FICNIJ  I  -  CNH1I  651,652,658 
CT  =  CM (  i  ) 


ABS(TP(J)))  607,609,607 

1 . 


CNU  )  =  CN(J  ) 

C N ( J  )  =  CT 
CT  =  TMU) 

TM(I)  =  T  M ( J ) 

TMU)  =  CT 
GOTO  657 

I  F ( A  B  S ( T  M ( I  ) )  -  A  0  S I TM ( J ) ) )  653,654,656 
CNU  )  =  CNU  )  ♦  0.01 
CN(J)  =  CNU )  -  0.01 
GOTO  651 


DO  655  K= J , NM 1 
C  N  (  K  )  =  CN ( K ♦ 1 ) 
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655  TM(  K )  t  T M ( K ♦ 1 | 

NC  =  NC  +  1 
GOTH  657 

656  CN(  I  )  =  C Ni (  1  )  -  O.Ol 
CMJJ  =  CN1J)  ♦  0.01 
GOTO  657 

65«  IHABSITMU))  -  ABS(TM(J)))  657,650,657 
650  TM(  I  )  =  TM  I  )  ♦  i. 

T  M  (  J  >  =  r  M  t  J  J  -  1. 

657  CONTINUE 

NM  »  NM  -  NC 
GOTO  (1451,1471)  ,  NO 
C 

1000  FORMAT  (2(7,  10F4.0) 

1001  FORMAT  (317) 

1002  FORMAT  (1117) 

1003  FORMAT  (?CA4) 

3000  FORMAT  (  1017) 

3001  FORMAT  (•  X  VALUE  =  1F8.0) 

3002  FORMAT  ( •  TS'  ,  10F7. 1  ) 

3003  FORMAT  (•  CS',  1 0  F  7 . 1  ) 

3004  FORMAT  ( ■  MS' ,  418) 

3005  FORMAT  (‘  517) 

3006  FORMAT  (•  ',  17) 

3010  FORMAT  ('!»> 

3011  FORMAT  (•  20A4) 

3012  FORMAT  (•  ',  4X,  ' S  T  AR  T=  ',  17,  4X,  'INTERVAL*  17, 

C  4X,  'CONTOURS  ',  1GF4.0) 

3013  FORMAT  (•  ',  4 X ,  'STOP  =  17,  4X, 

C  'X,Y  OFFSET  OF  2ND  RUNWAY*  ',  17,  «,  17) 

3020  FORMAT  (•  •) 

3021  FORMAT  (  •  1  •  ,  16X, 

C  'DISTANCE  TO  MERGED  CONTOUR  (IN  FEET)  FROM  FIRST  RUNWAY') 

3022  FORMAT  (•  ',  'CONTOUR  VALUES  ',  10F7.0) 

3023  FORMAT  ('  *x  VALUF*  ',  l  F  F .  0 ) 

3024  FORMAT  ('  ',  6X,  'TP  RIGHT  ',  1017) 

3025  FORMAT  («  6X,  'BETWEEN',  3X) 

3028  FORMAT  («  6  X ,  'RUNWAYS',  3Xf 

302°  FORMAT  (•  ',  6X,  'TO  LEFT  ',  1017) 

3031  FORMAT  ('s',  23X) 

3032  FORMAT  (•<■',  30X) 

3033  FORMAT  ('♦',  37X) 

3034  FORMAT  (•«■',  ^4X) 

3035  FORMAT  (•  +  •,  51X) 

3036  FORMAT  (•'<•,  FAX) 

3037  FORMAT  (  '♦'  ,  65X ) 

3041  FORMAT  (•♦',  1 6  X ,  17) 

3042  FORMAT  ('«•',  23X,  17) 

3043  FORMAT  Cm',  30X,  17) 

3044  FORMAT  (•♦',  37X,  17) 

3045  FORMAT  (•♦',  44X,  17) 

3046  FORMAT  (•♦',  5 1 X ,  17) 

3047  FORMAT  ('♦',  58 X ,  (7) 

3048  f  (J  RM  A  T  (  '  X  VA  L  IJF  =  ',  IFB.O,'****  CLOSURE  3F  ',F7.0, 

1'  CONTOUR  NFFUED.') 

END 
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